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Abstract 

We consider the time-dependent Schrodinger equation on a Rie- 
mannian manifold A with a potential that localizes a certain subspace 
of states close to a fixed submanifold C. When we scale the potential 
in the directions normal to C by a parameter e <C 1, the solutions con- 
centrate in an e-neighborhood of C. This situation occurs for example 
in quantum wave guides and for the motion of nuclei in electronic 
potential surfaces in quantum molecular dynamics. We derive an ef- 
fective Schrodinger equation on the submanifold C and show that its 
solutions, suitably lifted to A, approximate the solutions of the origi- 
nal equation on A up to errors of order e^\t\ at time t. Furthermore, 
we prove that the eigenvalues of the corresponding effective Hamil- 
tonian below a certain energy coincide up to errors of order with 
those of the full Hamiltonian under reasonable conditions. 

Our results hold in the situation where tangential and normal ener- 
gies are of the same order, and where exchange between these energies 
occurs. In earlier results tangential energies were assumed to be small 
compared to normal energies, and rather restrictive assumptions were 
needed, to ensure that the separation of energies is maintained during 
the time evolution. Most importantly, we can allow for constraining 
potentials that change their shape along the submanifold, which is the 
typical situation in the applications mentioned above. 

Since we consider a very general situation, our effective Hamilto- 
nian contains many non-trivial terms of different origin. In particular, 
the geometry of the normal bundle of C and a generalized Berry con- 
nection on an eigenspace bundle over C play a crucial role. In order 
to explain the meaning and the relevance of some of the terms in the 
effective Hamiltonian, we analyze in some detail the application to 
quantum wave guides, where C is a curve in ^ = M^. This allows us 
to generalize two recent results on spectra of such wave guides. 
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1 Introduction 

Although the mathematical structure of the linear Schrodinger equation 

idt^ = -A^ + V^=:Hilj, ^P\t=o^L^iAdT) (1) 

is quite simple, in many cases the high dimension of the underlying configu- 
ration space A makes even a numerical solution impossible. Therefore it is 
important to identify situations where the dimension can be reduced by ap- 
proximating the solutions of the original equation ([1]) on the high dimensional 
configuration space A by solutions of an effective equation 

idt(P = H,^cP, (t)\t=^eL\C,dii) (2) 

on a lower dimensional configuration space C. 
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The physically most straight forward situation where such a dimensional 
reduction is possible are constrained mechanical systems. In these systems 
strong forces effectively constrain the system to remain in the vicinity of a 
submanifold C of the configuration space A. 

For classical Hamiltonian systems on a Riemannian manifold {A, G) there is 
a straight forward mathematical reduction procedure. One just restricts the 
Hamilton function to T*C by embedding T*C into T*A via the metric G and 
then studies the induced dynamics on T*C. For quantum systems Dirac [12] 
proposed to quantize the restricted classical Hamiltonian system on the sub- 
manifold following an 'intrinsic' quantization procedure. However, for curved 
submanifolds C there is no unique quantization procedure. One natural guess 
would be an effective Hamiltonian H^^ in ([2]) of the form 

i/eff = -Ac + V\c , (3) 

where Ac is the Laplace-Beltrami operator on C with respect to the induced 
metric and V\c is the restriction of the potential V : A-^M.io C. 
However, to justify or invalidate the above procedures from first principles, 
one needs to model the constraining forces within the dynamics ([T]) on the full 
space A. This is done by adding a localizing part to the potential V . Then 
one analyzes the behavior of solutions of ([T]) in the asymptotic limit where 
the constraining forces become very strong and tries to extract a limiting 
equation on C. This limit of strong confining forces has been studied in 
classical mechanics and in quantum mechanics many times in the literature. 
The classical case was first investigated by Rubin and Ungar [39l|, who found 
that in the limiting dynamics an extra potential appears that accounts for the 
energy contained in the normal oscillations. Today there is a wide literature 
on the subject. We mention the monograph by Bornemann |3] for a result 
based on weak convergence and a survey of older results, as well as the book 
of Hairer, Lubich and Wanner [T9], Section XIV.3, for an approach based on 
classical adiabatic invariants. 

For the quantum mechanical case Marcus [29] and later on Jensen and 
Koppe |23] pointed out that the limiting dynamics depends, in addition, 
also on the embedding of the submanifold C into the ambient space A. In 
the sequel Da Costa [S] deduced a geometrical condition (often called the 
no-twist condition) ensuring that the effective dynamics does not depend on 
the localizing potential. This condition is equivalent to the flatness of the 
normal bundle of C. It fails to hold for a generic submanifold of dimension 
and codimension both strictly greater than one, which is a typical situation 
when applying these ideas to molecular dynamics. 

Thus the hope to obtain a generic 'intrinsic' effective dynamics as in ([3]), 
i.e. a Hamiltonian that depends only on the intrinsic geometry of C and 
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the restriciton of the potential V to C, is unfounded. In both, classical 
and quantum mechanics, the limiting dynamics on the constraint manifold 
depends, in general, on the detailed nature of the constraining forces, on 
the embedding of C into A and on the initial data of ([T]). In this work we 
present and prove a general result concerning the precise form of the limiting 
dynamics ([2]) on C starting from ([T]) on the ambient space A with a strongly 
confining potential V. However, as we explain next, our result generalizes 
existing results in the mathematical and physical literature not only on a 
technical level, but improves the range of applicability in a deeper sense. 
Da Costa's statement (like the more refined results by Froese-Herbst |17] . 
Maraner [27] and Mitchell |32], which we discuss in Subsection II. 2p requires 
that the constraining potential is the same at each point on the submani- 
fold. The reason behind this assumption is that the energy stored in the 
normal modes diverges in the limit of strong confinement. As in the classical 
result by Rubin and Ungar, variations in the constraining potential lead to 
exchange of energy between normal and tangential modes, and thus also the 
energy in the tangential direction grows in the limit of strong confinement. 
However, the problem can be treated with the methods used in P EZl E2] 
only for solutions with bounded kinetic energies in the tangential directions. 
Therefore the transfer of energy between normal and tangential modes was 
excluded in those articles by the assumption that the confining potential has 
the same shape in the normal direction at any point of the submanifold. In 
many important applications this assumption is violated, for example for the 
reaction paths of molecular reactions. The reaction valleys vary in shape 
depending on the configuration of the nuclei. In the same applications also 
the typical normal and tangential energies are of the same order. 
Therefore the most important new aspect of our result is that we allow for 
confining potentials that vary in shape and for solutions with normal and 
tangential energies of the same order. As a consequence, our limiting dy- 
namics on the constraint manifold has a richer structure than earlier results 
and resembles, at leading order, the results from classical mechanics. In 
the limit of small tangential energies we recover the limiting dynamics by 
Mitchell [52] • 

The key observation for our analysis is that the problem is an adiabatic limit 
and has, at least locally, a structure similar to the Born-Oppenheimer ap- 
proximation in molecular dynamics. In particular, we transfer ideas from 
adiabatic perturbation theory, which were developed by Nenciu-Martinez- 
Sordoni and Panati-Spohn-Teufel in [SHlEIlEllESlllllllI, to a non-flat 
geometry. We note that the adiabatic nature of the problem was observed 
many times before in the physics literature, e.g. in the context of adiabatic 
quantum wave guides [7] , but we are not aware of any work considering con- 
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straint manifolds with general geometries in quantum mechanics from this 
point of view. In particular, we believe that our effective equations have 
not been derived or guessed before and are new not only as a mathematical 
but also as a physics result. In the mathematics literature we are aware of 
two predecessor works: in [H] the problem was solved for constraint mani- 
folds C which are dimensional subspaces of M'^+*^, while Dell' Antonio and 
Tenuta [TT] considered the leading order behavior of semiclassical Gaussian 
wave packets for general geometries. 

Another result about submanifolds of any dimension is due to Wittich [15] . 
who considers the heat equation on thin tubes of manifolds. Finally, there 
are related results in the wide literature on thin tubes of quantum graphs. A 
good starting point for it is [TS] by Grieser, where mathematical techniques 
used in this context are reviewed. Both works and the papers cited there, 
properly translated, deal with the case of small tangential energies. 

We now give a non-technical sketch of the structure of our result. The de- 
tailed statements given in Section [2] require some preparation. 
We implement the limit of strong confinement by mapping the problem to 
the normal bundle of C and then scaling one part of the potential in 
the normal direction by e^^ . With decreasing e the normal derivatives of 
the potential and thus the constraining forces increase. In order to obtain 
a non-trivial scaling behavior of the equation, the Laplacian is multiplied 
with a prefactor . The reasoning behind this scaling, which is the same 
as in [TTl [32], is explained in Section With q denoting coordinates on 
C and V denoting normal coordinates our starting equation on NC has, still 
somewhat formally, the form 

\^^^^ = -e^ANci^' + V^{q,e-'^iy)iP' + W{q,iy)iP^ =: H'ljj' (4) 

for V'^|t=o ^ L'^{NC). Here Ajvc is the Laplace-Beltrami operator on NC, 
where the metric on NC is obtained by pulling back the metric on a tubular 
neighborhood of C in ^ to a tubular neighborhood of the zero section in 
NC and then suitably extending it to all of NC. We study the asymptotic 
behavior oi as e goes to zero uniformly for initial data with energies 
of order one. This means that initial data are allowed to oscillate on a 
scale of order e not only in the normal direction, but also in the tangential 
direction, i.e. that tangential kinetic energies are of the same order as the 
normal energies. More precisely, we assume that lleVVolP = (V'o I -^^^hi^o) 
is of order one, in contrast to the earlier works [TTl [32] . where it was assumed 
to be of order e^. Here V'^ is a suitable horizontal derivative to be introduced 
in Definition [H 

Our final result is basically an effective equation of the form ([2]). It is pre- 
sented in two steps. In Section 12.11 it is stated that on certain subspaces of 
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L'^{NC) the unitary group exp{—iH^t) generating solutions of (jl]) is unitarily 
equivalent to an 'effective' unitary group exp(— iif^jjt) associated with ([2]) up 
to errors of order e^\t\ uniformly for bounded initial energies. In Section [2l2] 
we provide the asymptotic expansion of H^^ up to terms of order e"^, i.e. we 
compute HeSfi, Hes,i and ifefr,2 in H^s = Hes,o + £Hes,i + e'^H^s,! + O^e^). 
Furthermore, in Section 12.31 and 12.41 we explain how to obtain quasimodes of 
from the eigenfunctions of -ffeff.o + ^H^^ i + e^ifeff,2 and quasimodes of 
He^fi + eHc^^i+e^HeK,! from the eigenfunctions of if^ and apply our formulas 
to quantum wave guides, i.e. the special case of curves in M^. As corollaries 
we obtain results generalizing in some respects those by Friedlander and 
Solomyak obtained in [16] and by Bouchitte et al. in |6]. In addition, we 
discuss how twisted closed wave guides display phase shifts somewhat similar 
to the Aharanov-Bohm effect but without magnetic fields! 

The crucial step in the proof is the construction of closed infinite dimen- 
sional subspaces of L'^{NC) which are invariant under the dynamics (jlj) up 
to small errors and which can be mapped unitarily to L'^{C), where the effec- 
tive dynamics takes place. To construct these 'almost invariant subspaces', 
we define at each point g G C a Hamiltonian operator H{{q) acting on the 
fibre NgC. If it has a simple eigenvalue band Ef{q) that depends smoothly on 
q and is isolated from the rest of the spectrum for all q, then the correspond- 
ing eigenspaces define a smooth line bundle over C. Its L^-sections define a 
closed subspace of L'^{NC), which after a modification of order e becomes the 
almost invariant subspace associated with the eigenvalue band E{{q). In the 
end, to each isolated eigenvalue band E{{q) there is an associated line bundle 
over C, an associated almost invariant subspace and an associated effective 
Hamiltonian H^^. 

We now come to the form of the effective Hamiltonian associated with a band 
E{{q). For -ffeff,o we obtain, as expected, the Laplace-Beltrami operator of 
the sub manifold as kinetic energy term and the eigenvalue band Ef{q) as an 
effective potential, 

HeSfi = — f^^Ac + Ef. 

We note that {Vc + W)\c is contained in E{. This is the quantum version 
of the result of Rubin and Ungar [3^ for classical mechanics. However, the 
time scale for which the solutions of (jl]) propagate along finite distances are 
times if: of order e^^. On this longer time scale the first order correction 
eHcs^i to the effective Hamiltonian has effects of order one and must be 
included in the effective dynamics. We do not give the details of -ffcff.i here 
and just mention that at next to leading order the kinetic energy term, i.e. 
the Laplace-Beltrami operator, must be modified in two ways. First, the 
metric on C needs to be changed by terms of order e depending on exterior 
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curvature, whenever the center of mass of the normal eigenfunctions does 
not he exactly on the submanifold C. Furthermore, the connection on the 
trivial line bundle over C (where the wave function takes its values) must 
be changed from the trivial one to a non-trivial one, the so-called generalized 
Berry connection. For the variation of the eigenfunctions associated with the 
eigenvalue band Ei{q) along the submanifold induces a non-trivial connection 
on the associated eigenspace bundle. This was already discussed by Mitchell 
in the case that the potential (and thus the eigenfunctions) only twists. 
When Ei is constant as in the earlier works, there is no non-trivial potential 
term up to first order and so the second order corrections in ifeff,2 become 
relevant. They are quite numerous. In addition to terms similar to those 
at first order, we find generalizations of the Born-Huang potential and the 
off-band coupling both known from the Born-Oppcnhcimer setting, and an 
extra potential depending on the inner and the exterior curvature of C, whose 
occurence had originally lead to Marcus' reply to Dirac's proposal. Finally, 
when the ambient space is not flat, there is another extra potential already 
obtained by Mitchell. 

We note that in the earlier works it was assumed that —e'^Ac is of order and 
thus of the same size as the terms in Hcs,2- That is why the extra potential 
depending on curvature appeared at leading order in these works, while it 
appears only in i?eff,2 for us. And this is also the reason that assumptions 
were necessary, assuring that all other terms appearing in our i^eff.o ^-nd i?efr,i 
are of higher order or trivial, including that Ei{q) = Ef is constant. 

We end this section with some more technical comments concerning our result 
and the difficulties encountered in its proof. 

In this work we do not assume the potential to become large away from the 
submanifold. That means we achieve the confinement solely through large 
potential gradients, not through high potential barriers. This leads to several 
additional technical difficulties, not encountered in other rigorous results on 
the topic that mostly consider harmonic constraints. One aspect of this is the 
fact that the normal Hamiltonian Hf{q) has also continuous spectrum. While 
its eigenfunctions defining the adiabatic subspaces decay exponentially, the 
super adiabatic subspaces, which are relevant for our analysis, are slightly 
tilted spectral subspaces with small components in the continuous spectral 
subspace. 

Let us finally mention two technical lemmas, which may both be of indepen- 
dent interest. After extending the pull back metric from a tubular neigh- 
borhood of C in ^ to the whole normal bundle, NC with this metric has 
curvature increasing linearly with the distance to C. As a consequence we 
have to prove weighted elliptic estimates for a manifold of unbounded curva- 
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ture (Lemmas P fc [TII|) . Moreover, since we aim at uniform results, we need 
to introduce energy cutoffs. A result of possibly wider applicability is that 
the smoothing by energy cutoffs preserves polynomial decay (Lemma fT2l) . 

1.1 The model 

Let {A, G) be a Riemannian manifold of dimension d+k {d, A; G N) with asso- 
ciated volume measure dr. Let furthermore C C ^ be a smooth submanifold 
without boundary and of dimension (i/co dimension k, which is equipped with 
the induced metric g = G\c and the associated volume measure dfj,. We will 
call A the ambient manifold and C the constraint manifold. 
On C there is a natural decomposition TA\c = TC x NC of ^'s tangent 
bundle into the tangent and the normal bundle of C. We assume that there 
exists a tubular neighborhood B d AoiC with globally fixed diameter, that 
is there is (5 > such that normal geodesies 7 (i.e. 7(0) e C,7(0) e NC) of 
length 5 do not intersect. We will call a tubular neighborhood of radius r an 
r-tuhe. Furthermore, we assume that 

A and C are of hounded geometry (5) 

(see the appendix for the definition) and that the embedding 

C ^ A has globally hounded derivatives of any order., (6) 

where boundedness is measured by the metric G\ In particular, these as- 
sumptions are satisfied for A = M'^"'"'^ and a smoothly embedded C that is (a 
covering of) a compact manifold or asymptotically flatly embedded, which 
are the cases arising mostly in the applications we are interested in (molecular 
dynamics and quantum waveguides). 

Let be the Laplace-Beltrami operator on A. We want to study the 
Schrodinger equation 

idtiJ = -e'^A^ + , ^\t=o e L\A, dr) , (7) 

under the assumption that the potential localizes at least a certain sub- 
space of states in an e-tube of C with e <^ 6. The localization will be realized 
by simply imposing that the potential is squeezed by in the directions 
normal to the submanifold and not by assuming to become large away 
from C, which makes the proof of localization more difficult. To ensure proper 
scaling behavior, we have multiplied the Laplacian in ([7]) by e^. The physical 
meaning of this is explained at the end of the next subsection. Here we only 
emphasize that an analogous scaling was used implicitly or explicitly in all 
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other previous works on the problem of constraints in quantum mechanics. 
The crucial difference in our work is, as explained before, that we allow for e- 
dependent initial data ip^ with tangential kinetic energy of order one instead 
of order e^. 

In order to actually implement the scaling in the normal directions, we will 
now construct a related problem on the normal bundle of C by mapping NC 
diffeomorphically to the tubular neighborhood S of C in a specific way and 
then choosing a suitable metric ^ on NC (considered as a manifold). On the 
normal bundle the scaling of the potential in the normal directions is straight 
forward. The theorem we prove for the normal bundle will later be translated 
back to the original setting. On a first reading it may be convenient to skip 
the technical construction of 'g and of the horizontal and vertical derivatives 
V'^ and and to immediately jump to the end of Definiton [H 

The mapping to the normal bundle is performed in the following way. There 
is a natural diffeomorphism from the 5-tube B to the (5-neighborhood Bs of 
the zero section of the normal bundle NC. This diffeomorphism corresponds 
to choosing coordinates on B that are geodesic in the directions normal to C. 
These coordinates are called (generalized) Fermi coordinates. They will be 
examined in detail in Section 14. 2[ In the following, we will always identify C 
with the zero section of the normal bundle. Next we choose any diffeomor- 
phism $ G C°°(]R, (—(5,5)) which is the identity on (—5/2,5/2) and satisfies 

VjGN 3Cj-<oo VrGM: \¥^\r)\ < Cj (1 + r\^^+^^'^ (8) 

(see Figured]). Now a diffeomorphism $ G C°°{NC, B) is obtained by first ap- 
plying $ to the radial coordinate on each fibre NqC (which are all isomorphic 
to M*^) and then using Fermi charts in the normal directions. 
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Figure 1: $ converges to ±5 like 1/r. 



The important step now is to choose a suitable metric and corresponding 
measure on NC. On the one hand we want it to be the puUback $*G of G 
on Bs/2. On the other hand, we require that the distance to C asymptotically 
behaves like the radius in each fibre and that the associated volume measure 
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on NC \ Bs is dfi (g) du, where du is the Lebesgue measure on the fibers of NC 
and dfi du is the product measure (the Lebesgue measure and the product 
measure are defined after locally choosing an orthonormal trivializing frame 
of NC; they do not depend on the choice of the trivialization because the 
Lebesgue measure is isotropic). The latter two requirements will help to 
obtain the decay that is needed to translate the result back to A. 
A metric satisfying the latter two properties globally is the so-called Sasaki 
metric which is defined in the following way (see e.g. Ch. 9.3 of [1]): The 
Levi-Civita connection on A induces a connection V on TC, which coincides 
with the Levi-Civita connection on (C, g), and a connection on NC, which 
is called the normal connection (see the appendix). The normal connection 
itself induces the connection map K : TNC NC which identifies the 
vertical subspace of T(^q^^)NC with NgC. Let tt : NC ^ C be the bundle 
projection. The Sasaki metric is then given by 



It was studied by Wittich in [?Sj| in a similar context. The completeness of 
{NC, g^) follows from the completeness of C (see the proof for TC by Liu 
in [26]). C is complete because it is of bounded geometry. But {NC,g^) is, in 
general, not of bounded geometry, as it has curvatures growing polynomially 
in the fibers. However, {Br C NC, g^) is a subset of bounded geometry for 
any r < oo. Both can be seen directly from the formulas for the curvature 
in p]. Now we simply fade the pullback metric into the Sasaki metric by 
defining 



g^,,4v,w) := e(|H)G^fe.)(D$^,D$^) + (l-0(|H))4,)(t;,u;) (10) 

with \u\ := ^G$(g,o)(D<l'z/,D$i/) and a cutoff function 9 G C°°([0, oo), [0, 1]) 
satisfying = 1 on [0, 5/2] and B = on [S, oo). Then we have 



The Levi-Civita connection on {NC,'g) will be denoted by V and the vol- 
ume measure associated with 'g by dJI. We note that C is still isometrically 
imbedded and that ^ induces the same bundle connections V and V"*" on TC 
and NC as G. Since A is of bounded geometry and {Bs,g^) is a subset of 
bounded geometry, {Bs,g) is a subset of bounded geometry. Furthermore, 
{NC,g) is complete due to the metric completeness of {Bs, $*(?) (implied by 
the bounded geometry of A) and the completeness of {NC,g^). 
The volume measure associated with g^ is, indeed, d/x ® dv and its density 
with respect to the measure associated with G equals 1 on C (see Section 6.1 



5'fg,i.)(^'^) ■= gq{^T^v,T)-nw) + G^q^o){Kv,Kw). 



(9) 




(11) 
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of [IS]). Together with the bounded geometry of {Bs,g) and {Bs,g^), which 
imphes that all small enough balls with the same radius have comparable 
volume (see [S]), we obtain that 



djj, du 



eC^{NC), -r^^>c>0, (12) 



{NC\Bs/2)uc d^®dv djj,0du 

where C^{NC) is the space of smooth functions on NC with all its derivatives 
globally bounded with respect to g. 

Since we will think of the functions on NC as mappings from C to the func- 
tions on the fibers, the following derivative operators will play a crucial role. 

Definition 1 Denote by T{S) the set of all smooth sections of a hermitian 
bundle £ and by T\^{£) the ones with globally hounded derivatives up to any 
order. 

i) Fix q E C. The fiber {NqC,g(^g Q^) is isometric to the euclidean R'^ . There- 
fore there is a canonical identification l of normal vectors at q E C with 
tangent vectors at {q, v) G NqC. 

Let V? G C^{NqC). The vertical derivative G N*C at v e NqC is the 
pullback via l of the exterior derivative of if E C^{NqC) to N*C. i.e. 

for ( G NqC. The Laplacian associated with — Jj^ ^gf^^Q^{'V^ip,'V^ip)di' is 
denoted by and the set of bounded functions with bounded derivatives of 
arbitrary order by C^{NqC). 

a) Let Sf := {{q,^p) \ q E C, ip E C^{NqC)} be the bundle over C which is 
obtained by replacing the fibers NqC of the normal bundle with C^{NqC) and 
canonically lifting the action ofSO{k) and thus the bundle structure of NC. 
The horizontal connection V*^ on Sf is defined by 

{V'^cp)iq,u) := ^ Viwis),v{s)), (13) 

us s=0 

where t G r(TC) and {w,v) G ^^([-1, l], NC) with 

w{0) = q, w{0) = T{q), & v{0) = u, V^v = 0. 
Furthermore, Ah is the Bochner Laplacian associated with V^: 

1 ifj* Ai^ifj djj, ® du = - g{V\\V\)dn® dv, 
Jnc Jnc 

where we have used the same letter g for the canonical shift of g from the 
tangent bundle to the cotangent bundle of C. 
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Higher order horizontal derivatives are inductively defined by 

m 

V"^ m := V'' ^ - y^V^ V7 ^ 

" ri,...,Tmr ri " T2,...,Tmr / ^ " T2,...,Vri7j,...,Tm r 

for arbitrary Ti,. . . ,Tm € T{TC). The set of bounded sections if of Sf such 
that T V '^^ '^^^o a bounded section for all ri, . . . , G Ty^iTC) is denoted 
byC^{C:ct{N,C)). 

Coordinate expressions for and V'' are given at the beginning of Section HI 

In the following, we consider the Hilbert space Ti := L'^(^{NC,g),dJi) of 
complex- valued square-integrable functions. We emphasize that the elements 
of Ti. take values in the trivial complex line bundle over NC. This will be 
the case for all functions throughout the whole text and we will omit this 
in the definition of Hilbert spaces. However, there will come up non-trivial 
connections on such line bundles! In addition, we notice that the Riemannian 
metrics on NC and C have canonical continuations on the associated trivial 
complex line bundles. 

The scalar product of a Hilbert space H will be denoted by ( . | . and the 
induced norm by || . H^^. The upper index * will be used for both the adjoint 
of an operator and the complex conjugation of a function. 

Instead of ([7]) we now consider a Schrodinger equation on the normal bundle, 
thought of as a Riemannian manifold (NC,^). There we can immediately 
implement the idea of squeezing the potential in the normal directions: Let 

V%q,u) = V,{q,e-'u) + W{q,u) 

for fixed real-valued potentials Vc,W E C^{C, C^{NgC)). Here we have split 
up any Q G NC as (g, u) where g G C is the base point and z/ is a vector in 
the fiber NqC at q. We allow for an 'external potential' W which does not 
contribute to the confinement and is not scaled. Then e <^ 1 corresponds to 
the regime of strong confining forces. The setting is sketched in Figure [2j 

So we will investigate the Schrodinger equation 

idt^ = H'i, := -e^^NC^ + V'^ , z^|i=o = V^o ^ (14) 

where /S.^c is the Laplace- Beltrami operator on {NC,g), i.e. the operator as- 
sociated with — J^^g{dip, dip)djl. The operator if^ will be called the Hamil- 
tonian. We note that H'^ is real, i.e. it maps real-valued functions to real- 
valued functions. Furthermore, it is bounded from below because we assumed 
Vc and W to be bounded. In Section 1.3 of [H] is shown to be selfadjoint 
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Figure 2: The width of is e but it varies on a scale of order one along C. 

on its maximal domain T>{H'^) for any complete Riemannian manifold Ai, 
thus in particular for {NC,g). Let W'^'^{NC,g) be the second Sobolev space, 
i.e. the set of all L^-functions with square-integrable covariant derivatives up 
to second order. We emphasize that, in general, W'^'^{NC,g) C V{H^) but 
W'^''^{NC,g) 7^ V{H^) for a manifold of unbounded geometry. 

We only need one additional assumption on the potential, that ensures lo- 
calization in normal direction. Before we state it, we clarify the structure of 
separation between vertical and horizontal dynamics: 

After a unitary transformation H'^ can at leading order be split up into 
an operator which acts on the fibers only and a horizontal operator. That 
unitary transformation Mp is given by multiplication with the square root 
of the relative density p := of our starting measure and the product 

measure on NC that was introduced above. We recall from (fT2l) that this 
density is bounded and strictly positive. After the transformation it is helpful 
to rescale the normal directions. 

Definition 2 Setn:= L'^{NC, dp (g) du) and p := . 

i) The unitary transform Mp is defined by Mp : H ^ H, ip ^ p^'^ip. 

ii) The dilation operator is defined by {Deil)){q,h') := e^^l"^ "^{ci^v j e). 
Hi) The dilated Hamiltonian Hg, and potential are defined by 

H, := D*M;H'MpD,, V, := D*M;V'MpD, = V. + DIWD^. 

The index e will consistently be placed down to denote dilated objects, while 
it will placed up to denote objects in the original scale. 

The leading order of ife will turn out to be the sum of — Ay-I- ^(g, ■) + W{q, 0) 
and — Ah (for details on Mp and the expansion of see Lemmas [1] & [6] 
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below). When — e^Ah acts on functions tliat are constant on eacli fibre, it 
is simply the Laplace-Beltrami operator on C carrying an e^. Hereby the 
analogy with the Born-Oppenheimer setting is revealed where the kinetic 
energy of the nuclei carries the small parameter given by the ratio of the 
electron mass and the nucleon mass (see e.g. [35]). 

We need that the family of g-dependent operators — A^ + Vc{q, ■) + W{q, 0) 
has a family of exponentially decaying bound states in order to construct a 
subspace of states that are localized close to the constraint manifold. The 
following definition makes this precise. We note that the conditions are 
simpler to verify than one might have thought in the manifold setting, since 
the space and the operators involved are euclidean! 

Definition 3 Let Hfiq) := L'^{NgC,du) andVo{q,u) ■=V^{q,u) + W{q,0). 
i) The selfadjoint operator {H{{q), H'^{NqC,di')) defined by 

H({q) := -A. + Vo{q, .) (15) 

is called the fiber Hamiltonian. Its spectrum is denoted by cr(iff(g)). 
a) A function E{ : C ^ W is called an energy band, if Ef{q) G a(^Hf{q)) for 
all q ^ C. Ef is called simple, if E{{q) is a simple eigenvalue for all q EC. 
Hi) An energy band E^ : C ^ is called separated, if there are a constant 
Cgap > and two bounded continuous functions /-t : C — > M defining an 
interval I {q) := [f-{q), f+{q)] such that 

EM) = /(g)na(iJf(g)) , inf dist(a(i/f(g)) \Ef(g), E^)) = Cgap- (16) 

iv) Set {v) := a/T+TH^ = ^1 + g{qfi){y, J^)- A separated energy band E{ is 
called a constraint energy band, if there is Aq > such that the family of 
spectral projections Pq : C ^ £(7if(g)) corresponding to Ef satisfies 

Remark 1 Condition Hi) is known to imply condition iv) in lots of cases 
(see I21j for a review of known results), in particular for eigenvalues below 
the continuous spectrum, which is the most important case in the applications. 
Besides, condition Hi) is a uniform but local condition (see Figure\B\). 

The family of spectral projections Pq : C ^ C((Hf{q)) associated with a 
simple energy band E{ corresponds to a line bundle over C. If this bundle 
has a global section (pf : C ^ (o') of normalized eigenfunctions, it holds for 
all g G C that {Poip){q,u) = Mq, ■)\i^iq, ■))Hi{g) ^liQ,^)- Furthermore, ipi 
can be used to define a unitary mapping between the corresponding subspace 
PoH and L'^{C,dfiy. 
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Figure 3: E{{q) has to be separated by a local gap that is uniform in q. 



Definition 4 Let the eigenspace bundle corresponding to a simple constraint 
energy band Ef admit a smooth global section y^f : C — *■ Ti.f{q) of normalized 
eigenf unctions. The partial unitary operator Uq : Ti ^ L'^{C,dfi) is defined 
by (f/o^)(g) := {ifiiq, ■)mq, Then U*Uo = Pq and UqU* 

given by {U^^/j){q, i^) = (pf{q, i^)^/j{q). 



1 with Uq 



So any ip G PoH has the product structure ip = iUoip)ipi. Since Vb and there- 
fore yjf depends on g, such a product will, in general, not be invariant under 
the time evolution. However, it will turn out to be at least approximately 
invariant. For short times this follows from the fact that the commutator 
[if£,Po] = [~^^'^h5-Po] + C'(e) is of order e. For long times this is a conse- 
quence of adiabatic decoupling. 

On the macroscopic scale the corresponding eigenfunction -D^v^f is more and 
more localized close to the submanifold: most of its mass is contained in the 
e-tube around C and it decays like e~'^°'^'/'^. This is visualized in Figure |H 



fVo(q,v) 



Vo(q,v/8) 




0(1) ^ 0(8) 

Figure 4: On the macroscopic level yjf is localized on a scale of oder e. 
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Our goal is to obtain an effective equation of motion on the submanifold for 
states that are locahzed close to the submanifold in that sense. More pre- 
cisely, for each subspace PqH corresponding to a constraint energy band Ei 
we will derive an effective equation using the map Uq. However, in order to 
control errors with higher accuracy we will have to add corrections of order e 
to PqH and Uq. 

1.2 Comparison with existing results 

Since similar settings have been considered several times in the past, we want 
to point out the similarities and the differences with respect to our result. 
We mostly focus on the papers by Mitchell [52] and Froese-Herbst [I?] , since 
|32j is the most general one on a theoretical physics level and [T7| is the 
only mathematical paper concerned with deriving effective dynamics on the 
constraint manifold. Both works deal with a Hamiltonian that is of the form 

= -Anc + + W. (17) 

The confining potential is chosen to be the same everywhere on C up to 
rotations, i.e. in any local bundle chart (g, u) there exists a smooth family of 
rotations R{q) G SO (A;) such that 

V:{q,u) = V,iq,6-'u) = V^o. e'^ R{q)u) 

for some fixed point go on C. As a consequence, the eigenvalues of the result- 
ing fiber Hamiltonian Hi{q) = — + Vc{q,-) are constant, E{{q) = Ef. As 
our Theorems [1] and O the final result in [32] and somewhat disguised also 
in [T7] is about effective Hamiltonians acting on L'^{C) which approximate 
the full dynamics on corresponding subspaces of L'^{NC). In the following we 
explain how the results in [TT . 32j about f[T7l) are related to our results on the 
seemingly different problem f[T^ . It turns out that they indeed follow from 
our general results under the special assumptions on the confining potential 
and in a low energy limit. 

To see this and to better understand the meaning of the scaling, note that 
when we multiply H"^ by e^, the resulting Hamiltonian 

e^H' = -£^Anc + + 

is the same as H'^ in (IT^ . however, with very restrictive assumptions on the 
confining part K and with a non-confining part of order e^. As one also has 
to multiply the left hand side of the Schrodinger equation f|T4l) by e^, this 
should be interpreted in the following way. Results valid for times of order 
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one for the group generated by would be valid for times of order e"^ for 
the group generated by e^H'^ . On this time scale our result still yields an 
approximation with small errors (of order e). Thus the results in [T71 132] are 
valid on the same physical time scale as ours. 

We look at ([Hj) for initial data with horizontal kinetic energies (V^ol "^^^hV'o) 
of order one. This corresponds to horizontal kinetic energies {iPq \ — AhV'o) 
order in (1171) . i.e. to the situation where potential and kinetic energies 
are of the same order. However, in [T7[ [271 E2] it is assumed that horizontal 
kinetic energies are of order one, i.e. smaller by a factor than the potential 
energies. And to ensure that the horizontal kinetic energies remain bounded 
during the time evolution, the huge effective potential e~^E{{q) given by 
the normal eigenvalue must be constant. This is achieved in [T71 [23, [32] by 
assuming that, up to rotations, the confining potential is the same everywhere 
on C. 

Technically, the assumption that (in our units) {ipol — e'^A^iPq) is of order 
simplifies the analysis significantly. This is because the first step in proving 
effective dynamics for states in a subspace PqH. for times of order is to 
prove that it is approximately invariant under the time evolution for such 
times. Now the above assumption implies that the commutator [H^, Pq] is 
of order e^, and, as a direct consequence, that the subspace PqH. is approxi- 
mately invariant up to times of order e~^, 

\\[e-'''^\Po]\\ = Oie'\t\). 

To get approximate invariance for times of order one needs an additional 
adiabatic argument, which is missing in [32]. Still, the effective Hamilto- 
nian in [32j is correct for the same reason that the textbook derivation of 
the Born-Oppenheimer approximation is incomplete but yields the correct 
result including the first order Berry connection term. In [1^ it is observed 
that one either has to assume spherical symmetry of the confining potential, 
which implies that [H^, Pq] is of order e^, or that one has to do an additional 
averaging argument in order to determine an effective Hamiltonian valid for 
times of order £~^. For our case of large kinetic energies the simple argument 
just gives 

||[e-'^^*,Po]|| = 0(^|t|). 

Therefore we need to replace the adiabatic subspaces PqH. by so called super- 
adiabatic subspaces PgH, for which || [e"^-'^^*, P^] || = C(e^|t|), in order to pass 
to the relevant time scale. 

We end the introduction with a short discussion on the physical meaning of 
the scaling. While it is natural to model strong confining forces by dilating 
the confining potential in the normal direction, the question remains, why 
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in f[T7|) there appears the factor e^"^ in front of the confining potential, or, in 
our units, why there appears the factor in front of the Laplacian in (1141) . 
The short answer is that without this factor no solutions of the corresponding 
Schrodinger equation would exist that remain e-close to C. Any solution 
initially localized in a e-tube around C would immediately spread out because 
its normal kinetic energy would be of order allowing it to overcome 
any confining potential of order one. Thus by the prefactor e'"^ in (1171) the 
confining potential is scaled to the level of normal kinetic energies for e- 
localized solutions, while in f|T^ we instead bring down the normal kinetic 
energy of e-localized solutions to the level of the finite potential energies. 
The longer answer forces us to look at the physical situation for which we 
want to derive asymptotically correct effective equations. The prime exam- 
ples where our results are relevant are molecular dynamics, which was the 
motivation for [271 EH [32], and nanotubes and -films (see e.g. [7j). In both 
cases one is not interested in the situation of infinite confining forces and 
perfect constraints. One rather has a regime where the confining potential is 
given and fixed by the physics, but where the variation of all other potentials 
and of the geometry is small on the scale defined by the confining potential. 
This is exactly the regime described by the asymptotics e <C 1 in f|T^ . 

2 Main results 

2.1 Effective dynamics on the constraint manifold 

Since the constraining potential Vc is varying along the submanifold, the 
vertical and the horizontal dynamics do not decouple completely at leading 
order and, as explained above, the product structure of states in PqH is 
not invariant under the time evolution. In order to get a higher order ap- 
proximation valid also for times of order e^"^, we need to construct so-called 
superadiabatic subspaces PeH. These are close to the adiabatic subspaces 
PqH. in the sense that the corresponding projections Ps have an expansion in 
e starting with the projection Pq. 

Furthermore, when there is a global orthonormal frame of the eigenspace 
bundle defined by Po{q), the dynamics inside the superadiabatic subspaces 
can be mapped unitarily to dynamics on a space over the submanifold only. 

Remark 2 Let Ef be a simple constraint energy band (see Definitionl^. If 
C is contractible or if E{ is the ground state energy of Hi, then the associated 
eigenspace bundle has a smooth global section ipi : C ^ 'Hi{q) of normalized 
eigenf unctions. 
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To see this we notice that, on the one hand, all bundles over a contractible 
manifold are trivializable. On the other hand, E{ has to be an eigenvalue for 
all q due to the gap condition and the eigenfunctions of Hi{q) can be chosen 
real- valued because Hi{q) is a real operator for all q & C. Furthermore, the 
groundstate of a Schrodinger operator with a bounded potential can always 
be chosen strictly positive (see [M])- This defines an orientation on the real 
eigenspace bundle and thus a trivialization on the real line bundle. 

In the following, we restrict ourselves here to a simple energy band, i.e. with 
one-dimensional eigenspaces because we do not want to overburden the result 
about the effective Hamiltonian (Theorem [2]). To circumvent the possible 
non-triviality of the eigenspace bundle, we simply assume the existence of a 
trivializing frame for our main results. Generalizations to non-trivial bundles 
are in preparation. 

Theorem 1 Let K, W G C~(C, C^{NqC)) and let Ef be any associated sim- 
ple constraint energy band whose eigenspace bundle has a smooth section of 
normalized eigenfunctions. 

Fix E < oo. Then there are C < oo and Eq > which satisfy that for 
all e < 6o there are 

• a closed subspace P'^Ti C with orthogonal projection , 

• a Riemannian metric gl^ on C with associated measure dfi^^, 

• U' ■.n^n^s := L\C, djjl^) with U'*U' = P' and U'U'* = 1, 
such that [Hl^ := U^H^U^*, U^V{H^)) is self-adjomt on and 

||(e-'^'*-f/"*e-'^-ff*t/")P"x(^') < Ce^\t\ (18) 

for allt eM. and each Borel function x : K ^ [— 1, 1] with supp x C {—oo,E]. 
Here x{H'^) is defined via the spectral theorem. 

The proof of this result can be found in Section 13.11 The estimate f|T8|) 
means that, after cutting off large energies, the superadiabatic subspace P^Ti 
is invariant up to errors of order e^\t\ and that on this subspace of Ti the 
unitary group e~^^^^ is unitarily equivalent to the effective unitary group 
e~'^cff* on L'^lCjdfilg) with the same error. In particular, there is adiabatic 
decoupling of the horizontal and vertical dynamics. 

The energy cutoff x{H'^) is necessary in order to obtain a uniform error esti- 
mate, since the adiabatic decoupling breaks down for large energies because 
of the quadratic dispersion relation. It should be pointed out here that, while 
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P'^x{H'^) is not a projection, it is really the difference of the unitary groups 
that is small and not P^x(if^). More precisely, P^x(if^) can be replaced 
by U^* x{H^q)W , which is a projector for any characteristic function x, as is 
done in the following corollary. 

Before we come to the form of the effective Hamiltonian, we state our result 
about effective dynamics for A, which follows from the one above. 

Definition 5 Set Aip := (^|;)~^^^ {ipo^) with $ : NC Bs as constructed 
in Section [7771 and $*(ir the pullback of dr via $. This defines an operator 
A G C{L\A,dT),n) with AA* = 1. 

The stated properties of A are easily verified by using the substitution rule. 
Of course, the choice of our metric (ITU]) changes the metric in a singular way 
because it blows up a region of finite volume to an infinite one. However, 
it will turn out that the range of consists of functions that decay faster 
than any negative power of \C\/£ away from the zero section of the normal 
bundle. Then leaving the metric invariant on Bs/2 is sufficient; due to the 
fast decay the error in the blown up region will be smaller than any power 
of e ioT e <^ 5. 

Corollary 1 Fix 5 > and E < oo. Let := — £:^A_4 + be self-adjoint 
on L'^{A,dT). Assume that D*^AVXA*D, = K + D*^WD, for some V^,W 
satisfying the assumptions of TheoremUi Then there are C < oo and £o > 
such that 

||(e-i^5.*_^*f/-e-^^^«*f/e^)^*f/-^(/f^^^)f/^^||^^^^ < Ce'\t\ 

for all < e < 6o, t E M, and each Borel function x : M ^ [—1, 1] with 
suppx C {-oo,E]. 

The proof of this result can be found in Section 13.21 We note that for 
any charateristic function x the operator A* f/^*x(if|g)f/^y4 is a projector 
whose image is canonically identified with x{H^s)L'^{C,dfilQ). The assump- 
tion D*AVXA*D^ = Vc + D*WD^ means that should be of the same 
form as in the 5-tube of C. Moreover, the assumptions on Vc and W in 
Theorem [1] translate into assumptions on on the (5-tube only. 

2.2 The effective Hamiltonian 

Here we write down the expansion of the effective Hamiltonian i^eff- We 
do this only for states with high energies cut off. Then the terms in the 
expansion do not depend on any cutoff, which is a non-trivial fact, since we 
will need cutoffs to construct i^eff! 
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Theorem 2 In addition to the assumptions of TheoremUl assume that the 
global family of eigenfunctions (pf associated with E{ is in (C,H{{q)) . 

(2) 

For all e small enough there is a self-adjoint operator H^^ on Ti.cs such 
that for each Borel function x '■ ^ ^ [~1> 1] '"^^^^ supp x C {—oo,E], for 
every ^ G {U''x{H'')U^*,x{HlB),x{His)}, and for all e H^s satisfying 
^ = x(-^'Ac + Ei)^lj zt holds that \\ {Hlg - H^^) ^ \\c[n,^) = 0{e^) and 

[aUiPl^^t^r.Plsi^) + cl>*{E, + e{if,\{V''W)^,U+e^W^^^)ij 

-e'M {^*{eWes<P, Pls<P, 0), ^i^Wesi^, Plsi^, ^))) df^l^, 
where for Ti,T2 G T(T*C) 

9ts{n,r2) = g{Ti,T2) + e {ipi\2ll{.){Ti,T2)(p{)'Hf 

pIqiP = -iedij - Im('£((/^f|VVf)wf - / i'/'f R(V>f,i/)i/c/z/ 

^ JNqC 

+ e'{^i\2{Wi.) - (v9f|W(.)v'f)7i,)VVf>„Jv^, 

with W the Weingarten mapping, II the second fundamental form, R the 
curvature mapping, TZ the Riemann tensor, and Tq*^C and Ng*^C canonically 
included into t'^^^-^NC. The arguments ' . ' are integrated over the fibers. 

Furthermore, W^^^) = {ip^\\{yiW)ipi)H, + 14h + V^geom + Kmb and 







- Po) VVf) du, 












- i (Tt + trcRic + trc^), 






V^V^f , u) du, 




JNqC 


M{^\,^2) = 


{^^\{l-Po){Hi 


-E,y\i-Po)cp2,)^^ 




-iff trc(W(z/)A) 





with f] the mean curvature vector, k, 7t the scalar curvatures of C and A, and 
tic Ric, tie TZ the partial traces with respect to C of the Ricci and the Riemann 
tensor of A (see the appendix for definitions of all the geometric objects). 
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This result will be derived in Section 13.31 One might wonder whether the 
complicated form of the effective Hamiltonian renders the result useless for 
practical purposes. However, as explained in the introduction, the possibly 
much lower dimension of C compared to that of A outweighs the more com- 
plicated form of the Hamiltonian. Moreover, the effective Hamiltonian is of a 
form that allows the use of semiclassical techniques for a further analysis. Fi- 
nally, in practical applications typically only some of the terms appearing in 
the effective Hamiltonian are relevant. As an example we discuss the case of 
a quantum wave guide in Section [231 At this point we only add some general 
remarks concerning the numerous terms in Hes and their consequences. 

Remark 3 i) If C is compact or contractible or if Ef is the ground state 
energy of H{, the assumption Vq G C^(C,C^{NgC)) implies the extra 
assumption that ipf G C^{C,TC{) (see Lemma [77] in Section \J7S^ . We 
do not know if this implication holds true in general, but expect this for 
all relevant applications. 

a) Vf^ip •= (iPefrV')(''") o metric connection on the trivial complex line 
bundle overC where ip takes its values, a so-called Berry connection. It 
is flat because iff can be chosen real-valued locally which follows from Hf 
being real. The first order correction in p^^ is a geometric generalization 
of the Berry term appearing in the Born- Oppenheimer setting. When 
the constraining potential is constant up to rotations, the first- order 
correction reduces to the Berry term discussed by Mitchell in l3^ . 

Hi) The correction of the metric tensor by exterior curvature is a feature 
not realized before because tangential kinetic energies were taken to be 
small as a whole. Its origin is that the dynamics does not take place 
exactly on the submanifold. Therefore the mass distribution of iff has 
to be accounted for when measuring distances. 

iv) The off-band coupling M. and V^n, o,n analogue of the so-called Born- 
Huang potential, also appear when adiabatic perturbation theory is ap- 
plied to the Born-Oppenheimer setting (see J3^). However, M. contains 
a new fourth order differential operator which comes from the exterior 
curvature. Both M. and Vbh can easily be checked to be gauge-invariant, 
i.e. not depending on the choice of iff but only on Pq. 

v) The existence of the geometric extra potential V^eom has been stressed 
in the literature, in particular as the origin of curvature-induced bound 
states in quantum wave guides (reviewed by Duclos and Exner in [T3f). 
In our setting, these are relevant for sending signals over long distances 
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only (see Remark\^ below). A simple example where the inner curvature 
of the ambient manifold plays a role was given by Freitas and Krejciiik 
in Section 5 of 115]. The potential V^mb wcls also found in 



vi) If H^g was defined by the expression in the theorem, the statement 
would be wrong for = x{H^q ) because the fourth order term in Ai 
would be dominant. Therefore Ai is modified in the definition of H^^ 
so that the associated operator is bounded (see 137]) below). However, 
when energies of H^^ are approximated by perturbation theory or the 
WKB method, that modification is of lower order as the leading order 
of a quasimode ip satisfies ip = xi—^'^^c + Ef)ip + 0{e) for some x- 

Using Theorem [2] we may exchange if^g with Hj^^ in Theorem [H After 
replacing and If^ by their leading order expressions, which adds a time- 
independent error of order e, it is not difficult to derive the following result. 

Corollary 2 Fix E < oo and set Uq := UqD*. Under the assumptions of 
Theorem\^ there are C < oo and Eq > Q such that for all e < Eq, t & and 
each Borel function % : R — [— 1, 1] with supp x C {—oo^E] it holds 

(e-^^*-t/o^*e-^-*f/o^)f/o"x(ifSVo^ < Ce{e'\t\ + l). (19) 

Corollary [2] will also be proved in Section [3^ While (|T9|) is somewhat weaker 
than (1181) . it is much better suited for applications, since Uq is given in terms 
of the eigenfunction ip^ and depends on e only via the dilation D^. So, in view 
of Theorem [21 all relevant expressions in f|T9l) can be computed explicitly. 



2.3 Approximation of eigenvalues 

In this section we discuss in which way our effective Hamiltonian allows us 
to approximate certain parts of the discrete spectrum and the associated 
eigenf unctions of the original Hamiltonian. The following result shows how 
to obtain quasimodes of from the eigenf unctions of and vice versa. 

Theorem 3 Let Ef be a simple constraint energy band whose associated 

(2,) 

eigenvalue bundle is smoothly trivializable, and let U^,H^g be the operators 
associated with E{ via TheoremsU\ & \E 

a) Let G M. Then there are Eq > and C < oo such that for any family 
(E^) with limsup^^gE'g < E and all e < Eq the following implications hold: 

i) i/SVe = E,^, =^ II [H^ - E,) U'*^e \\n < Ce' WU'^^sWh, 

ii) H^r = E,r =^ iK^iff^ - E,)u'r\u, < ce^wrWn- 
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b) Let E{{q) = inf o"(iff (g)) for some (and thus for all) q E C and define 
Ei{q) := inf cr(iJf(g)) \ Ef{q). Let {ip^) he a family with 

limsup (V^"|(-e2MpAvM; + -l/o(g,'^/e))^'> < mf^i- (20) 

Then there are eg > 0, c > such that \\U^il)'^\\-H^g > c W'ip'^W^ for all e < Eq. 

We recall that for any self-adjoint operator H the bound ||(i/ — A)-?/^!! < S\\ip\\ 
for A G M implies that H has spectrum in the interval [A — 5, A + 5] . So a) i) 
entails that if^ has an eigenvalue in an interval of length 2Ce^ around E^, if 
one knows a priori that the spectrum of if^ is discrete below the energy E. 
The statement b) ensures that a) ii) really yields a quasimode for normal 
energies below inf^gc-^i, i.e. that 

H'^' = E,r =^ \\{H^;^ - E,)U'i,'\\n^^ < ^ \\U' \\^^^ . 

Remark 4 If the ambient manifold A is flat, —e^MpA^M* is form-hounded 
by —e'^Afyfc + Ce"^ for some C < 00 independent of e (this follows from 
Lemma\J\ below and the expression (5.5) for —e^A^c W^)- Then, since 
= -e^Ajvc + Vo{q, u/e) + W{q, u) - W{q, 0), ^ follows from 

limsup(V^"|i/"V'') < inf El -sup (Vr(g,0) - l^(g,z/)) =: E,. 

Therefore Theorem\^ in particular, implies that at least for flat A there is a 
one-to-one correspondence between the spectra of and below E^,. 

One may ask whether a family {E^) of energies of if^ with limsup-E^ < E^, 
exists at all. A sufficient condition is that sup(q_^) {W{q, 0) — W{q, u)) is 
strictly smaller than inf i^i — inf^^c i?f . For this implies mfg^cEf < E^, 
and the spectrum of H^g in the interval [inf^gc E{, E^] has either a continuous 
part or the number of eigenvalues is at least of order because H^^^s 
leading order term —e'^Ac + -Ef is a semiclassical operator. Then by a) i) this 
is also true for H'^. 

(2) 

The eigenvalues of H^^ can be approximated by the WKB construction, 
which is quite standard (see e.g. [20]). In the simplest case one obtains: 

Corollary 3 Assume that A is flat and let E{ is a constraint energy band 
with inf Ef < E^ and Ef{q) = inf a[H{(q)) for all q E C. Let there be qo such 
that Ei{qo) < Ef{q) for all q qo and (V| , g .EA{qo) is positive definite. 
Denote by Ef {A) the i-th eigenvalue of a semi-bounded operator A, counted 
from the bottom of the spectrum. Then for any £ G N 

Ee{H') = Eiiqo) + eEe{Hno) + 0{e^), 

where Hhq := —A^d + ^(Vg . q ,E{){qQ)x^x^ is a harmonic oscillator on M'^. 
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We want to relate this to results by Friedlander and Solomyak in and 
Borisov and Freitas in [3J. In both works the spectrum of the Dirichlet 
Laplacian —Ad on the thin domain {{x,y) : x G Q, —eh^{x) < y < 
for positive functions /i± : f2 — >■ M is considered, where Q is any interval in M 
in [16] and a bounded domain in M'' in p]. It is shown that, if /i := /i+ — /i_ 
has a global maxmimum at Xq E Q which is non-degenerate, then 

Eei-e'A^) = Ef(xo) S^.^. + |(V^ ^,«.Ef)(a;o)xV) +0(5), 

where Ef{x) := n'^/h'^{x) is the lowest eigenvalue of the Dirichlet Laplacian 
on the interval [0, h{x)]. So our result generalizes this to an Q that is curved 
and of arbitrary codimension, but with the Dirichlet Laplacian replaced by a 
constraining potential. For a set with smooth boundary, however, we do not 
see a problem in using the Dirichlet Laplacian instead of the constraining 
potential. The strict localization to an e-tube around C would even simplify 
many steps in our proof considerably. 

2.4 Application to quantum wave guides 

In this section we look at the special case of a curve C in ^ = M'^ equipped 
with the euclidean metric. Such curves may model quantum wave guides 
which have been discussed theoretically for long times but are nowadays also 
investigated experimentally. We will provide the expression for our effective 
Hamiltonian in this case and discuss which terms remain, if we add restric- 
tions on the constraining potential or the geometry. Furthermore, we will 
apply Theorem [3] to obtain a statement about the spectrum of a twisted 
wave guide. For the sake of brevity, we assume that W, the non-constraining 
part of the potential, vanishes. Its contributions could be trivially added in 
all formulas and as long as supj-^ {W{q,0) - W{q,u)) is small enough also 
in the statements (see the preceding subsection). 

We first look at an infinite quantum wave guide. So let the curve C be given as 
a smooth injective c : R — >■ M^, x ^— c{x) that has bounded derivatives of any 
order and is parametrized by arc length (|c| = 1). The mean curvature vector 
of cis rj = c and its (exterior) curvature is \ri\. Denoting by ■ the usual scalar 
product in we define n{i>) := ty-rj/lril where rj ^ and n(z/) := elsewhere. 
By the Frenet formulas the Weingarten mapping satisfies VV(?7) = |?7p (see 
e.g. [Sj) and W = on the orthogonal complement of 7] (which is meant to 
be NgC if r]{q) = 0). 

A normalized section of TC is given by r := c. We extend this to an or- 
thonormal frame of TC x NC in the following way: We fix g G C, choose an 
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arbitrary orthonormal basis of NqC^ and take z/i, z/2 to be the parallel trans- 
port of this basis with respect to the normal connection V"*" along the whole 
curve. This yields an orthonormal frame of NC. Together with r we obtain 
an orthonormal frame of TC x NC, which is sometimes called the Tang frame. 
We denote the coordinates with respect to r, i/i, and 1^2 by a:, ni, and n2 
respectively. In these coordinates it holds V'^ = (as can for example be 
seen from the general coordinate formula (1701) below). 

Now let Ei be a simple constraint energy band and yjf a global family of 
eigenf unctions associated with it. We start by spelling out the formula for 
ifgg from Theorem [2J Of course, all terms containing the inner curvature 
of C and ^ = M'^ vanish due to the flatness of C and A with the euclidean 
metric. Since C is one-dimensional and contractible, yjf can be chosen such 
that = —iedx globally. Then the effective Hamiltonian is 

^qwg = -edx{l + e\T]\{ifi\n^i) +?,e^\r]\^{ifi\n^^i))edx + Ei 

+ (Aedx {dx(Pi\RHf{Ef)d:,(Pf) ed^, 

+ 4:\r]\ Reed,j: {dx(pi\RHtiEi)nipi) e^dl^ 

+ H'^e^dl, {n^,\RH,{E,)n^,)e^dl)i (21) 

with RnAEi) ■■= {1 - Po){Hf - Eiy\l - Po) and := J^, (j)* ^ dnidn2. 

We emphasize that formula (^T\i is only valid when applied to states with 
high energies cut off because this is required for the application of Theorem [21 
In particular, this explains why the differential operator of fourth order is 
not to be thought of as the dominant term but only as of order e"^. But still 
ll^i'^ll ~ is possible for a of finite energy! Before we consider some 
special cases, we want to make the following crucial remark about sending 
signals through wave guides with constant 'cross section'. 

Remark 5 For highly oscillating states ip, i.e. with {ip\ — e'^d'^^ijj) ~ 1, the 
only term of order one besides —s^d^^ is E{. In particular, if E{ is constant, 
the dynamics is free at leading order and, even more, the potential terms are 
of order e"^. So they only become relevant for times of order e'"^. However, a 
semiclassical wave packet ip covers distances of order e^^ on this time scale. 
Hence, for such ip note-worthy trapping occurs only for very long wave guides! 

If we consider a straight wave guide, i.e. ?7 = 0, the formula we end up with is 
a complete analogue of the one derived by Panati, Spohn, and Teufel in |3H| 
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in the case of the Born-Oppenheimer approximation: 

+ Aed^, {d,ipi\RH,{Ef)d,^i) ed,. (22) 

We note that, although rj = 0, the x-dependence of the constraining potential 
still allows us to model interesting situations, e.g. a beam splitter [22j. 
Now we drop the assumption r] = and assume that the constraining poten- 
tial Vc is parallel with respect to V*^ instead. This means V^K = 0. Then 
we obtain a global family of eigenfunctions ipf with V^ipf = by taking it 
to be the parallel transport with respect to of (p{{qo) for any go ^ C. In 
addition, since Vc does not change its shape, E{ is constant and thus may be 
removed by redefining zero energy. Therefore we have 

^qwglvhVc^o = -ed:,{l + e\r]\{(p{\n(pf) +3e'^\r]\'^{(p{\n'^(p{))ed^ 

-e'lvWi + e'\v\'e'dl{n^,\RHAEi)ny,,)e'dl. (23) 

There is a wide literature on quantum wave guides where the effects of bend- 
ing and twisting on the spectrum of the Dirichlet Laplacian on an e-tube 
with a fixed cross section is investigated (see the review [21] by Krejcink). 
If we consider the corresponding situation that Vc does not change its shape 
but is allowed to twist, Ef is the only term in fl2T]) that may be neglected. 
Since the remaining potential terms are, however, of order e'^, the kinetic 
energy operator —e^d^^ will also be of order for low eigenvalues. So H^„^ 
may be devided by e"^. Keeping only the leading order terms we arrive at 

^twi.t := -dL - \V\V^ + {d.^ild.^i) - \{^f\d.^i)\' (24) 

The twisting assumption means that there is Vc G C^(M^) and a G C{f (R) 
such that the constraining potential has the form: 

(y^{x)){ni,n2) := K(ni cosq;(2;) — ^2 sina(a;), rii sina(x) + ^2 cosa(a;)). 

Then the family of eigenfunctions ipf may be chosen as 

(yjf (x)) (ni, 712) := $f (ni cos a(x) — ^2 sin q;(x), ni sin q;(x) + n2 cos q;(x)) 

for an eigenfunction $f of —Ar2 + Vc{x) with eigenvalue E^. A simple calcu- 
lation yields 

{dxVi\dx<^i) - \{Vf\dx<^i)\'^ = dM \nidn.2^i - n2dn^^i\^dnidn2. 

We note that the integral is the expectation value of the squared angular 
momentum of $f and thus vanishes for a rotation- invariant $f. So fl24p 
shows that bending is attractive, while twisting is repulsive. Now Theorem [3] 
together with Remark [2] and Remark H] implies the following result. 
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Corollary 4 Let C gM.^ be an infinite curve, W = 0, and E{ be a constraint 
energy band with E{{q) = inf cr(iff(g)) for all q E C. Let L G NoU{oo} be the 
number of eigenvalues of //twist below the continuous spectrum, where i/twist 
is the following operator on M; 

^twist := -C - 1^7174 + C($f)«2 

with C($f) := /jg2 |ni(9„2<l>f - n2dn^(^i\^dnidn2. 

Denote by Ei{A) the £-th eigenvalue of a semi-bounded operator A, counted 
from the bottom of the spectrum. If Vc only twists, i.e. Vc = for some a 
as above, then for any i < L 

Ee{H') = Ei + ^'^.(i/twist) + 0{e'). 

This is an analogue of the result by Bouchitte, Mascarenhas and Trabucho 
in [6] for £-tubes twisted by a, extended by Borisov and Cardone in [2]. In [2l] 
it was posed as an open problem to generalize this result to an infinite tube. 
Corollary H] achieves this for a constraining potential that twists instead of 
the Dirchlet tube. 

Up to now we have considered an infinite wave guide which is topological 
trivial. The only possible non-trivial topology for a one-dimensional manifold 
is that of a circle. So let C now be diffeomorphic to a circle. We refer to 
such a C as a quantum wave circuit. Because of the non-trivial topology our 
choices of the family iff made above are only possible locally but in general 
not globally. Therefore we rewrite fl^2]) without those choices and ignoring 
the terms of order for the moment: 

^qwc = p:{l + e\v\{^i\ny,,))p, + E, + 0{e') (25) 

with = —iedx + E (^(pi\idx(pf) ■ Although the curvature of the connection ip^ 
always vanishes, it may have a non-trivial holonomy over the circle, which 
we will discuss next. 

For the sake of simplicity we consider a round circle, i.e. with constant |?7|. 
Let a; be a 27r-periodic coordinate for it. The eigenfunction v^f(x) can be 
chosen real-valued for each fixed x because Hf is real. This associates a real 
line bundle to Ei. From the topological point of view, there are two real 
line bundles over the sphere: the trivial one and the non-trivializable Mobius 
band. In the former case the global section iff can be chosen real everywhere. 
This implies {<fi\dxipi) = which results in ip^ = ed^. Thus the holonomy 
group of ipe is trivial in this case. We will now provide an example for the 
realization of the Mobius band by a suitable constraining potential and show 
that the holonomy group of ipe becomes Z/2Z! 
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Let Vc G C{^(]R^) have two orthogonal axes of reflection symmetry, i.e. in 
suitable coordinates 

K(-iVi,iV2) = K(iVi,iV2) = K(iVi,-iV2). (26) 

Then the real ground state $o of — A]fj2 + Vc with energy Eq is symmetric 
with respect to both reflections, 

<l>o(iVi,iV2) = $o(-iVi,iV2) = $o(iVi,-iV2), 

while the first excited state $i, also taken real- valued, with energy Ei is 
typically only symmetric with respect to one reflection and anti-symmetric 
with respect to the other one, e.g. 

$i(iVi,iV2) = -$,(-iVi,iV2) = $i(iVi,-iV2). (27) 

This is true in particular for a harmonic oscillator with different frequencies. 
As the potential constraining to the round circle we choose the twisting 
potential introduced above with a{x) = x/2, i.e. 

{V^^'^{x)) (7^1,^2) := V"c( cos(x/2)?2i — sin(x/2)n2, sin(x/2)ni + cos(x/2)ri2) • 

We note that due to ([26D this deflnes a Vc^^ E C^{C,C^{NC)). Then 

[ipj{x)) (ni, 712) := $j ( cos(a;/2)ni — sin(a;/2)n2, sin(x/2)r2i + cos(a;/2)n2) 

is an eigenfunction of Hf{x) := — Av + Vc{x) with eigenvalue Ej for every x 
and j G {0, 1}. However, while (po is a smooth section of the corresponding 
eigenspace bundle, ipi is not. For by (1271) it holds (pi{x) = — (^i(a; + 27r). Still 
the complex eigenspace bundle admits a smooth non-vanishing section. A 
possible choice is ^pi{x) := e'^/^<^i(x). Using (l27j) we obtain that for the flrst 
excited band the effective Hamiltonian (l25i) reduces to 

= E^ + i-ie^. + e/2f + Oie'), 

while for the ground state band it is 

-f^qwc.o = - e^dl,^ + 0{e^) 

This shows that depending on the symmetry of the normal eigenfunction the 
twist by an angle of vr has different effects on the effective momentum operator 
in the effective Hamiltonian. With respect to the connection appearing in 
-^qwc.i holonomy of a closed loop 7 winding around the circle once is 
given by h{^) = e'-^o 1/2^2; _ Hence, the holonomy group of ipe is indeed 
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Z/2Z and the 1/2 cannot be gauged away. Furthermore, a wave packet which 
travels around the circuit once accumulates a vr-phase. This can be seen as 
an analogue of the Aharanov-Bohm effect, though with the only possible 
phase 71. 

The effect of this phase can also be seen in the level spacing of Hj and thus, 
with Theorem [31 also in the spectrum of H"^. The arguments that led to fl24p 
for an infinite wave guide may be applied here, too, except that, of course, 
—dl^ has to be replaced by {—idx + 1/2)^ for i?q„c,i- Since \ri\ and a are 
constant, the eigenvalues of -f^q„c,i ^-^e 



C(-^i)-|'?l 



while for H^^^ q we find 



qwc,o) 



+ 



Although a constraining potential that twists along a circle was investigated 
by Maraner in detail in [27| and by Mitchell in [32] , the effect discussed above 
was not found in both treatments. The reason for this is that they allowed 
only for whole rotations and not for half ones to avoid the non-smoothness 
of (fi. Finally, we note that it easy to generalize the statements above to a 
circuit whose curvature and potential twist are non-constant. 



3 Proof of the main results 

In the following, Y) is the Banach space of bounded operators between 
two Banach spaces X and Y. T>{A) will always denote the maximal domain 
of an operator A, equipped with the graph norm. For convenience we set 
J){H^) := 7i. A* will always be used for the adjoint of A onTi. if not stated 
differently. We recall that we have set (z/) := a/1 + A = (z/)' is meant 
to be the multiplication with (z/)'. Finally, we write a < 6, if a is bounded 
by b times a constant independent of e, and a = 0{e^), if ||a|| < 

Throughout this section we assume that Vc,W E C^{C, C^{N.C)) and that 
E{ is a constraint energy band as defined in Definition [31 

3.1 Proof of adiabatic decoupling 

As explained in the introduction the first step in proving Theorem [T] is the 
unitary transformation of if^ by multiplication with the square root of the 
relative density p := , ^ , of the volume measure associated with 'g and the 
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product measure on NC. This transformation factorizes the measure, which 
will allow us to easily split the integral over NC later on, but it also yields 
an additional potential term. The abstract statement reads as follows: 

Lemma 1 Let {Ai,g) be a Riemannian manifold. Let dai,da2 be two mea- 
sures on A4 with smooth and positive relative density p := Define 

Then Mp is unitary and it holds 

M;(-A,.JM,^ = -A,.,^ - (i(7(d(lnp),d(lnp))-iA,.,(lnp))V' 

= : -Arf^jV + Vpi^^ 

with A^o-. := diVdai grad-?/', where gia-dip is the vector field associated with dip 
via g and div^o-, is the adjoint o/grad on L'^{A4,dai). 

The proof is a simple calculation, which can be found in the sequel to the 
proof of Theorem [H We recall from f|T2|) that p = -j^^ is in C^{NC) and 
strictly positive. Therefore Vp is in C^{NC) for our choice of p. Since p is 
equal to 1 outside of Bs, Vp is even in C^{C, C^{NqC)) which coincides with 
C^{NC) inside Br for any r < oo. 

The heart of Theorem [T] is the existence of a subspace PsTi. C 7i that can be 
mapped unitarily to L'^{C, dp) and approximately commutes with Hs'. 

Lemma 2 Under the assumptions of Theorem Q] there is Eq > such that 
for all e < Eq there are an orthogonal projection G C{TC) and a unitary 
lie e C{n) with Pe = U*PoUe and 

Ue-lWdH) = 0{e) , WPeWdViH^)) ^ 1, 

{i^ypeii^ywciH) < 1, m'Ps{iyy\\civiH.)) < i, 

[H,,P,]xm\\cin,viHr)) = Oi^') (28) 

for all j,l,m G Nq and each Borel function % : M ^ [^1)1] satisfying 
suppx C {-oo,E]. 

The construction of and Us is carried out in Section 14.31 There is a 
heuristic discussion at the beginning of that section that the reader may find 
instructive to get an idea why and U^ exist. When we take its existence 
for granted, it is not difficult to prove that the effectice dynamics on the 
submanifold is a good approximation. 
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Proof of Theorem [T] (Section 12.11) : 

Let d/igg be the volume measure associated with g^^ which we define by 
the expression in Theorem [2j For any fixed E < oo, Lemma [2] yields some 
unitary Us for all e below a certain Eq- Since we assumed that the eigenspace 
bundle associated with Ef is trivializable, there is ?7o : ?i — > L'^{C,dfi) as in 
Definition |H We define '■= UqUs- Using Definition H] and Lemma [2] we 
have U;Us = U;U*UoUs = U^PqUs = Pe and 

UeU; = UoUsU:U* = UoU* = 1. (29) 

In view of Lemma [H we next set := UsD*M* with p := , , and 
p ■= In view of ((29]), the unitarity of Mp, Mp, and implies U^U'* = 1. 
Furthermore, we simply define by P^ := U^*U'^ . Then is unitary from 
Pm to L'^iC.dpil^). Finally, we set 

Hl^ := U'H'U'* = MpUeHsU*M*~. (30) 

We notice that H^^ is symmetric by definition. Since Mp is unitary and 
f/g is unitary when restricted to P^Ti. due to Lemma [2], the self-adjointness 
of f/^P(if ^)) on TCes '■= L'^{C,dfilfi) is implied by the self-adjointness 

of [PsH^Pe, PeT>{Hs)) on PeH, which is in turn a consequence of the self- 
adjointness of {PeHsPe + il-Pe)Hs{l-Pe),V{He)) OU H. FoT E Small CUOUgh 

this last self-adjointness can be verified using Lemma [2] and the Kato-Rellich 
theorem (see e.g. [37|): 

- [PeHsP, + (1 - Ps)Hs{l - P,)) 

= (1 - Pe)HsPs + PMl - Ps) 
= {l-Pe)[He,Pe] " Pe[He^ Pe] 

= (l-2P,)[iJ„P,]. 

Lemma [2] entails that [ifg, Pe] is operator-bounded by eH^. Hence, for e 
small enough (we adjust Eq if nescessary) the difference above is operator- 
bounded by Hs with relative bound smaller than one. Now the Kato-Rellich 
theorem yields the claim, because (i^e, V(Hs)) is self-adjoint (as it is unitarily 
equivalent to the self-adjoint H^). 

We now turn to the derivation of the estimate f|T8|) . To do so we first 
pull out the unitaries Mp, Mp and D^. Using that D*^M* x{H^) MpD^ = 
x{D*M*H^MpDs) = x{Hs) due to the spectral theorem we obtain by a 
straight forward calculation that 

(e"^^'* - f/^*e-'^°ff*?7") P'x{H') 
= MpDe (e-'^=* - U*e-'^^"^^'-'Ue) U*Ue x{He) D*M*p. 
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Since Mp and are unitary, we can ignore them for the estimate and con- 
tinue with the term in the middle. Next we use Duhamel's principle to 
express the difference of the unitary groups as a difference of its generators. 
Because of UeU* = 1 and Pg = U*Ue we have that 



AHA 



_ lf*Q-iUeH,U*t 



Ue)U:UeXiHe 



(P,-f/;e-'^^^^^^^**?7,e^^^*)e-'^^*x(^.) + [e"^^^*, P.] 



i / u:e 



*-iU.HM*s 



+ e- 



Pe]x{He) 



J* ^—iUeHsU* s 



+ e 



-\HA 



Pe]xm, (31) 



where we used that [e xl-f^s)] = for any s due to the spectral theorem. 
Now we observe that fl28|l implies that 



-iHrt 



Pe]x{He 



\C{H) 



0{e'\t\), 



(32) 



as it holds 

[e-^^*,P,]x(i^. 



^p^ _ ^iHAp^^-iH^^ ^^^^ 



il e'''^-'{H,P,-P,H,)e-'''^'dsx{H, 



= -e-'^=''i / e^^^^[/7„P,]x(^.)e-^^^^rfs ^ Oie'\t\) 
Jo 

because of Lemma [2] and ||e~'^='^||£(-H) = 1 for any s. So, in view of fl3Tl) . 



(e-^^^* - f/>-i^^^='^^*[/e) Pe xm 



\C{H) 



< 



u:e 



*-iUeHeU; 



^'U,[H,,P,]xme'''^'ds 



< \t\ u:e 



*-iUeH^U*S 



\C{H) 



C(H) 



< 1 



0{e'\t\). 

This proves the error estimate (fTSl) . 



□ 
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Proof of Lemma [TJ 

Mp is an isometry because for all ipjip E L'^{Ai, d(j2) 

M pip* M pip da I = I ijj* (f dai = / ipipda2. 
M Jm Jm 

Therefore it is clear that 

Mpifj = p^ip 

which is well-defined because p is positive. One immediately concludes 

MpM*p = 1 = M*pMp 

and thus Mp is unitary. Now we note that [grad, p"^] = — i p^i gradlnp . 
So we have 

M;{-Ad^,)Mpi; = -p5divrf,,grad(p-^) 

= - p^divd^^ p"5 (grad?/' - Kgradlnp)?/^) 

= -p^divrf^i p~5 grad^/) + p^div^^, (^p"5i(gradlnp)?/^j 

On the one hand, 

p^divd,^^ p~5 gradT/) = pdiv^,^^ p"^ grad?/) + ^ 5f(gradlnp, grad^)) 
and on the other hand, 

p^divd^^ (^p"^(gradlnp)^j = - i ^((gradlnp, gradlnp)?/) 

+ i (divdai grad lnp)V' 
+ (grad In p, grad ip) ■ 

Together we obtain 

M;{-Ad^,)MpiP = -pdivrf^^p-igradV- 

I g (grad In p, grad In p) - | div^^^ grad In p j ip 

= -Ada^ip - (^i5((gradlnp,gradlnp) - ^ A^^^ Inpjv, 

which is the claim. □ 

3.2 Fullback of the results to the ambient space 

In this section we show how to derive Corollary [T] about effective dynamics 
on the ambient manifold A from Theorem [H To do so we first state some 
immediate consequences of Lemma [2] for and from Theorem [H 
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Corollary 5 For e small enough and from TheoremUl satisfy 

P''\\c{V{H^"')) ^ 1, 

{i^/eyPHi^/eyWcCH) < 1, U^/eyP^r^/sYWciviH^)) < 1, 
[-^^5 -P^]ll£(x>(/i'E'"+\x'(ii'=™)) ^ ^(^)' (33) 
[H',nxm\\cCHMH^-) = 0{e\ (34) 

/or a// j,l,iTL G No one? eac/i Borel functions x '■ ^ ^ satisfying 
suppx C (-00, E]. 

The proof can be found at the end of this subsection. Now we gather some 
facts about the operator A defined in ([5]) and its adjoint. 

Lemma 3 Let A be defined by Aip := (^^)~^ {ip o $) with $ : NC ^ B as 
constructed in Section [7T71 

z) It holds A e C{L'^{A,dT),n) with 

|l2 (7vc,4r) < 1 1 ^ 1 1 L2 {A,dT) y e L'^{A,dT). 



a) For Lf G H the adjoint A* G C(T-C, L'^{A, dr)) of A is given by 



onA\B. 



It satisfies \\A*(p\\L2(A,dr) = WvW L^NC,dji), A* A = xb, and AA* = 1. 
It holds A*P^ G C{V{H^),V{H%)) and 

\\{H\A* - A*H')P'\\c^T,^H^)^LHA4r)) < e\ (36) 

The last estimate is crucial for the proof of Corollary [H It results from the 
two facts that Hj^A* = A*H'^ on Bs/2 by construction and that is 'small' 
on the complement. Lemma [3] will be proved at the end of Section W7\\ The 
following estimates for energy cutoffs will be useful not only for the proof of 
Corollary [T] but also in the following sections. 

Lemma 4 Assume that (if, V{H)^ is self-adjoint on some Hilbert space Ti. 
Let Xi ^ C*^(]R) and X2 '■ ^ ^ ^ a bounded Borel function, 
a) Let A G C{n). If \\[H, A]x2{H)\\^^^^^,^^^^^^_,^^ < 6 for some l,meN, 
then there is C < oo depending only on xi such that 

\\[Xi{.H),A]x2{,H)\\c(j)(^Hi'^),v{H^)) < C5. 
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b) Let (^H,V{H)^ be also self-adjoint on Ti.. If there are l,m E N with 
\\iH — X2(-^) 25(j|^m-i)) ^ ^! then there is C < oo depending only 

on Xi such that 

\\iXiiH) - XiiH))x2m\civim-.)MH-)) < C5. 

a) Let Ti. be another Hilbert space and B G C{T-C,T-C) such that BB* = 1 and 
[H := BHB*,V{H)) is self-adjoint on H. A ssume that there zs m G N such 
that B G C{V{H^),V{H^)) and B* G C{V{H^),V{H^)) for all I < m. 

i) Ifx2 e Co°°(M) and \\[H, B* B]x2{H)\\^^.^^^^^^^^ < S, then there is 

C <oo depending only on xi, X2, \\B\\c(v{hi),v{hi))> \\B*\\c{v{Hi),viHi)) 
for I < m such that 

\\{xiiH)-BxiiH)B*)BxliH)\\^^^^^^^^^^ < C 5. 

ii) If \\[H, B* B]\\c{viH^),v{H-^-^)) < 5, th^n there is C < 00 depending only 
on Xi, \\B\\c{v{Hi),v{Hi))' ^'^^ \\B*\\c{v{Hi),viH')) fori <m such that 

These statements can be generalized in many ways. Here we have given ver- 
sions which are sufficient for the situations that we encounter in the following. 
We emphasize that the support of X2 in a) and b) need not be compact, in 
particular X2 = 1 is allowed there. We now turn to the short derivation of 
Corollary [T] Lemma S] will be proved afterwards. 

Proof of Corollary [T] (Section 12.11 ): 

Let X '■ ~^ ["Ij 1] be a Borel function with suppx C (— oo,£']. Let 
:= min{inf cr(if^),inf cr(if|g)} and x,x e C^(M) with x\ie.,e] = 1 and 

xlsuppx = 1- The spectral calculus implies xi.H'^) = X^{H'^)x{H'^)- We 
recall from Theorem [1] that = U^H'^U'^* and = U'^*U'^. In view of 
Corollary \5[ satisfies the assumptions on B in Lemma H] c) i) with 5 = 
and c) ii) with 5 = e. Thus in the norm of C[L'^{C,dfilff),V{H^fi)) it holds 

x\h:^) = u^xmu'^'mi^) + {miK)-u'xmu'*Y 
+ {x{Hi^) - u^xmu'*) u^xmu'-* 

= WxmU'*x{Hl^) + 0{e') 

+ {x{Hl^) - U^xmU'*) U^~x\H')x{H')U'* 
= U^xmU^*x{H!s) + 0{e'). (37) 
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Then Corollary O yields that in the norm of C(TC,V{H'^)^ it holds 

U^*x{H:s)U' = P^x{H'W*x{Hl^W + 0{e') (38) 

because the spectral calculus also implies x(iy^g) = By 
Lemma [3] we have AA* = 1. Therefore 

(e"'^^* - A* U'*e-'^'^s'U'A) A* P'x{H') 

= (^(e-'^^*A* - A*e-'^'*) + A*{e-'"'' -W*e-'^'^«'U')^P'x{H') 

Since A and A* are bounded independent of e by Lemma [21 Theorem [T] 
implies that the second difference is of order e^\t\. So it suffices to estimate 
the first difference. The estimate (EH) implies [e"'^'*, P''] x(^^^) = 0{£^\t\) 
analogously with the proof of fl32l) . So 



[e-'^A^A* - A* e-'^'')P'x{H') 

= [A*P' - q'^a^A* P'e-'^'')xiH')A + A* [e"'^'*, P'] x{H') 

= ie-^^5i* f e'^A' {A*P'H' - H'^A*P' x{H')) e''^^' ds + 0{e^\t\) 
Jo 

® ie-'^^t f e'^A^A*H' -H'j^A*)P'x{H')e-'^^'ds + 0{e^\t\) 
Jo 

= 0{e'\t\) 

due to (136|) and \\x{H^)\\c{h,v{h^)) ~ 1- The latter holds because if'^ is 
bounded from below and the support of x is bounded from above, both 
independent of e. Because of ( 138|) and ||y4|| < 1 we have shown that 

||(e-i^^5.*_A*f/-e-^cVf/-A)A*[/-x(^eff)f^'^L(^2(^,^^^^ < C \tl 

which was the claim. □ 

Proof of Lemma HJ 

We want to apply the so called Helffer-Sjostrand formula (see e.g. chapter 2 
in [To]) to xi. It states that for any x ^ 

X{H) = - [ d,x{z)RH{z)dz, 

Jc 

where Rh{z) := {H — z)^^ denotes the resolvent and x : C ^ C is a so-called 
almost analytic extension of x- We emphasize that by dz we mean the usual 
volume measure on C. With z = x + iy a possible choice for x is 



Xix + iy) := riy) ^x^^\x 



3=0 ■' 
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with arbitrary r G C°°(]R, [0, 1]) satisfying t|[_i^i] = 1 and suppr C [—2,2] 
and I > 2. Then obviously x = x ^oi y = and there is < oo depending 
only on x such that 

d^xi^) := d,x + idyx < Cx\lmz\\ (39) 

which is the reason why it is called an almost analytic extension. We choose 
such an extension Xi ^ C^{C) of Xi with / = 2. Next we observe that for 
all J e No 



\R (.)\\ < v^l + 2|Im^P + 2kP 

\^H[Z)\\c{V{H^),ViH^+^)) < l4Uj 

because for all ip E H 

\\H^+^RH{z)^ljf + \\RHiz)^^ = \\HRHiz)W^p\\^ + \\RHiz)^^ 

< \\{1 + zRh{z))H^^\\^ + \\Rh{z)^^ 



|ImzP/ |ImzP 
|Im2;| 



Now by the Helffer-Sjostrand formula 

[Xi{H),A]x2{H) = - [ d,Xi{z)[RH{z),A]dzx2{H) 



chx.iz) RHiz)[A, H] Rh{z) dzx2{H) 
d^Xi{z) Rh{z)[A, H] X2{H) Rh{z) dz, 



where in the last step we used that [Rh{z),X2{H)] = due to the spectral 
theorem. Using the assumption ||[y4, i/] X2(-f^) [[^(^^(ho i'(-f/™)) ^ "we obtain 

\\[Xi{H),A]x2m\c^^^H 

< - \^Xl{z)\\\RH{z)\\ciViH"^),ViH"^+^)} 

X \\[H,A] X2iH)\\^^^^^i^^^^^^^^\\RHiz)\\c{v{H'-^),v{Hi)) dz 
mm ^ r ,^ ,2l + 2|Imz|2 + 2|^|2 

^ ^-V.P,J^"^' — — 

< C6, 
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with a C depending only on C^^ and the support of Xi- This shows a). The 
proof of b) can be carried out analogously because 

{xi{H)-xi{H))x2{H) = - [ d,Uz){RH{z)-Rfj{z))dzx2iH) 

Jc 

= - [ chxi{z)RH{z){H-H)Rfj{z)dzx2{H) 

Jc 

= - [ d,Xi{^)RH{z){H~H)x2iH)RH{z)dz. 

due to the Helffer-Sjostrand formula. For c) the formula yields: 

Xi{H) - Bxi{H)B* = - [ ckU^){R^{z)-BRH{z)B*)dz. (41) 

^ Jc 

So we have to estimate Rfj{z) — BRh{z)B* . We set A := B*B and note 
that BB* = 1 implies that BA = B, AB* = B* and = A. By definition 
H = BHB*. Therefore 

Rh{z) - BRh{z)B* = Rh{z){1-{BHB* -z)BRh{z)B*) 

= Rh{z) {1- B{H - z)ARh{z)B*) 

= Rh{z) {1- BAB* - B[H,A]Rh{z)B*) 

= -R^{z)B[H,A]Rh{z)B*. (42) 

For the first part of c) we compute 

B*{R^{z) - BRh{z)B*)B xliH) 

= -B*R^iz) B[H,A]RH{z)Ax2iH)x2iH) 

= -B*Rfj{z)B[H,A]RH{z){x2{H)A+[A,X2m)x2{H). (43) 

We will write Cb for a constant depending only on H-BH^^p^j^i-) ^^^5^;-)-, and 
\\^*\\civ{H''),v{H'')) f*^^ ^ — We note that the estimate fHOj) holds true with 
H replaced by H because H is assumed to be self-adjoint. Then, on the one 
hand, B e C{V{H"'-^),V{H"'~^)) implies 

\\R^{z)B[H,A]Rh{z) X2{H)Ax2{H)\\ 

= \\Rh{z)B[H,A]x2{H)Rh{z)Ax2{H)\\ 

Jl + 2\lmz\^ + 2\z\^ , , ,„ , , 1 

- 'Ih^l —\\{H,A]x2m\cin,viH--^)) llm^r' 



^ ^^l + 2\lmz\^ + 2\z\^ 

~ |ImzP 
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by the assumption on the commutator term. On the other hand, the assump- 
tions on B and B* imply that 



\Rf,{z) B[H, A]Rh{z) [A, X2{H)]X2{H)\\ 

= WR^iz) B{HB*B - B*BH)Rh{z) [A, X2{H)]x2{H)\\ 



2(1 + 2|Im^| 


^ + 2\z\ 






Im^; 


2 



2(l + 2|Imz| 


^ + 2\ 


\z\ 






\lm.z 


2 



< Cb rr:::::i2 \\[-A,x2{h)]x2{h)\\c('h,v{h^''-^)) 

where Cb,x2 depends also on X2 because in the last step we used that H,A, 
and X2 satisfy the assumptions of a). After plugging fH5]) into fl^Tl) the latter 
two estimates allow us to deduce the first part of c) analogously with a). 

For the second part of c) we observe that A"^ = A entails A[H, A] A = 0. 
Then we may derive from (1421) that 

Rh{z) - BRh{z)B* = -R^{z)BA[H,A]{1-A)Rh{z)AB* 

= -Rf,{z)BA[H,A]{l~A)[RH{z),A]B* 

= Rsiz)BA[H,A]RHiz)[H,A]RHiz)B* 

= Rfj{z)B[H,A]RH{z)[H,A]RH{z)B\ 

Hence, with B e C{V{H"'-^),V{H'^'^)) and B* e C{V{H"'-^),V{H'"'-^)) 
we obtain 

W^Hi^) ~ -^-^^^(^)-^*||£(D(i/m-i),X)(i?m)) 

= \\RHi^)B [H,A]RHiz) [H,A]RBiz)B*\\^^^^^^_,^^^^^^^^ 

2 



(1 + 2|Im;z| 


\^ + 2\z\ 


2)3/2 




Im^; 


3 



< Cb ' ' , ' ' ' ' ' \\[H,A]\\ 



< Cb6'^ IT I, — 

|lm2|'^ 

by assumption. Together with fHTl) this yields the claim as in a) when we 
put / = 3 in the choice of the almost analytic extension. □ 

Proof of Corollary O 

We will only prove that fl35l) is a consequence of the other statements. These 
follow directly from Lemma [2] by making use of the unitarity of Mp and Z)e 
as well as of D^{u)D* = {u/e), when we recall that P" = MpD^P^D*M* from 
the proof of Theorem [H 
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We prove fl5^ by induction. For m = both statements are clear. Now 
we assume that it is true for some fixed m G Nq. Theorem [T] yields that 
= U'*U^ and i/f^ = WH^W*. On the one hand, this implies 



Then \\P'^\\c{v{h^"^+'^)) ~ induction assumption immediately imply 

WU'^W r,^,zj.m+i. n^u£'"+l^^ 5, 1- Ou the othcr hand, we have 



= H'"' U'* His + H'"^[H',P']U'*. 
By the induction assumption and fl33p it holds for all ip that 

\\H'"'^'U'*ij\\ < \\H'''^U'*Hlsij\\ + \\H'''^[H',P']U'*^\\ 



where we used that lower powers of a self-adjoint operator are operator- 
bounded by higher powers. For e small enough, we can absorb the term with 
the e on the left-hand side, which yields ||f/^*|L,n,r7e'"+i\ -n^Lre'"+i\\ 5, 1- D 

' " ll£(D(_ff^^jj ),T>{H^ )) ~ 

3.3 Derivation of the effective Hamiltonian 

The goal of this section is to prove Theorem [2l We first take a closer look at 
the horizontal connection (see Definition [T]) : 

Lemma 5 It holds {V^^\ip)nf + {(l)\V^ip)nf = (d(0|V^)7^j) (r) and 

r,)^P := (v^^^ Vjl, - V^^.Vjl^ - Vf,^,,,])^ = -V^x(.,,.,).^, (44) 

where R"*" is the normal curvature mapping (defined in the appendix). 

The proof of this result can be found at the beginning of Section |H In order 
to deduce the formula for the effective Hamiltonian we need that ifg can be 
expanded with respect to the normal directions when operating on functions 
that decay fast enough. For this purpose we split up the integral over NC into 
an integral over the fibers NgC, isomorphic to M.^, followed by an integration 
over C, which is always possible for a measure of the form ® du (see e.g. 
chapter XVI, §4 of [21]). 
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Lemma 6 Let m G Nq. If a densely defined operator A satisfies 

\\Mi^y\\c(v(H^),n) < 1, \\{^yMc{viH^+^)MHe)) ~ ^ 

for every I E N, then the operators H^A^AH^ G C(V(H^~^^),'H) can be 
expanded in powers of e: 

H,A = [Ho + sHi + e^H2) A + 0{e^), 
AH, = A{Ho + eH, + e^H^) + 0{e^), 

where Hq, Hi, H2 are the operators associated with 

{(PlHo^U =11 g{eV''<P*,eV''^)dud^^ + {(PlHi^U, (45) 

JC JNqC 
JC JNqC 

{<P\H2^lj)n =11 3(7(>V.£VV,>V.£VV) +7^(£VV,^^,^VV,^^) 

JC JNqC 

+ |7^(£vV,^^,v>,z/) + |7^(w,^^,£vV,^^) 

where II is the second fundamental form, W is the Weingarten mapping, and 
TZ is the Riemann tensor (see the appendix for the definitions) . Furthermore, 
/or/ G {0,1,2} 

This will be proved in Section 14.21 Definition [3], Lemma [21 and the following 
lemma imply that Lemma E] can be applied to the projectors Pq and with 
m = 0. In the next lemma we gather some useful properties of the spectral 
projector Pq, the unitary 11, satisfying = U*PoUe from Lemma [21 and the 
global family of eigenfunctions ipf associated with Pq (see Definition [1]) : 

Lemma 7 It holds E{ G C^{C), as well as: 

t) Wl,j G No : \MPo{i^y\\cmH.)) < 1, \\[~e'A,„Po]\\civiH,m < e. 

ii) There are t/f,f/| G £(7^) H C{V{Hi;)) with norms bounded indepen- 
dently of e satisfying PqU^Pq = and f/|Po = Pof/|Po = Pof^l ■^'^^/i 
thatUe = l + eUl + e'^Ul 

Hi) \\PoUl{vy\\c{v{H^)) ^ 1 for all I G No and m G {0, 1}. 
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%v) For B, := PoU,x{He) and all u e {1, (f/f )*, (f/|)*} it holds 

II [-e^Ah + Ei,uPo] Be ||^(„) = 0{s). (47) 

t;; For RH^Ei) ■= (1 - Po)(^f - ^f)"\l - ^o) holds 

\\UrB, + Rh,{E,) ([-£Ah,Po] +^i)Po5.L(^,^(^^)) = 0{e) (48) 

and there is with sup^ ||e'^°^'^Vf('?)llwt(g) ~ '^'^^ 

sup,,cl|e'°<'^>V::^,..„,,V^^,...,,^y.f(g)||«,(,) < 1 
for all ui, ... ,1^1 G rb(A^C) and ri, . . . , G rb(TC). 



The proof of this lemma can be found in Section I4l3l Since f/| does only effect 
PeTi but not the effective Hamiltonian, we have not stated its particular form 
here, as we did for Uf in v). Now we are ready to derive the theorem about 
the form of the effective Hamiltonian. We deduce its corollary concerning 
the unitary groups before. 

Proof of Corollary [2] (Section 12.21 ): 
In order to check that 

(e-^^^* - t/o^*e-^-Vo^) f/o^*x(^^e?Vo^ < ^(l + ^'l^l), (49) 

with Uq = UqD*, indeed, follows from Theorem [1] and Theorem [2] we start 
by verifying that \\U^ - U^\\cCh,h,«) = ^(^)- 

We recall that we defined p := ^ as well as := M^UoUeD*^M* in the 
proof of Theorem [H Since dfi^^ is the volume measure associated with g^^, 
which is given by the expression in Theorem [2], we have ||p — 1 ||oo = C>{e) and 
thus II Mp — l\\c(L^{c,dfi)) = C^(^)- Using in addition that ||f4 — 1||£(h) = C^(^) 
by Lemma [7] and MiMp = 1 we obtain that 

\\U'-mci%H.,) = \\M;{UoUsD:m; - M,UoDt)\\cCH,n.,) 

= WUoUeD^M; - MpUoD:\\^^j^^^,^c4,^) 

= ||t/o/^:(M;-l)||^(^,^.(,^,^)) + Oie) 

= \\UoPoD:{M; - l)\\cCH,LHc4,)) + Oie) 

< ||(^)-'/^:(M,-l)||^(^,^) + 0{e) 
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because Uq = UqPq and the projector Pq associated with the constraint energy 
band Ef satisfies ||-Po(i^) ||£(W) ^ 1 by assumption (see Definition [3]). In view 
of (|T2|) . a first order Taylor expansion of p in normal directions yields that 
D*{M* — 1) is globally bounded by a constant times e{h'). Hence, we end up 
with \\U^ — f^o ll£(w,Wcff) ^ ^(^) ™-^y ^^^^ replace Uq by in (H9|) . 
Now let X • ["1; 1] be a Borel function with supp x C {—oo, E]. Using 
the triangle inequality and U'^U^* = 1 we see that 



-iHH 



< 



(2)^ 



+ 



-iHH 



cm 

{2)^r^£ 



eff- _ e 



C(H) 
C(H) ' 



(50) 



The second term is of order e^\t\ because 



P Gtl p Git 



(2)^ 
cfT / 



r(2) 



^(2)^ 



X(i/eff)e^^^«^ds = 0{eM) 



(2) 



by Theorem H Let x e 
and Lemma mb) we have 



with suppx|[ij^j:^(j:^(2)) = 1- By Theorem [2] 



S^) = u^*x{h'§)x{h'§) 



We recall from Theorem [T] that if^g = U'^H'^U'^* and P*^ = U'^*U'^ . In view of 
Corollary O f/^ satisfies the assumptions on B in Lemma Hie) ii) with 5 = e. 
Therefore 



u'*u'x{h')U'*x{h'§) + 0(^2) 



After plugging this into the first term in (1501) we may apply Theorem [T] to it. 
This yields the claim. □ 
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Proof of Theorem [2] (Section 12.21) : 

Let X : M ^ be a Borel function with suppx C {—oo,E]. We re- 

call that T>{A) always denotes the maximal domain of an operator A (i.e. 
all ip with WAipW + WipW < oo) equipped with the graph norm. A differ- 
ential operator A of order m will be called elliptic on {C,g), if it satisfies 
[ ... [A, /]..., /] > c|d/|™ for some c > and any /. 

m— times 

We set ifgg := —e'^Ac + Ei with Ac the Laplace-Beltrami operator on {C,g). 
Since E{ E C^{C) due to Lemma [TJ all powers of H^!^ are obviously elliptic 

operators of class C^{C) on Hcs- This implies that (-ffgff ' ^(-^is )) 
adjoint on Hcs because C is of bounded geometry (see Section 1.4. of [H]; 
in particular, this entails that 'D{H^) is the Sobolev space W^^'^(C), but 
equipped with an e-dependent norm). Let E^ := min{inf a{H'^), inf cr{H^'^)} 
and XjX ^ C!^(R.) with x\[e.,e] = 1 and xlsuppx = 1- Then we define if^g 
for 0,V^GP(if^°))by 

{<I>\h!:^^P) := ^(^?:ff((p^fj0)*,p^fjV^) +0*(Ef + e(y,f|(v:'l^VfK)^ 

+ _ £2^($*(0),$(^(i7W)^)) (51) 

-e'M{^*mHi^)cj^), $(^ - ^(ifSV))) d^il, 

where $(V') := ^(^Vp^gV'jPeff'^; V') a^id all the other objects are defined by 
the expressions in Theorem El Because of ^{H^'^)x{H^s) = x{H^s) this 
definition immediately implies that H^^ operates on ip with ip = x{H^)il) 
as stated in the theorem. 

The rest of the proof will be devided into several steps. 
Step 1: {H^q , T'(iJ^g )) is self-adjoint on Hcs and 

W^is -Hi^\\ciV(Hi°J),H.s) = ^(^)- 

It easy to verify that H^g is symmetric. Then it suffices to prove the stated 
estimate because by the Kato-Rellich theorem (see e.g. |37j) the estimate 
implies that ((-f^es ' ^(-^iff )) self-adjoint on Hes for e small enough. 
Since C is of bounded geometry, maximal regularity estimates hold true there 
(see Appendix 1 of [H]), in particular, differential operators of order m G N 
with coefficients in C^{C) are bounded by elliptic operators of same order 
and class. 
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The operator M associated with ($(0), is a fourth order 

differential operator which, in view of Lemma [7] vi), has coefficients in C^{C). 
Hence, it is bounded by [H^^Y with a constant independent oi e because all 
derivatives carry an e. We notice that [[^(ifg) || -n^H■(o)™^^ ~ ^ m G Nq 

because the support of x is bounded independently of e. Thus we obtain that 
M^{H^'^) is bounded. The same is true for ^{H^'^)M{1 - ^{H^'^)) because 
it is operator-bounded by the adjoint of MxiH^'^). Therefore the A^-terms 
in (15T|) correspond to bounded operators! All the other terms are associated 
with differential operators of second order whose coeffcients are in C^{C) by 
Lemma [7] vi) and whose derivatives carry at least one e each. Therefore they 
are bounded by the elliptic H^g. 

So we obtain that \\H^^ — H^s\\c(v{h^°^)h g) ~ ^(^) observing that the 
leading order of H^^ is indeed H^q. 

Step 2: V{Hl^) = V{H^^) and \\Hl^ - H'§\\c^j,^Hi,),n.,) = 0{e). 

We recall that we defined Ue := UoUe, := M? UeD*M;, := U^*U', and 
:= V^H^U^* in the proof of Theorem [H which implied = U*Uf.. 

Since \\Ue\\c{v{H,)) < 1 and l|f^il£(25(//^),i5(i^(o))) < 1 by LemmalU it also holds 
\\^''\\c{v(H-),v{H^^)) ~ ^- Using, in addition, that \\U''*\\c(v(h-),v(hi^)) < 1 due 
to Corollary [5] and U'^U'^* = 1 we conclude that for all ip G V{H^fi) 



On the other hand. Lemma [7] and Corollary [S] imply via the analogous argu- 
ments that for all ip G ^'(iJgg ) 

Hence, = 



Using Hlff = U^H^U^* = WP^H^P^W* and again Corollary we get 

(0), 



\\r n r- u n^^u \\c{v{h^),h) 
= \\P,H,P,-U:MpHl^M;U,\\civ(H,),n) 
= \\P,H,Po~U*MpH^SMW4cmH.),n) + 0{e) 
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because = U*Us, = UqUs, and by Lemma [7] ii) it holds f/^ — 1 = 0{e) 
both in C{n) and in C{V{H^)). Lemma Elimphes that Po{H,-Hq)Pq = 
in C{V{He),n). Hence, 

= \\PoH,P,-U;M-,Hf^M;Uo\\cmH^),H) + 0{e) 

< \\U,H,U; - M,i/S^M;||,(^(^,o))^,,(,_,^„ + (52) 

where in the last step we used Pq = UqUq and ||?7o|| /■/'•nm ^-n,'H■(o)^^ ~ 1 due 
to Lemma [7] vi). It holds Uoip = ipfip by definiton of Uo and Hq = — Ah + -fff 
by Lemma [61 In view of Definition [H we have 

eV^ipifii = (fiiedip + V^eVVf, (53) 

where d is the exterior derivative on C. We note that sup^ HeV'^v^f ||-Hj.(g) and 
supg ||e^Ah</?f ||>if(g) are of order e and e"^ respectively by Lemma[71 Therefore 

UoHoU*^ = Uo{~e'A^ + Hf)U*^ = {^i\{-e'A^ + Ei)<fi^)n, 

= HjS^ + Oie). (54) 

We recall that p = dpl^/dp with dpl^ the measure associated with g^f^. Since 
dfi and dpl^ coincide at leading order and p G C^{C, g) due to Lemma [7] vi), 
we have \\M*~ - '^W^jj^jjW^) = 0{e) and \\Mp - iWciL^cdf^)) = C>{e). So we 

obtain that ||Mpi7^(°^Mi - /f^J^ ||^(^(^(o))_^2(c,d^)) = ^i^)- Together with §^ 
and m this yields - f^iJ = Ois). 

Step 3: It holds \\{H^^s - ^j:^) U'x{H^)U'*\\cin^^) = 0{e^). 

This step contains the central order-by-order calculation of H^g and is there- 
fore by far the longest one. For any ip we set ip := Mpip, ip^ := U'^x{H'^)U'^*ip, 
and 'ip^ := UeX{H^)U*i). Of course, we have = 'ijj^, ||^||L2(c,d/.) = ll^lkeff' 
and ||V'xllL2(c,dM) ^ Ml^ca^) for all V e Kff- 

We first explain why the cut off in the definition of H^^ does not matter here. 
We note that and f/^ satisfy the assumption on A and B in Lemma H] a) 
and c) ii) with 5 = e by Corollary O In addition, Hj:^ and H^g satisfy the 
assumption of Lemma H] b) with the same 6 by Step [2j Therefore 

^{Hi^)U^xmU'* = mi^W'xmU'* + 0{e) 

= W^{H')P'x{H'W* + 0{e) 

= U^P'^iH^xmU'* + 0{e) 

= U'xiHnU'* + 0{e), 
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which shows that 



+ 0{e'mnjmU,). (55) 

So now we aim at showing that the same is true for ( | In the fol- 

lowing, we omit the e-scripts of Uf, [/|, and Us and set Hh '■= L'^{C, dfi). 
Next we will show that 

{4>\H,s^P^)h.^ = {4>\Uo{Ho + eH^ + e'H,)U*tP^)n, 

+ £ ( I t/o (f/i (Ho + eH,) + {Ho + sH,) U^) U* ^ 

+ 6'{4>\Uo{UiHoU* + U2H0+ HoU;)U*4>^)n, 
+ OiE'\\<l>\\njmn.,)- (56) 

By definition of Hes it holds 

If we could just count the number of e's after plugging in the expansion of 
Hs from Lemma M and the one of U from Lemma [TJ the claim (1561) would be 
clear. But the expansion of H,. yields polynomially growing coefficients. So 
we have to use carefully the estimate fH6l) . 

By Lemma[7]it holds \\uPoJJ\\c(v(h,)) ^ 1 for each u e {U*, 1, U^, U^}. Since 
U*Po = PeU*Po and U^Pq = PqU^Pq by Lemma El uPqU satisfies the as- 
sumptions on A in Lemma [6] with m = for all those u due to the decay 
properties of P^, Pq, and U*Po from Lemma [2] and Lemma [71 We notice that 
\\x{He)\\ c{'H,v{H^)) ^ 1 because is bounded from below and the support of 
X is bounded from above, both independently of e. Hence, using UqU,. = PqU 
we may conclude from (|3H|) that 

WhuU^i^Jn = \\huPoUxmu:4,\\n < UWh., (57) 

for each h G {H^, Hq, Hi, H2}. Furthermore, Lemma [HI implies in the same 
way that 

\\{H, - {Ho + sHi + e^H2))U*U*^^\\^^^ = 0{e'). 
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So we have 

= {H^ + eHi + e^H2){l + eUl + e'u;)U*i,^ + ^(^l^ll) 

+ e{Ho + eHi)U* + e^HoU;^U*4^^ + 0{e^ij\\). 

For the rest of the proof we write 0{e') for bounded by £^'||0||w^g||^||7i^ff times 
a constant independent of e. The above yields 

{(j)\H,sij) = {4>\UoUH,U*U*i^^) 

= (01 UoUiHo + sH^ + e'^H^) U* ^ ) 
+ e{^\UoUiHo + sH,)U;u;^^) 
+ e^{U*4>\UHoU2U*i^^) + 0{e^), 

After plugging U = l+£Ui+e^U2 we may drop the terms with three or more 
e's in it because of (157|) . Gathering all the remaining terms we, indeed, end 
up with (l56i) . 

Now we calculate all the terms in fl56l) separately. By Definition H] 

(0|f/oAf/*^;,)w, = {ip,^\A^,i,^)n. (58) 

for any operator A. Furthermore, the exponential decay of and its deriva- 
tives due to the Lemma [7] guarantees that, in the following, all the fiber 
integrals are bounded in spite of the terms growing polynomially in v. 
We observe that ifj^ = UoUx{,He)U*ilj implies that 

ll^eff < WH^xmWMn + WxmWMn < 1 

because l|f^o||£(2,(j^^)^x,(jy(o)) ^ 1 and \\U\\c(v{h,)) < 1 by Lemma [71 As ex- 
plained in Step [1] every differential operator of second order with coefficients 
in C^{C) on Ties is operator-bounded by H^'^ . Therefore derivatives that 
hit ipx pose any problem, either. These facts will be used throughout 

the computations below. We write down the calculations via quadratic forms 
for the sake of readability. However, one should think of all the operators 
applied to (p as the adjoint applied to the corresponding term containing ip. 
Since Hv^f ||?Yf(g) = 1 for all q EC, Lemma [5] implies 
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Thus {ip{\V^{pf)T-j^^ = Im(y3f I V'^v^f)>^f . Therefore the product rule imphes 

{(pi(j)\Ho(piip^)n 

m 



Jc Jc JNqC 

I 4>*Eii>^dfi + I I |(/^fp(?(£d0*,£d4) +£(7(¥^f*£d0*,4vVf) 
Jc Jc JNqC 

+ £(?(0*VVf,¥'f^dV^x) + £'^7(0*VVf,V^xVVf) dudfi 
fi'((Pefr0)*,Pefr^x) + 0*^fV^x + £^^0*Vbh^x ^/^ 



5^ 



with 



p^ffT/; = -iedip - Im (e((^f|VVf)wf - / | <^f R(V>f, z/)z/ dz/ 

+ 52(v^f|2(W(.) - (^f|W(.)^f)7.,) VVf>^JV', 

as well as ri := ImRi for i?i := ( | 2(W( . ) " ( '/'f I W( . )(/?f )n,) VVf >7^^ 
and r2 := Imi?2 for R2 := /^y | R(V^v^f, z/)z/(ii/. When we split up Ri 
into real and imaginary part for i G {1, 2}, an integration by parts shows 



g{{-ied(j)y,ilj^R^) +g{(j)*R*, -iedip^) 

g[{-ied(i))\i)xri) + g{(i)*r*,-ied^x) c//i + 0{£). 
Therefore the ri-terms are cancelled by terms coming from Hi. 

# / / 2II(z/)(eVVf0*,^VVfV^x) + 0*(V:H^)|^fNxC^i^t^/i 

Jc JNqC 

= [ [ |<^fp2II(z/)(£dr,£dV'x) + £2II(z/)(^*£d0*,4vVf) 

Jc JNqC 

+ £2II(z/)(0* VVf,</'f^d4) + (f>* (yiW)\^i\^ du dfi + Oie^) 
(v9f|2II(.)((Peff0)*,Peff4)<^f)HfCi/x + J^^*{ipt\iV':W)ifi)nJxdi^ 

+ e J^g{{-ied4>yJxRi) + g{4>* R*i, -iedi^^) dfi + 0{e^), (60) 

50 



where we used that 5'(ti, >V(z/)r2) = II(z/)(ri, = 5'(W(z/)ri, (see the 
second appendix). At second order we first omit all the terms involving the 
Riemann tensor: 

{ificj) I H2 (p{ip^)n — 'Riemann-terms' 

# / / 3^?(W(z/)£VVf*0*,W(z/)£VVf4) 

JC J NqC 

ipi\3g{W{ . )(Pefr^/'x)*' ■ )Pcsi'x) 'Pi)^ 

+ jj* ((^f|(iV?;.W^)^f)w, + rgeom)V^x + (61) 

where we used that —iedip^ = Pesi^x + ^(^) ^^^^ step. Now we take 

care of the omitted second order terms. Noticing that V'xV^f ~ "^x ^^V'f 
we have 

' Riemann-terms ' 

m 




= [ [ |</^f|'^(£d0*,z/,£d4,z/) + |^((/^f*£d0*,i/,4VVf,z/) 

+1 ^(0*V^(^f*, z/, (^f edV-;^, I/) + i 0*^(V^<^f*, u, V>f , I/) V-;^ di/ rf/i + 0{e) 
= J {ipf\n{ed4>*, . ,edijj^, .)(f{)^^dfi + J 4>*V^,nhi'xdf^ 

+ J^g{{-ied4>y,^xR2) + 9{4>*R*2, -i^d^) rf/i + 0{e) (62) 

with Vamb = In cl '^(^^V^f ' ^^V^f) Again replacing — ied with 

and (? with Qcs yields errors of order e only. In view of (l58l) - (l62l) . we have 

(0 I Uo {Ho + eH, + e^H2) U* ^^)n, 

9ts ( (Peff0) * , Pea^x) + 

+ 0* {e{^f\V':W^f)Ht + e'VT^'O V^x dfi + 0{e^) (63) 
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with 



g{n,T2) + e ((^f |2II(.)(n,T-2)<^f)wj 



Hi 



We define Pf^ := (1 — Pq)- Before we deal witli tlie corrections by Ui and U2 
in ([56]), we notice tfiat due to Pq = ^0^0 PttU^ = 

Po^([-£Ah,Po] + //i)f/o>x 



= Po^[{VlW) - trc(2(VVf)f/o + £V^^W(z/)) ^V^j U*^^ + 0{e) 

= Po^(<^f(V:iy)V'x - 2^7(VVf,^d^x) - '^ftrcW(z/)£2vd^^) + 0{e) 
= Po^v]>(eVd4,d^^,7Aj + (64) 

with = -v9ftrc(W(z/)A) - 2(7,^jj(VVf*,p) + <^f(V^W^)0. 

We note that U^iIj^ = B'U*i) with = PoUxiH,). So we may apply dS]) 
und fHHj) in the following. Since Uo = UqPo by definition and we know from 
Lemma[7]that PqUiPq = 0, the first corrections by Ui are an order of e higher 
than expected: 



Uo (Ho + sHi) U* + Ui (Ho + sHi) U* 



uJ{[Po,Ho]+eH,) U* + U, {[Ho, Po] + eH,))U* i> 



t/o(([£Ah,Po] + f/i {[-eA^,Po]+H,)Po)u*i>. 



Uo ([^Ah, Po] + i/i) P/^,(Pf) ([-5Ah, Po] + i/i) U* 



-£(0|f/of/i(//f-Pf)f/;f/o*^^)Hb 



Hb 



-5(01 Wl(i^f-i5^f)f/l*f/oV^x)Hb 

-e yA^(^*(£Vd0,ed0,0),^(£2vd^^,ed^^,V^;^))) rf/x 
- £ (0 I t/ot/i {Hi - E,) U*U* ij^)n,. 



(65) 
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with M{^1, ^2) = ( </^i I (1 - Po) (^f - ^f) (1 - Po) ^2 >„j. Furthermore, 

{4>\Uo{U2Ho + HoU;)u;i^^)n, 

= {Uf4> I Po{U2 (-e'Ah + H,) + (-£2 Ah + H,) U;) PoU*Mn, 
= {U*4> I {P0U2 (-e'Ah + Ei)Po + Po(-^r'Ah + ^f) U;Po) U*i^^)n^ 

^ {U*4> I Po (f/2 + f/2*)Po(-£'Ah + E() U*tP^)n, + 0{e) 
= - I Po UiU*Poi-6^Ai, + El) UMn, + 0{e), (66) 

because U = l + eUi+ e^U2 imphes via PqUU* Pq = Pq and PqUiPq = that 
Po(^2 + ^^2)-Po = -Pof^it^iPo + C){e). Finally, the remaining second order 
term cancels the term from (l66l) and the second term from (l65l) : 

(0|f/o?7ii/of/rf/o4>„, 

= (0 I f/of/i (iff - ^f) UlU; 4 + f/of/i(-£'Ah + Ei) U^PoU* ^^)h^ 

© (0if/ot/i(ijf-Pf)t/rf/o*4)7Y, 

+ (0 I f/of/if/rPo (-e'Ah + ^f) U* i^^)H, + 0{e). (67) 

We gather the terms from fl63|) to fl67|) and replace d-^^ by Peff'^x 
argument of which only yields an error of order e^. Then we obtain 
that ( (p I ii|g ) equals the right-hand side of flSSl) up to errors of order e, 
only with d/i instead of dfies- Here ip = Mpip enters. By Lemma [T] Mp 
interchanges the former with the latter but may add extra terms. However, g 
and Qes coincide at leading order and so do their associated volume measures. 
Therefore d(ln p) and Ac In p are of order e. This shows that the extra 
potential from Lemma [H given by — ^ (7(d(lnp), d(lnp)) + YAc(lnp), is of 
order e^. Exploiting d(lnp) = 0{e) we easily obtain that all the other extra 
terms are also only of order e^, which finishes the proof of Step [31 

Step 4: It holds WiHI^ - H^;S)x{Hl^)\\c(H..) = 0{e^). 

The spectral calculus implies xiHcs) — 'X^{Hls)x{.Hls)- i'^ (|37|) this 
implies that 

\\x{Hl^)-U^miU^*x{Hl^)\\cin.,MH^^,)) = 0{e^). 
Now Step a follows from Step [3] and V^H^j^) = V{Hl^) due to StepE 
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step 5: It holds - H^S)x{H^S)\Wh..) = 0{e^). 

We note that Step EfcE] imply that II - i/^g^ II £(o(H-g),Hefi) = C'(e). So 
in the norm of ^(T-^efr, ^^(-f^cfr)) holds that 

+ [xiH'S) - X{HI^)Y X{HI^) 

+ (x(^S^) - x{Hi^) x{Hi^) miS) 

by Lemma H] b) and Step [2] & HI Hence, Step [5] can be reduced to Step H] in 
the same way as we reduced Step S] to Step [31 

Theorem [2] is entailed by Step [3] to [5] and the remark preceding Step [1] □ 

3.4 Proof of the approximation of eigenvalues 

With Theorem [21 Corollary [5l and Lemma [H we have already everything at 
hand we need to prove Theorem [31 which relates the spectra of H'^ and Hj:^ . 

Proof of Theorem [31 (Section 12.31) : 

We fix E < oo and set E_ := minjinf cr(if^), inf a{H^^)} — 1. Let x be the 
characteristic function of [E_,E] and x ^ C*^(]R) with x\[e_,e] = 1- 

(2) 

To show a) i) we assume we are given a family of eigenvalues (Es) of H^g 
with limsupE'e < E and a corresponding family of eigenf unctions (ipe)- 
Since ipi; is an eigenfunction of -f^^g , we have that tjj^ = x{H^g)ip£ for 
e small enough. By Theorem [2] and Lemma [H b) it holds in the norm 

of /:(L2(C,c//ieff),I^(i/eff)) 

xiHi?) = x\hS)x{h'S) 

= x\HI^)x{h'S) + 0{e') 

© U^x{H''W'x{Hl^)x{H^:S) + 0{e^) 
= U^xmU'*x{H!:l^) + OiE'). (68) 
Therefore with = P^U^*, U'*U' = P^ and i/^jj = U'H'U'* 
U^'*^s = (P^ + (1 - n)H^'W*x{H!§)^e 

= U^*H!sX{H^^)i^e + (1 - n[H\ PI U^*x{Hi^)A 

# t/^*i/effV. + (1 - niH', n xm w*x{H!S)i^e 
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where we made use of the assumption and Corollary [S] in the last step. This 
proves a) i) because U'^U^* = 1 and thus II^AellTi^g- = ||f/^*^/'e||^. 
To show a) ii) we now assume that we are given a family of eigenvalues [E^) 
of H'^ with limsup < E and a corresponding family of eigenf unctions {ip'^). 
Here this implies ip"^ = x{H'^)ip'^ for e small enough. With = U^P'^ and 
lje*ue ^ pe obtain 

where we used Lemma H] a) & c) in the two last steps. In view of Theorem [2], 
we get 

Hi^U^i;^ = Hl^x{Hl^)U^x{H^W + 0{e') 

= U'H'U'*U'x{H')P'x{H')^' + 0{e^). 

Using again Lemma H] a) & c) and the assumption we end up with 

H^^U'^' = U' x{H'W + 0{e^) = H'x{H')^' + 0{e^) 

= E.U'il)' + 0{e^). 

This finishes the proof of a) ii). 

For b) we set ipe := D*M*ip'^ and observe that — e^A^ = D^A^D* by Defini- 
tion □ and thus -e^MpA^M* + Voiq,u/e) = MpD,HiD*M*. Therefore the 
statement is equivalent to 



limsup (v^,i//f^,) < infEiiiv^.r ^ wuei^eW > \m 

because := MiUeD*M* by definition in the proof of Theorem[TJ We have 



> m^Ei\\Po^Per + inf^ill(l-Po)^, 



= inf llV^.f + (inf El - inf Ef) ||(1 - Po)^ef 
qec qec qec 

= infEfll^.f + (infill -infi^f) 11(1- P,)^,f + 0{e), 

qGC q&C qGC 

where we used that — Pq = C'(e) by Lemma [2] in the last step. Since 
E{ is a constraint energy band, hence, separated by a gap from Ei, and 
limsup ('?/'£ I iff ■j/'e) < infggc -^illV'elP by assumption, we may conclude that 

limsup 11(1 - Pe)i'e\\'^ < limsup \\ipeV- 

Because of Pe = U*Us this implies > \\ipe\\ for all e small enough. □ 
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4 The whole story 



In Section [3] we proved our main theorems with the help of Lemmas [T] to HI 
We still have to derive Lemmas [2] to [71 which is the task of this section. Before 
we can start with it, we have to carry out some technical preliminaries. 

Remark 6 Since C is of hounded geometry, it has a countable covering 
of finite multiplicity (i.e. there is Iq such that each Qj overlaps with not 
more than Iq of the others ) by contractable geodesic balls of fixed diameter, 
and there is a corresponding partition of unity (^j G C^{Qj)) whose deriva- 
tives of any order are bounded uniformly in j (see e.g. App. 1 of fJJ^ ). 
We fix j ^ By geodesic coordinates with respect to the center q G Vtj we 
mean to choose an orthonormal basis {vi)i ofTqC and to use the exponential 
mapping as a chart on flj. Let (a;*)j=i,,,,,rf be geodesic coordinates on Qj. The 
bounded geometry of C that we assumed in ^ yields bounds uniform in j 
on the metric tensor gu and its partial derivatives, thus, in particular, on all 
the inner curvatures of C and their partial derivatives. For the same reason 
the inverse of the metric tensor gu is positive definite with a constant greater 
than zero uniform in j . 

We choose an orthonormal basis of the normal space at the center ofQj and 
extend it radially to A^CI^^ = NQj via the parallel transport by the normal 
connection V"*" (defined in the appendix). In this way we obtain an orthonor- 
mal trivializing frame (z/q,)q, overQj. Let {n°')a=i,...,k be bundle coordinates 
with respect to this frame. The connection coefficients T]^ of the normal 
connection are given by Vg^ Va = Ylj=i ^la^-r- ^'^^ the smooth embedding 
ofC assumed in ^ the exterior curvatures ofC, the curvature of NC, as well 
as all their derivatives are globally bounded. This implies that all the partial 
derivatives of Fj^ and of the exterior curvatures of C are bounded uniformly 
in j in the coordinates {x^)i=i^...^d (^nd {n°')a=i,...,k- 

From now on we implicitly sum over repeated indices. The vertical derivative 
in local coordinates is given by 

{V:j){x,n) = 9„Xx,n). (69) 

and the horizontal connection is given by 

{V\^i:){x,n) = d^.ij{x,n) - Tln-dnM^,n). (70) 

The former directly follows from the definition of (see Definition [T]). To 
obtain the latter equation we note first that for a normal vector field v = rfVa 
over C it holds 

{Vivy = d.M^ + F7„n". (71) 
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Now let {w,v) e Ci([-1,1],A^^) with 

w(0) = X, w{Q) = d^^, k w(0) = n, V^v = 0. 
Then by definition of V*^ we have 

= dxiip{x,n) + {d,j.in'^)dn-t4'{x,n) 
= ^^^^jJ{x,n) - T]^n°' dni^ix,!!.), 

where we used ( ITT]) and the choice of the curve v in the last step. 

With the formulas fl69|) and flTOj) it is easy to derive the properties of V'' that 
were stated in Lemma O 

Proof of Lemma [5] (Section 13.31) : 

Let T,Ti,r2 E T(TC) and tjj,tjji,tp2 G C^(C,7if(g)). We fix a geodesic ball 
Q E C and choose {x'^)i=i^...^d and {n°')a=i,...,k as above. We first verify that 

is metric, i.e. (d (V'll^sJwJ (r) = {\/'^ii\^2)H, + {H'^^Mur Since 
is a metric connection, T]^ is anti-symmetric in a and 7, in particular = 
for all a. Therefore an integration by parts yields that 

Therefore we have 

(d(V'i|V^2))(r) = r^(5.,V^i|V^2) + r^(z^i|9.,^2) 

= (V^V^i|V'2) + (^i|V!l^2). 
To compute the curvature we notice that 

Using the commutator identity 
we obtain that 

which was the claim. □ 
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4.1 Elliptic estimates for the Sasaki metric 

In the following, we deduce important properties of differential operators 
related to the Sasaki metric defined in the introduction (see ([9])), in particular 
we will provide a-priori estimates for the associated Laplacian. 

In bundle coordinates the Sasaki metric has a simple form. Here we keep 
the convention that it is summed over repeated indices and write a^^ for the 
inverse of aij. 

Proposition 1 Let be the Sasaki metric on NC defined in Choose 
Q G C where the normal bundle NC is trivializable and an orthonormal frame 
(z/q,)q, of NC\q. Define T]^ by Vg^ Ua = ^la^-y- ^''^ corresponding bundle 
coordinates the dual metric tensor gs € Tl{TNC) for all q & Q is given by: 

= {c^ i) (o b) (o i*) ' 

where for i,j = 1, d and a, 7, 5 = 1, .., k 

A'\q,n) = g'^iq), B^\q,n) = 
Cj{q,n) = -n"C(g). 

In particular, {det{gs)ab){(l,n) = {det gij){q) for a,b = 1, ...,d + k. 

The proof was carried out by Wittich in [15] . From this expression we deduce 
the form of the associated Laplacian. 

Corollary 6 The Laplace- Beltrami operator associated with gs is 

As = Ah + A,. 

Proof of Corollary [S 

We set fj, := det gij and fis '■= det{gs)ab- Since {1^0)^=1 orthonormal 
frame, we have that g(qfi){h'a, up) = 6"'^. So fl69|) and fITOl) imply that 

A, = 5„.(5"^9„, & Ah = fi-' {d,. - dn.)fig'' [d,, - T]^ . (72) 

Now plugging the expression for g^'' and det g^^ from Proposition [1] into the 
general formula As = Yltj>=iifJ-s)~^da fJ.s ds'db yields the claim. □ 

Next we gather some useful properties of A^, Ah, and V*^. We recall that in 
Definition [2] we introduced the unitary operator D,. for the isotropic dilation 
of the fibers with e. 
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Lemma 8 Let f e C^iR) and r G T{TC) be arbitrary. Fix A G M. It holds 
%) DA.Dl = £2 A,, D,/\^Dl = Ah, D,V,Dl = V , 
ii) [V^A.] = 0, [Ah, A.] = 0, [V^/((Az/))] = 0, 
m; [A.,/((Az.))] = Af((Az.))(A ^<^-/;-riV p|yVg + AV"((Az.))i^ . 

In the following, we write A ^ B when A is operator-bounded by B with a 
constant independent of e, i.e. if V{B) C and < + ||^|| 

for all tl) G V>{B\ We will have to estimate multiple applications of 
and V*^ by powers of if^, which was defined as if^ := DlMpH'^MpD^ with 
:= —e^/S.Nc + V'^ . Essential for our analysis, especially for the proofs of 
Lemmas [2] & [71 are the following statements: 

Lemma 9 Fix m G No and M G {0, 1,2}. For all I e Z, X e [0, 1] and 
nil + ^2 < 2m the following operator estimates hold true on Ti: 

i) -< (-e^Ah-Av + V,)™ -< H^, 
n) (-A.)'" (-5^ Ah)*' -< 

Hi) {XuY [H^^^^ , {Xu)^'-] -< H^^^^ with a constant independent of X, 

The last three estimates rely on the following estimates in local coordinates. 
Here we a use covering of C and coordinates {x'^)i=i^...^d and {n°')a=i^...^k 
as in Remark E] in the introduction to Section |H 

Lemma 10 Let a,(3,'j be multi-indices with \a\ < 21, \a\ + < 2m and 
\j\ = 2. Set fj, := det gij. For all smooth and compactly supported it holds 

< IK-e^Ah-A. + K)"^!! + ll^ll, 

< iK-^^^h-^^A. + i^O^^II + ll^ll- 
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We now provide the proofs of these three technical lemmas. 
Proof of Lemma [SJ 

We fix a geodesic ball Q G C. Let {1^0)0=1,.. .,k be an orthonormal trivial- 
izing frame of NQ with associated coordinates {n°')a=i,...,k and {x^)i=i,...,d 
be any coordinates on f2. Observing that Dsip{x,n) = e~^^'^il){x^n/ e) and 
D*ip(x,n) = e''^'^ip{x,en) we immediately obtain i) due to (!72l) . 

Since V"*" is a metric connection, T]^ is anti-symmetric in a and 7 and so 
(Trnjl implies 

Va^^^(g, z/) = ^^^^p{x,n) - ^T]^{n°'dni - n'^dnc')ip{x,n). 

Using that = 5"^dn'^dnP by (!72|) we obtain that for any r = r^d^i 

[V^A.] = rT7"(a„.9„. = 0. 

We recall that (z/) = ^/l + g(q,o){i^, z^)- Since is an orthonormal frame, 

we have that g(^qfi){iya, up) = 5°-^. This entails that {v) = ^Jl + dapW^n^- 
With this the remaining statements follow by direct computation. □ 

Proof of Lemma M 

We recall from Definition [2] that Ve = Vc + D*WDs and that we assumed 
that Vc and W are in C^{C,C^{NgC)). These facts together imply that 
VeeC^{C,C^{N,C)). 

Since and Mp are unitary. Lemma [S]i) yields that Lemma Oi) is equivalent 
to 

{HT ^ Mp(-£2Ah-£'A, + l^^)"M; ^ {HT (73) 

for all m G N. By choice of 'g it coincides with the Sasaki metric outside 
of Bs and, hence, so do A^c and As. In addition, this means p = 1 outside 
of Bs and so Mp is multiplication by 1 there. Then Corollary [6] implies 
= Mp{- £2 Ah - e^A^ + V^)M* on NC \ B5. Hence, by introducing 
suitable cutoff functions it suffices to prove (173|) for functions with support 
in B25 n NQj. The set B25 H NQj is easily seen to be bounded with respect to 
both g and g^ and thus relatively compact because NC is complete with both 
g and g^ as explained in the sequel to the definition of in ([9]). Furthermore, 
on B2S n N^j both {H')"' and Mp{ - e^A^ - e^A^ + V')"^ M; are elliptic 
operators with bounded coefficients of order 2m. Therefore fl73l) follows from 
the usual elliptic estimates. These are uniform in j because B25 is a subset 
of bounded geometry of NC with respect to both 'g and g^, which was also 
explained in the sequel to 

In the following, we prove the estimates only on smooth and compactly sup- 
ported functions, where we may apply Lemma fTOl Then it is just a matter of 
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standard approximation arguments to extend them to the maximal domains 
of the operators on the right hand side of each estimate. In this context one 
should note that the mamixal domains V(H^) and V{{—e^Ah — Av + Vs)"^) 
coincide for all m G N by i). 

We recall that K e C^{C,C^{NgC)) and turn to ii). By i) we may replace 
He by -e^Ah - A^ + V;. We first prove the statement for M = inductively 
in m. In view of (1721) . Lemma fTOl ii) implies that — A^ -< —e^A^ — Av and 
thus also — e^Ah -< — e^Ah — Av. So due to the boundedness of Ve the triangle 
inequality yields the statement for m = as well as 

-^^Ah ^ -^^Ah-Av + K- (74) 



In the following, we will write A ^ B + C, ii \\Aip\\ < \\BiIj\\ + WCipW + 
We note that with this notation A ~< B implies AC -< BC -(- C. 
Now we assume that the statement is true for some m G Nq. Since V'^ G 
and NC with the Sasaki metric is complete, the operator —e^As + K 
is self-adjoint on Ti and so is — e^A^ — Av + as it is unitary equivalent 
to -e'^As + V via D^. Therefore by the spectral calculus lower powers of 
— e^Ah — Av + Ve are operator-bounded by higher powers. In addition, Av 
and Ah commute by Lemma [HI Then we obtain the statement for m + 1 via 

(_Av)'"+i -< (-£2 Ah - Av + Ve) (-Av)'" -f (-Av)'" 

= (-Av)'"(-£2Ah - Av + Ve) + [Ve, (-Av)™] + (-Av)'" 
-< (-^^Ah - Av + v;)'"+i + (-Av)"* 

.h - Av + Ve) 



-< (-^^Ah - Av + v;^'"+' 



Here we used K G C^{C,C^{NqC)), A^ = S'^^dr.cd^p locally, and i) of 
Lemmamto bound [K, (-Av)""] by (-Av)"". Using [Av, AJ = and ([ZID 
we have 

(-Av)'" (-^^Ah) = (-e^Ah) (-Av)"' 

-< (-e^Ah - Av + V) (-Av)"' + (-Av)"'. 

Continuing as before we obtain the claim for M = 1. Furthermore, 

(-Av)"'(-5^Ah)^ = (-£2Ah)(-Av)'"(-e^Ah) 

-< (-£2Ah-Av + V;)(-Av)'"(-£'Ah) 

+ (-Av)"'(-£'Ah) 
-< (-Av)'"(-£'Ah)(-£'Ah-Av + V,) 

+ [Ve, (-Av)"'(-£=^Ah)] + (-e^Ah - Av + K)™+' 
-< (-e^Ah - Av + 1^.)"'+' + [Ve, (-Av)"'(-£'Ah)] , 
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where in the last step we used the statement for M = 1 and again that lower 
powers of (— e^Ah — + V^) are operator-bounded by higher powers. To 
handle the remaining term on the right hand side we choose a partition of 
unity corresponding to the covering as in Remark[6]and orthonor- 

mal sections (T/)j=i^...^rf of TQj for all j. Then it holds 

i ' i i i ^3 i 

hi hi » 

The finite multiplicity of our coverings implies 
: ■ Jn.xR'' "^^ Jnc 



= (V'l-e'Ah^) 
< \\-6^A^^lj\\ + 

Therefore 



[-e'A^){-A^r + J2^jeVl,{-A^r + i-\ 



m 



^ {-e'Ai,)i-A,r + (-A 
^ (-e^Ah - A. + V;)"+'. 

We prove iii) only for M = 2 which is the hardest case. We notice that 
{Xiyr[H^,{Xu)-"^] = {\vr[H,,{\v)-^]Hl + {\vrH,[H,,{\v)-^]H, 

We also only treat the hardest of these summands which is the last one. 
The arguments below also work for the other summands and for M G {0, 1}. 
Inside of B25 the estimate iii) can be reduced to standard elliptic estimates 
as in i). Therefore we may replace if^ by —e'^A^ — A^ + because both 
operators coincide outside Bs. In view of ii) of Lemma [8], we have 

A-i (Az/)™ {-E^A^ - A. + [-e^Ah - A. + {Xum 
= (Az/)™ (-e^Ah - A. + V,)^ [-A., (Az/)-™] 

= ((Az/)™(-A. + K)M-Av,(Az/)-™] + (Az/)™[-A.,(Az/)-'"](-eX)2 



+ 2 (Az/)™ (-A. + Ve) [-A., (Az/)-'"] (-£^A 
+ (Az.)" [-£2 Ah, K] [-A., (Az.)-'"]) A-i 



62 



Because of = d^^^dnAp the operator (Az/)™(-A^ + K)' [-A^, (Az/)-'"]A~i 
contains only normal partial derivatives. It has coefficients bounded inde- 
pendently of A for any /, as the commutator [— Av, (Az/)"™] provides a A 
due to Lemma [8] iii). So by i) of Lemma [10] it is bounded by (— Av)^^^. 
Then ii) of Lemma [9] immediately allows to bound the ffist three terms by 
(— e^Ah — Ay + V^)^. The last term can be treated as follows. In the proof 
of ii) we saw that [— e^Ah, K] -< — e^Ah. Therefore 

h^'Ah, K] [-A., (Az/)-™]A-i 
= [-£2Ah,K](Az/)™[-A.,(Az/)-™]A^i 

-< -e^A^ {Xu)"' [-A., {Xuy^] + {Xv^ h^v, (Az/)-"] A'^ 
= {Xur [-Av, (Az/)-'"] A-i (-£2 Ah) + {Xur [-Av, (Az/)-'"] A-^ 
-< (-AO (-^^Ah) + (-£^Ah) + (-AO 

which is bounded independently of A by (— e^Ah — Av + K)^ again due to ii). 
Here again [— A^, (Az/)"™] has provided the lacking A. 

In view of (1691) and (1701) . the estimate iv) follows directly from i) of this lemma 
and iii) of Lemma [TUJ A polynomial weight is nescessary because here the 
unbounded geometry of {NC, g) really comes into play. In i) we could avoid 
this using that the operators differ only on a set of bounded geometry, while 
in ii) and iii) the number of horizontal derivatives was small! □ 

Proof of Lemma [TOt 

The first estimate is just an elliptic estimate on each fibre and thus a con- 
sequence of the usual elliptic estimates on R^. To see this we note that 
Av = S'^^dn^dnp is the Laplace operator on the fibers by (172|) and that the 
measure d^® dv = dn fi{x)dx is independent of n. 

To deduce the second estimate we aim to show that 

d'^^l'^ dn ixix)dx (76) 

< V / / \dl^\{\i-e^Ai,-A,)^\ + \eV''ij\ + \^\)dnfi{x)dx. 

with a constant independent of j. Then the claim follows from the Cauchy- 
Schwarz inequality and || |£:V'^-?/'| || = —e'^Ai.ip)^ < {ip\{—e'^Ai^ — A^)ip)2 
which is smaller than ||(— e^Ah — Av)\E'|| + ||^E'||. We note that here and in 
the sequel there is no problem to sum up over j because the covering 
has finite multiplicity! 
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On the one hand, there are a, (3 E {1, . . . , k} such that 




dnllljj* dn^dn^dnfli' du fl{x)dx 
dnfidnfi^J* dn'^dn^lp du fi{x)dx 

In, JW" 



dnpdn/sip* A^ip dn fi{x)dx. 





On the other hand, 



< 




/ / g^eV^dniiip* jSV^dniiip) dn n{x)dx 

g'^^eipxi +Tfi^n''dn'^)dnPip* e[d^i +T^^n^dnv)dni3t/j dn jj,{x)dx 
[ I + Tln^dn'^)d^,d^,re{d^i + r>^9„.)^ 

— e g''''T1^dncdn0'ip* £{d^i + T'^^n^dnv)ip dnfi{x)dx 
- 2elm(g'^T'^f^dnC'd^^tp*€{d^i +Tln^dnv)tp) dnfx{x)dx 



with Im(a) the imaginary part of a. When we add the last two calculations 
and sum up over all multi-indices 7 with I7I = 2, we obtain the desired 
(— e^Ah — Av)-term. However, we have to take care of the two error terms 
in the latter estimate: 

g''^V%dnai}j* T^gdnv^ dn fj,{x)dx 





—g^ Tlf^dnvdn^ijj* T'^f^ijj dn n{x)dx 
< supl/rf^r^'^l V / [ \d^,d^c.r\\^\dnfi{x)dx 

7|=2 J 
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and 




'X. 



This yields f l7^ because (7*' and can be bounded independently of j in 
our coordinates due to the bounded geometry and the smooth embedding 
of C assumed in ([5]) and ([6]) as explained in Remark [61 
To see that iii) is just a reformulation of iv), we replace n by iV = en in iii), 
put in V = I^;^, and use that {-e^Ai,-A^+V,)D* = D;{-e'^Ah-e'^A^+V^) 
by Lemma m 

So we immediately turn to iv). We notice that the powers of e on both sides 
match because all derivatives carry an e. Therefore we may drop all the 
e's in our calculations to deduce the last estimate. Since we have stated the 
estimate with a non-optimal power of (u) , there is also no need to distinguish 
between normal and tangential derivatives anymore. So the multi-index a 
will be supposed to allow for both normal and tangential derivatives. We 
recall that As = + A^. We will prove by induction that for all m G Nq 



with a constant independent of j. Applying this estimate iteratively we 
obtain our claim because as explained before —As + is self-adjoint and 
thus (—As + Vy is operator-bounded by (—As + V)"^ for I < m due to the 
spectral calculus. 

Before we begin with the induction we notice that, in view of Proposition [T], 
g'^ is positive definite with a constant that is bounded from below by 
times a constant depending only on the geometry of C. More precisely, the 
constant depends on sup and the inverse constant of positive definiteness 
of which are both uniformly bounded in our coordinates again due to ^ 
and 

We start the induction with the case m = 0. For |a| = there is nothing to 
prove. Since = detg^^ by Proposition [1], it holds As = fJ'^^da fJ^ g^^ db. So 




(77) 
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for lal = 1 we have 



J2 [ [ {i^y^ld^^l^dNfidx < [ [ gtdai^*db^dN fidx 

„,i_i -^K* J^i "'K* 

ip* da figs'' dbi' dN^ dx 
[ [ ip* {{-As + V -V)ip)dN fx dx 



< ||^||(||(-As + \/)V^||+sup|V| 

< u-As+vw + 

< (\\i-As + V)^\\ + im?. (7J 



Taking the square root yields the desired estimate in this case. For |q;| = 2 
we have 



\a\=2 



< 



y2 f f {^y^'^ld'i^^dNfidx 

I {y)-^^gtdad,rdbdctpdNfidx 

V / / -{i^ygtdadcd.rdb^ 

- /i"^ {dcij{u)-^^gf)dadcij* db-^ dNjj dx 

-(fi'\d,fi{u)-'^gt)dadcr - {da{i^)-'^) gtdcdcr) db^dNfxdx 



< 



|o|=2 J 



which yields ( 1771) via (178|) when we apply the Cauchy-Schwarz inquality and 
devide both sides by the square root of the left-hand side. Here we used that 
both fi~^{dcfi{i^)~^^gf) and {da{y)~^'^) g^ are bounded by {y)^^"^ ■ This is 
due to the facts that the derivatives of /i are globally bounded due to the 
bounded geometry of C, that g'^ and its derivatives are bounded by (z/)^ 

due to Proposition [T|, and that any derivative of (z/)' = a/1 + dapW^nl^ is 
bounded by iyY . We will use these facts also in the following calculation. 

We assume now that fl77|l is true for some fixed m G Nq. Then it suffices to 
consider multi-indices a with lal = m -|- 3 to show the statement for m -|- 1. 
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We have 



\a\=m+3 



f3 ^ 



|a|=m+3 |/3|=m+l 



a|=m+3 |a|=m+2 ' 

where we used that [Ag, 9"] includes no terms with more than m + 3 partial 
derivatives and that its coefficients are bounded by (z/)^ times a constant 
independent of e. Using again the Cauchy-Schwarz inequality, deviding by 
the square root of the left hand side, and applying the induction assumption 
to the a-term we are almost done with the proof of (I77|l for m + 1. We 
only have to introduce V in the Laplace term. We recall that it follows from 
Vc,W e C^{C,C^{NgC)) that K e C^{C,C^{NgC)). When we put it in 
and use the triangle inquality we are left with the following error term: 

E / / wy^'^^^id^y^i^dNi^dx 

E / / {u)-^\''\\d''Vmu)-''\^\\d^^\UNi^dx. 

|a| + |p|=m+l -I 

In order to apply the induction assumption to this expression, we have to 
bound sup(z/)-8|"l|5"i/|2. To be able to use V e C^{C,C^{NgC)) we ffist 
replace the tangential derivatives in 9" by V'^ and afterwards the normal 
derivatives by V^. In view of (l69l) and (1701) . this costs at most a factor {i')^^ 
for each derivative. □ 
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We still have to give the proof of Lemma [3] from Section 1X21 It was postponed 
because it makes use of Lemma [TDl 

Proof of Lemma [3] (Section 13.21) : 

All statements in i) and ii) are easily verified by using the substitution rule. 
To show iii) we first verify that {H^A* - A*H^)P^ is in C(V{H^), L'^iA, dr)) 
at all. For A*H^P^ this immediately follows from ii) and Corollary O So we 
have to show that H%A*P'' E C{V{H^), L^{A, dr)). By Corollary [5] we have 

\\H^j^A*P''\\c(v{H^),L'2(A,dT)) ~ WH^j^A* {iy/E)'''\\c('D(H^),L^A,dr)) 

for any / G N. Now we again fix one of the geodesic balls Qj C C of a 
covering as in Remark [6] and choose geodesic coordinates (a;*)j=i,..._rf and 
bundle coordinates (/^")o=l,...,A: with respect to an orthonormal trivializing 
frame (t'^)a over Qj. When we write down A* and in these coordinates, 
we will end up with coefficients that grow polynomially due to our choice 
of the diffeomorphism $ and the metric g. However, this is compensated 
by {h'/e)~K Choosing / big enough allows us to apply Lemma [10] iii) to 
bound H%A*{u/e)-^ by -e^Ah - e'^A^ + V . The proof of Lemma [9] i) 
also shows that — e^Ah — e^A^ + V'^ H^. To sum up over j is once 
more no problem because the covering has finite multiplicity. Hence, 

H^A*P'^ G C{T>{H^),L?'{A,dT)). With the same arguments one also sees 
that \\A*{u/efA {H'^A* - A*H')P'\\^T,^H^),LHA4r)) < 1- 
Since g is by definition the pullback of G on Bs/2, the operators H^A* and 
A*H'^ coincide on functions whose support is contained in Bs/2- But outside 
of Bs/2, i-e. for > 6/2, we have that 

Hence, denotig by Xbs/2 characteristic function of NC \ Bs/2 we obtain 
that Wxsg 2('^/£^)^||oo ^ Using that A*iIj = on ^ \ ;B for all ijj and 
AA* = 1 by ii) we may estimate 

||(/7^y4* - A*H'')P''\\c{V{H^),L^{A,dr)) 

= WA^xl^^A {H:,A* - A*H^)P^\\civiH^)^LHA4r)) 

= \\A*XB,,M/er'AA*{u/efA {H^^A* - A*HnP^\\c(^(H^)^LHA4r)) 

^ \\A*XBs,,i^/^y'MciLHA,dr)) 

= \\Xbs/M/^)'^\\-- 

< 6^ 

which was the claim. □ 
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4.2 Expansion of the Hamiltonian 

In order to expand the Hamiltonian if^ in powers of e it is crucial to expand 
the metric 'g around C because the Laplace-Beltrami operator depends on it. 
The use of the expansion will be justified by the fast decay of functions from 
the relevant subspaces Pq and in the fibers. 

Proposition 2 Let 'g be the metric on NC defined in l[W\} . Choose Q G C 
where the normal bundle NC is trivializable and an orthonormal frame (z/a)a 
of NC\^ as in Remark\^ In the corresponding bundle coordinates the inverse 
metric tensor g G Tq{NC) has the following expansion for all q G Q: 

^ = {c^ ?) (o I) (I T) + 

where for i,j, l,m = 1, d and a, /5, 7, 5 = 1, .., k 

A'^{q,n) = g'Ko) + n-{^:i9'' + f^^,){q) 

S^^(g,n) = 6,s + \n^nPR\%{q), 
C]{q,n) = -n-rjjg) + | n°n^ R^,,^(g). 

Here R denotes the curvature tensor of A andWa is the Weingarten mapping 
corresponding to Ua, i-e. W(z/q) (see the appendix for definitions). The 
remainder term vi and all its derivatives are bounded by \n\^ times a constant. 

For the proof we refer to the recent work of Wittich |35]. He does not 
calculate the second correction to C but it is easily deducable from his proof. 
Furthermore, Wittich actually calculates the expansion of the pullback of G, 
which coincides with 'g only on Bs/2- Then vi is only locally bounded by \n\^. 
To see that the global bound is true for g we recall that outside of Bs it 
coincides with g^, which was explicitly given in Proposition [TJ Comparing 
the expressions for 'g and gs we obtain a bound by |np which is bounded by 
\n\^ times a constant for \n\ > 6. 

In addition, we need to know the expansion of the extra potential occuring 
in Lemma [H which is also provided in |45j : 

Proposition 3 For p := ctji/da with da = dfi ® dv it holds 

Vp{q,n) = -jg(q^o)iv,v) + - + trc Ric + trc R) (g) + r2{q,n) 

=■■ V"geom(g) + r2{q,n), 

where rj is the mean curvature normal, k,k are the scalar curvatures of C and 
A, tic Ric, trc R are the partial traces with respect to TqC C TqA of the Ricci 
and the Riemann tensor of A and r2 is bounded by \n\ times a constant. 
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Again there is only a local bound on r2 in [45J. In our setting the global 
bound follows immediately from the coincidence of dJI and da outside of Bs-, 
see (|T2|) . With these two inputs the proof of Lemma[6]is not difficult anymore. 

Proof of Lemma [6] (Section 13.31) : 

Let P with \\{j-'Y P\\c{v{H^+'^)v{H.^)) ~ ^'^^ / G No be given. The similar 
proof for a P with \\P {i'Y\\c{v{h^),h) ^ 1 fo^^ ^ ^ will be omitted. 
We choose a covering of C of finite multiplicity and local coordinates as at 
the beginning of Section [Hand start by proving ^ 1 for 

j e {0,1,2}. Exploiting that all the coefficients in Hj are bounded and 
have bounded derivatives due to the bounded geometry of A and C and the 
bounded derivatives of the embedding of C assumed in ^ and (jHD we have 

\\HjP\\c{V{HT+^),H) ~ \\Hj{yy^^\\c{V{H,),H) 

|a| + |/3|<2 

^ \\He\\c(V(H^^),H) = 1, (79) 

where we made use of Lemma [10] iii) and Lemma [9] for the bound by if^. 
Now we set ipp := Ptp. By definition of if^ and it holds 

(0|ff,7/;p) = {<f)\D;M;{-e^Ag + V')M,D,^p) 

= {(/) I D;M;i-E'Ag)MpD,^Pp) + {<P I (K + DlWD,)^ljp). (80) 

Due to ||(z/)'^P|| < 1 a Taylor expansion of D*WD^ in the fiber yields 
DlWD,{q, u)P = {Wiq, 0) + e{V:W)iq, 0) + le'iVl,W)iq, 0))P + 0(8^). 
Recalling that Vo{q, v) = K(q', i^) + W{q, 0) we find that 

(0 I (K + DlWD,)i,p) 

= {(l)\{V^ + e{V''W){q,Q) + \e\VlW){q,f)))i,p) + 0{e^). (81) 

The error estimate in Proposition [3] yields that \\Dlr2D fr{v)^^ilj\\ < and 
thus \\Dlr2Di;iljp\\ < So Lemma [Hand Proposition [3] imply that 

{^\DlM;{-e^Ag)MpD,^p) 

= [ [ e'g{dD,(l)*,dD,^p)diydfi + e^<j)\D;V,D,^p) 

JC JNqC 

= [ [ e^g{dD,(l)*,dD,^p)diydfi + E^(l)\Vg,,^iljp) + Ois'), (82) 

Jc JNgC 

where we used that V^eom does not depend on v. 
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Next we fix one of tlie geodesic balls C C of our covering and insert the 
expansion for 'g from Proposition [2] into ([H2D- Noting that d^iD^ = D^d^i and 
dncD^ = e'^D^dnc we then obtain that 

Jn J NqC 

+ B''^{q,n)d^0D,i/jp dndfi + 0{e^) 

+ B''^{q,en) d^0ip + (j)*V,{q,n)'^p dndfi + 0{e^) (83) 

because the bound on ri from Proposition [2] allows to conclude that the term 
containing DlriD^ is of order e^. To do so one bounds the partial derivatives 
by as in (1791) . After gathering the terms from (IHOj) to (l83i) and plugging 
in the expressions for A, B, and C from Proposition [2] the rest of the proof is 
just a matter of identfying and via and (17U]1 . When we sum up 
over the whole covering, the error stays of order because our covering has 
finite multiplicity and the bounds are uniform as explained in Remark [61 □ 

4.3 Construction of the superadiabatic subspace 

Let Ef be a constraint energy band. We search for G C{T-C) with 

i) P,P, = P„ 

ii) [H,,P,]xm = 0{e') 

The former simply means that is an orthogonal projection, while the latter 
says that P^x{H^)T-C is invariant under the Hamiltonian up to errors of 
order e^. 

Since the projector Pq associated with Ef is a spectral projection of Hf, we 
know that [H{, Pq] = 0, [Ef, Pq] = 0, and HfPo = E^Pq. Lemma [6] yields that 
= Hq + 0{e) with Hq = — e^Ah + Hf. So Pq satisfies, at least formally, 
[H„ Po] x{He) = [-e'Ah, Po] xiH,) + Oie) = 0(e). Therefore we expect P, 
to have an expansion in e starting with Pq: 

P, = Po + ePf + e^P2 + O(e'). 

We first construct P^ in a formal way ignoring problems of boundedness. 
Afterwards we will show how to obtain a well-defined projector and the asso- 
ciated unitary Us- We make the ansatz Pi := T^Pq + PqTi with Ti : TC ^ H 
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to be determined. Assuming that [Pi, — e^Ah + E{] = 0{6) we have 

[H,,P,]/e = [Ho/e + HuPo + eP,] + 0{e) 

= [Ho/e + HuPo] + [Ho, Pi] + 0{e) 

= [-£Ah + H,, Po] + [^f - Ef, Pi] + Ois) 

= [-^Ah + H,, Po] + (H, - E,)T;Po - PoTi{H, - E,) + 0(6) 

We have to choose Ti such that the first term vanishes. Observing that 
every term on the right hand side is off-diagonal with respect to Pq, we may 
multiply with Pq from the right and 1 — Pq from the left and vice versa to 
determine Pi. This leads to 

- {Hi-Eiy'{l-Po){[-eA^,Po] + Hi)Po = {l-Po)T*Po (84) 

and 

- Po ([Po, -eAh] + Hi) (1 - Po) {Hi - E,)'^ = PoT, (1 - Po), (85) 

where we have used that the operator Hf — E{ is invertible on (1 — Po)7Yf . In 
view of and (jHSD, we define Ti by 

Ti := -Po{[Po,-eA^]+H,) RH,{Ei) + RH,{Ei) {[-eA^,Po]+H,)Po (86) 

with Rhi{E{) = (1 — Po)[H{ — El) ^(1 — Po). Ti is anti-symmetric so that 
p(i) ._ _|_ _ _|_ ^(J'*Pq _|_ PoTi) automatically satisfies condition i) 
for Pe up to first order: Due to Pq = Pq 

p(i)p(i) = p^^^{t-p^^p^Ti + Po{TI + Ti)Pq) +0{e') 
= Po + e{T*Po + PoTi) + 0{e^) 
= P« +0{e^). 

In order to derive the form of the second order correction, we make the ansatz 
P2 = T* PqTi +T2*Po + PoT2 with some T2 : Ti ^ Ti. The anti-symmetric part 
of T2 is determined analogously with Ti just by calculating the commutator 
[Pg, H;.] up to second order and inverting Hf — Ef. One ends up with 

{T2-T;)/2 = -Po{[P^'\H'^'^]/e')RHAEf) + RHAEf){[H^'\P^'^/e')Po 

with H^"^^ := i^o + ^Hi + e^H2- The symmetric part is again determined by 
the first condition for P^. Setting P*^^) := P^^^ + £:^P2 we have 

p(2)p(2) ^ p(2) ^ e^i^p^TiTlPo + Po(T2* + r2)Po) + 0{e^), 

which forces T2* + T2 = —TiT^ in order to satisfy condition i) upto second 
order. 
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We note that Ti includes a differential operator of second order (and T2 
even of fourth order) and will therefore not be bounded on the full Hilbert 
space and thus neither P^. This is related to the well-known fact that for a 
quadratic dispersion relation adiabatic decoupling breaks down for momenta 
tending to infinity. The problem can be circumvented by cutting off high 
energies in the right place, which was carried out by Sordoni for the Born- 
Oppenheimer setting in [12] and by Tenuta and Teufel for a model of non- 
relativistic QED in ^3] . 

To do so we fix ii^ < 00. Since is bounded from below, := mia{He) 
is finite. We choose Xe+i £ C*^(M, [0, 1]) with Xe+i\{E--i,e+i] = 1 and 
supp Xe+1 C. {E_ — 2, E + 2]. Then we define 

,= p(2)_Pq = e{T*Po + PoT^)+e\T*PoT, + T;Po + PoT2) (87) 

and 

PXE+I ,= P, + P^j.j,^,(H,) + XE+imPe{l-XE+im) (88) 

with xe+i{Hs) defined via the spectral theorem. We remark that P^'^+^ is 
symmetric. 

We will show that p^'^+^ — Pg = 0{e) in the sense of bounded operators. 
Then for e small enough a projector is obtained via the formula 

where r = {z & C||2 — 1| = 1/2} is the positively oriented circle around 1 
(see e.g. [H|). Following here the construction of Nenciu and Sordoni [31] 
we define the unitary mapping : P^Ti — > P^Ti by the so-called Sz-Nagy 
formula: 

ij, := (PoPe + (1 - Po)(l - Pe)) (1 - {Pe " ^o)')"'^'. (90) 

We now verify that P^ and f/g have indeed all the properties which we stated 
in Lemmas [2] & [7| and state here again for convenience: 

Proposition 4 Fix E < oo. Let E{ he a simple constraint energy hand and 
Xe+1 e C°°(M, [0, 1]) with xe+i\{-oo,e+i] = 1 and suppxs+i C (-oo, E + 2]. 
For alle small enough P^ defined hy ^^Sl^j-^EW ^■^ a hounded operator on Ti and 
Us defined hy ^9^) is unitary from P^Ti to PqH. In particular, P^ = f/*Pof/e. 
For all m G No and Borel function x '■ ^ [—1) 1] with supp x C {—oo, E+1] 
it holds \\Ps\\c{v(Hp)) ^ 1 and 

Wl^e, Pe]\\c{V{H]r+^),V{Hp)) = ^(^)' \\[He, Pe] XiHe)\\c(nMHr)) = ^{^^)- 
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Furthermore, it holds Ef G C^{C), as well as: 

z; Vj,/G No, mG {0,1}: \M Psii^y WaviH^)) < 1. 
tt) Vj,/GNo: \MPo{i^y\\ciViH,)) < 1, \\he'A^,Po]\\civiH.),n) < e. 

in) There are t/f,f/| G C{l-C) fl C{V{Hi;)) with norms bounded indepen- 
dently of e satisfying PqUIPq = and t/|Po = P0U2P0 = P0U2 such 
that U, = l + ellf + e^f/f. In particular, WU^ - iWan) = 0{e). 

iv) \\PqUI{uY\\c{v(h^)) < 1 for all I eNq and me {0, 1}. 

v) For B, := PqU,x{H,) and all u G {1, (f/f)*, (f/|)*} it holds 

vi) For Rh,{Ei) := (1 - Po)(^f - Ei)~\l - Pq) it holds 

vii) Ifipi G C^{C,nf), It holds 

\\Uo\\c{V{H,),V{-e^Ac+Et)) ^ 1> II f^O II /:(I>(-e2Ac+£;t),»(^^6)) ~ 1' 

and there is Xq > I with sup^ ||e'^°^'^Vf(Q')ll?i;f(ij) ~ 1 ^'^'^ 

sup||e^«<'^V;,...,,,V:^^,...,,^^f(g)||^,(,)<l 
<? 

/or a// i/i, . . . , z/; G rb(A^C) and n, . . . , G rb(TC). 

The proof relies substantially on the following decay properties of Pq and the 
associated family of eigenfunctions. 

Lemma 11 Let Vq G C^{C,C^{NC)) and Ef be a constraint energy band 
with family of projections Pq as defined in Definition 
Define V^^^Pq := [Vi^^,Po] (ind, inductively, 

'^Ti,...,r„,Po ■= [V^^, V^2,...,r™^0] - Y.%2^T2,...,\7r,Tj,...,TmP0 

for arbitrary ri, . . . , G r(TC). For arbitrary z/i, . . . , z/; G r(A^C) define 
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i) Then E{ e C^{C), Pq e C^{C, CO<f)) , and there %s Xq > independent 
of e such that for all X G [— Aq, Aq] 

and 

11^ ^ri,...,r„-^o;e \\c{n) ^ 

for all ui, ... ,1^1 G rb(A^C) and Ti, . . . , G Ti,(TC). 

Let E{ be simple and (ff be a family of normalized eigenf unctions that define 
a smooth section of the associated eigenspace bundle . 

ii) Ififi G Cl^{C,Hi), then G Cl^{C,C^{NC)). Furthermore, 

sup||e^«<-Vf(g)lk,(.)^l, sup ||e^«<'^V;,.,„,^V:^^,...,,„^f(g)||^,(,) < 1 
(JSC qec 

for all ui, ... ,1^1 G rb(A^C) and ri, . . . , G Ti,(TC). 

Hi) If C is compact or contractable or if Ef{q) = inf cr(iJf(g)) for all q E C, 
then ipi can be chosen such that Lpi G C^(C,7if). 

In addition, we need that the apphcation of XE+i{,He) does not completely 
spoil the exponential decay. This is stated in the following lemma. We notice 
that we cannot expect it to preserve exponential decay in general, for we do 
not assume the cutoff energy E to lie below the continuous spectrum of Hg,] 

Lemma 12 Let x ^ C'^(]R) be non-negative and (^H,V{H)^ be self-adjoint 
on Ti. Assume that there are I ^ Z,m E N and Ci < oo such that 

\\{Xuy[W,{Xu)-']\\civiH^),n) < CiX (91) 

for all X G (0, 1] and 1 < j < m. Then there is C2 < oo independent of H 
such that 

\\{iyyxiH){u)-'\\cin,viH-)) < C[C2. 

This lemma can be applied to for m < 3 in view of Lemma O Now we give 
the proof of the proposition. Afterwards we take care of the two technical 
lemmas. 

Proof of Proposition [U 

We recall that V{H^) := H and E_ := ini a{H,). Let xe e C^{R, [0,1]) 
with Xe\[e.,e] = 1 and suppxe C [E_ — 1,E -\- 1]. Then by the spectral 
theorem XE{,H^)x{He) = xiHe) and XE+i{,He)xE{,He) = XsiHe) for x and 
Xe+1 as in the proposition. In the sequel, we drop all e-subscripts except 
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those of and write x,Xe, and xe+i for x{He),XE{H^), and xe+i{H^) 
respectively. 

The proof of the proposition will be devided into several steps. We will often 
need that an operator A G jC{H) is in C{V{Hl),'D{H^)) for some l,m e No- 
The strategy to show that will always be to show that there are li,l2 G N 
with h + I2 < 21 such that for all j G Nq 

(-^'Ah - A, + K)"^ -< H-^'(V")'n^V'^)''- (92) 
Then we can use Lemma [9] to estimate: 

WH^^A^I^W + \\Aij\\ < \\(-e^A^-A^ + V,rAij\\ + M 

< ||(z/)-^'-^'Hv^)'H^v^)'¥ll + ll^ll 

< WHi^W + IIV^II, (93) 
which yields the desired bound. 

Step i.- 3 Ao > 1 V A < Ao, m G No : He^^'^^ Pq e^^'^^ WcmHs^)) < 1 and 

Both statements hold true with e^^'^'^ replaced by (z/)' for any I eNq. 

Let Ao be as given by Lemma fTTl When we choose a partition of unity (^j)j 
corresponding to the covering as in Remark E] at the beginning of 

Section m and orthonormal sections {i^l,)a=i,...,k of NQj and (T/)j=i of TQj 
for all j, the coordinate formulas fl72|) imply 

In order to obtain the estimate flU^ for A = e'^"^'^^ Pq we first com- 
mute all horizontal derivatives to the right and then the vertical ones. Using 
Vq G C^[C,C^{NC)) and Lemma [8] we end up with terms of the form 
e^^"^ {VI, V^,. Po)e^<'^> (V^)'i {eV^Y^ times a bounded function with 

^1 + ^2 < 2m. By Lemma [TT] we have 

V^, Po)e^<''>(V")'H£V'^)'' ^ e-(^''-^)<''>(V")'H^V'^)'' 

which implies (192|) due to A < Aq. This yields the first claim of Step [1] 
via ([93]). 
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The second claim can easily be proven in the same way. For the last claim 
it suffices to notice that \\{i'Yq~^°^'^'^\\c{V(h^^)) ^ 1 for all l,m E Nq, which is 
easy to verify. 

Step 2: It holds V A < Aq, m G No, i G {1, 2} : 

In particular, V A < Acm G No : We'^^"'^ P ^'^^"'^ W c{v{h^+^) v(h^)) ~ ^• 

The last statement is an immediate consequence because by definition of P 

qX{u)P^x(u) ^ ^^x(u) (^i^T*Po + PoTi) + e{T*PoPoTi + T;Pq + PoT2))e^<^\ 

We carry out the proof of the first estimate only for T*Pq. The same argu- 
ments work for the other terms. To obtain fl92|) for A = e'*'^''^T*Po e'*'^''^ we 
again commute all derivatives in (— e^Ah — Av + K)™" and T*Pq to the right. 
In view of fl86l) . the definition of Ti, we have to compute the commutator of 
RHiiEf) with and V^. For arbitrary r G T^{TC) it holds 

[V^i?Hf(i?f)] = -{V'^Po)RH,{Et) - RH,{Ef){V'^Po) 

-RH,{Ef)[V'^,Hf-Ef]RH,iEf). 

with [V^,Hf - Ef] = {V^Vo - VrEf). The latter is bounded because of 
Vo G C^ {C,C^{NgC)) by assumption and E^ G C^{C) by Lemma [IH An 
analogous statement is true for V^. Hence, we end up with all remaining 
derivatives on the right-hand side after a finite iteration. These are at most 
2m + 2. After exploiting that We^^"^ RHtiEi) e-^^''^\\c(n) < 1 by Lemma [H] 
we may obtain a bound by H^~^^ as in Step [H 

Step <?.■ Vm G No : \\P^''+'\\c{V{hv-)) < 1 and 

Vj, /, G No, mG {0,1}: || H^P^-+^ < 1- 

We recall that P^^+i -^^g defined as 

Step d] implies that Po ^ C(V{H^)) for all m G No. So it suffices to bound 
the second and the third term to show that P^^^+i g C{V{H^)). Since 
is bounded from below and the support of Xe+i is bounded from above, 
\\Xe+i\\c{'H,v{h^)) ^ 1 for every m G No. So the estimate for P obtained in 
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Step [2] implies the boundedness of the second term. By comparing them on 
the dense subset V{H^) we see that xe+iP is the adjoint of Pxe+i and thus 
also bounded. This finally implies the boundedness of the third term, which 
establishes \\P^'^"^^\\c{v{h^)) ^ 1 for all m G Nq. 
We now address the second claim. We fix A with < A < Aq. Then 

+ {uye-^^'^"^ (e^<'^>PeA(.>) (e-M-)(^)') {u)-\e+M 
+ WXE+i{y)-'{{^Ye-^^''^) (e^<'^>PeA(.)) 
x(e-^<^)(z.)') {v)-\l-XE^i){vy 

It is straight forward to see that \\{i'yQ~^°^''''\\c(v(H]j^)) ^ 1 for all j, m G Nq. 
Therefore Step [T] yields the desired estimate for the first term. In addition, 
we know from Lemma [T^ that \\{i')^^XE+i{i-'y\\c{n,v{H'-^)) ^ 1 because 
satisfies the assumption of Lemma fT2] due to Lemma [9] iii) . So Step [2] implies 
the desired estimate for the second term. Then it also follows for the third 
term again by estimating it by the adjoint of the second one. 

Step 4: It holds V m G No, i G {1, 2} 

II [TiPo, -e^Ah + Ef] \\c(v(Hr+'+''),v{Hp)) = ^'(e), 
II [PoTi, — e^Ah + Ef] ||£(X)(//™+'+i),o(_f/|")) ~ ^(^)- 

We again restrict to T^Pq because the other cases can be treated in quite a 
similar way. 

We note that E{ commutes with all operators contained in T*Pq but eV^. 
Furthermore, \\[eV^, Ei]Po\\c(viHs-)) = £\\C^tE{)Po\\c{v{hp)) = O^e) for any 
r G rb(rC) by Lemma [TTl With this || [T^i Pq, -£'f]||£(x)(_H-^"+2),x>(_ffm)) = C^(^) is 
easily verified. 

We will obtain the claim of Step H] for T^Pq, if we are able to deduce that 
||[Ti*Po, -^^Ah]||£(p(j|^™+2) .p(j:/m)) = 0{e). Again we aim at proving by 
commuting all derivatives to the right. In Step [1] and Step [2] we have already 
treated the commutators of — e^Ah with Pq an RHf{E{). So it remains to 
discuss the commutator of eV^ and — Ah, which does not vanish in general! 
To do so we again fix a covering {Qj)j(zfq of C and choose a partition of 
unity corresponding to the covering as in Remark El as well as 

orthonormal sections {T-)i=i^,,,^d of TQj for all j. 
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Recalling from ^ that Ah = ^ti 0(V^ - ,) we have 

i 

j,i * 

+ V;:,R'^(r,r/) + V;:,Vj; + [V^V^ ,^,)]). 

i 

In view of the expression for R*^ in Lemma all these terms contain only two 
derivatives. So we have gained an e because, although R^ and its derivatives 
grow linearly, we are able to bound the big bracket as required in fl92|l using 
the decay provided by Pq. The estimate is independent of Qj because R^ is 
globally bounded due to our assumption on the embedding of C in (jH]). 

Step 5: For all m G No 

\\[Hs,P''^^mciV(Hr^^)MH-)) = ^(^)' me,P''^^']XE\\c(H,ViH-)) = 0{e'). 

We fix m G Nq. Due to the exponential decay obtained in Steps [T] & [2] for Pq 
and P we may plug in the expansion of from Lemma when deriving the 
stated estimates. The proof of Step [2] entails that P^^+i _ jg Qf order e 
in C{T>{H^)) for any m G Nq. Therefore 

\\[He-, P^'^'^^]\\c(V(HT'^^),V{H'^)) = \\[Hs, Po]\\c{V(HT'^^),V{H^)) + ^(^) 

= \\[Hq-,Po]\\c(v{ht"''-),v(h^)) + ^(^) 

by SteplU On the other hand we use [if^, xe] = and (1 — Xe+i)xe = to 
obtain 

\\[He,P^''+^]XE\\c{H,V{H^)) 

= II [He, P^^^] XElldHMHS')) 

= \\[He, Pq + P]XE\\c{n,v(H^)) 

= \\[Ho + eH^ + e''H2,Po + P]xE\\cin,viHp)) + Oi6') = ^(e'), 

where the last estimate follows from the construction of Ti and T2 at the 
beginning of this subsection (which were used to define P). To make precise 
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the formal discussion presented there one uses Step H] and once more the 
decay properties of Pq and P to bound the error terms by for some 
m G N as in ([92D and ([93]). 

Step 6: For e small enough P & U are well-defined, P^ = P, and UIp-h is 
unitary. \\P\\c{v{Hp)) < 1 and \\P - Po\\c{v{h^^)) = 0{e) for all m G Nq. 

Since Pq is a projector and ||P>^£;+i — Po\\c{h) = C^(^) by the proof of Step [3], 
we have 

||(px.+i)2_pxK+i||^^^j = 0(^). (95) 

Now the spectral mapping theorem for bounded operators implies that there 
is a C < cx) such that 

a{pxE+i) C [-Ce, Ce] U [1-Ce,l + Ce]. 

Thus P '■= §Y (-Pe"^^^ — z) ^ dz is an operator on Ti bounded independent 
of e ioT e < 1/2C and satisfies = P by the spectral calculus (see e.g. [S]). 
By the spectral theorem P = X[i-Ce,i+Ce]{.P^'^^'^) and so ||P — P^'^+^WciH) = 
0{e). With ||pxB+i _ Po||^(„) = d{e) this entails ||P - Po\\c{n) = 0{e). 
Hence, 1 — (P — Pq)^ is strictly positive and thus has a bounded inverse. 

Therefore U := (PqP + (1 - Po)(l - P)) {l - {P - Po)^)~^^^ is also bounded 
independent of e as an operator on H and satisfies 

U = Uo{P + 0{e^)). 

We set 5:= [l - {P - Pofy^^\ It is easy to verify that [P, 1 - (P - Po)^] = 
= [Pq, 1 - (P - Po)2] and thus [P, S] = = [Po,S]. The latter implies 
U*U = 1 = UU*. So U maps PH unitarily to PqH.. Since Uq is unitary when 
restricted to PqH, we see that U = UqU is unitary when restricted to PH. 
The combination of with Steps [3] and \5\ immediately yields 

\\{P''-^^Y-P''-^^\\civm) = 0{e). 

for all m G Nq. So for e < 1/2C and z G 9Pi/2(l) the resolvent (P^-^+i 
is an operator bounded independent of e even on T>{H^). In view of P's 
definition, this implies ||-P||z;(d(//'")) ^ 1 for all m G N . Then we obtain that 
ll-P ~ Po\\c{v{H^)) = C^(^) in the same way we did for m = 0. 

Step 7: \\[He, P]\\c(p(^HT+^),v{H^)) = ^'(e) & \\[Hs, P]xe\\c{h,v{hp)) = ^(e^) 
/or all m G Nq. 

We observe that 

[He, P] = — I [P^^+^ - z)~^[Hs, P^^+^] [P^^+' - z)'^ dz. 
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Since we saw that || (P'^^+i — z) '^\\c{v{Hr-)) ~ 1 the preceding step, the 
first estimate we claimed follows by inserting the result from Step O To 
deduce the second one we set RpXE+i (z) := (P^-^+i _ -.^ g^^^j _ -^^-^ 

to compute 

[He,P]x = ^ f^Rp-E+i{z)[H,,P^^+^]Rp.E+^{z)xEXdz 



^- j RpxE+i (z) [H„ P^^+^]xi? Rp-E+i (z) X 

+ RpxE+i (z) [He, P^^+^] [RpxE+i (z), xe] X dz. (96) 



Furthermore, 



[RpXE+i {z),xe]x = Rp^E+i (z) [px^5+i , Xe] Rp^e+i (z) xe X 
= RpXE+i {z) [P^B+^ , Xe] Xe Rp^e+i {z) x 

+ RpxE+i{z) [P^^+SXi?] [Rp^e+i{z),xe] X 
= RpXE+i {z) [P^E+^ , Xe] Xe Rp^e+i {z) x 

+ {Rpxe+,{z) [P^^+\xe])^ Rp-E+.{z)x. 

Since due to Step [5] we have || [-P^'^^S -^^e] ||£(x)(//'"+i) ©(/^m)) = and 
\\[P'''=+\H,]xE\\c{n,v{Hr)) = ^(^^)' LemmaSyields 

lll^'^'^^'Xi^llLd^CHDW^^)) = ^(^)' \\[P''''^'^XE]XE\\m,viHr)) = 0{e''). 

Applying these estimates, \\RpxE+i{z)\\c{v{H^)) ^ 1, and Step [5] to (p6l) we 
obtain \\[H,,P] x{He)\\c{HMHr)) = ^(^')- 

Step 8: yjj eN, me {0,1}: 11^)' P {iyy\\c(v(Hi^)) < 1- 

This can be seen by applying the spectral calculus to P'^^+i which we know 
to be bounded and symmetric. Let / : C — C be defined by := z and 
let : C — s> {0, 1} be the characteristic function of P2/3(l)- Then due to (p5|) 
the spectral calculus implies that for e small enough 

p = g^pXE+l) = f(^pXE^.)(^glf)^pXE+.)f^pXE+r) 

= P^^+' {g/f){P^^+')P^^+K (97) 

We note that {g I P){P^'^^^) € C-iTi) because (7 = in a neighborhood of 
zero. Since g/ p is holomorphic on Pi/2(1), it holds 

{g/f){P^-^^) = {g/f){z)RpxE^.{z)dz 

81 



by the Cauchy integral formula for bounded operators (see e.g. [H]). In the 
proof of Step [H] we saw that \\RpxE+i{z)\\xi{H^) < 1 for 2; G (9i?i/2(l), which 
implies that also \\{g/ f^){P^'^^'^)\\c{v{H^)) ^ 1- Then applying the result of 
Step [3] to (1971) yields the claim. 

Step 9: \f me No: \\iP - P''''^')x\\cinMHr)) = ^^^'^ 

By construction we have Ti = —T* and T2 + T2 = — TiTj* as well as PqTiPq = 

0. With this it is straight forward to verify that P^^-* = Pq + P satisfies 

||X.(P(V(^)-P(^))X||,(^,,(,., = Oie% (98) 

Since \\[P^''+\H,] 

x\\c(n,v(H-^-'^)) = ^(^^) by Step[5l Lemma H] yields 
\\[P''''^\XE]x\\cin,viH-)) = Ois'). 

Recalling that \\P^'^+'^\\c(v(Hp)) ^ 1 due to Step [3] we have that in the norm 
oiC{n,V{H^)) 

= (P>^^+i - l)P^^+ixi?X 

= (pXB+i _ l)^^j^pXE+i ^ + (pXB+1 _ l)[pXE+i^ ^ 

= Xe{P''^+' - l)P™x + [P^^+^ X + 0{e^) 

= XE{P^'^-l)P^'^X + 0{e') 

= xe (P(2)p(2) - p(2)) ^ + 0{e^) # 0(£=^). 

Since we know from the proof of Step E] that \\RpxE+i{z)\\c{v{H^)) < 1 for z 
away from and 1, the formula 

p_pXE.r = J_ / i^px... jz) + Rp.E,. (1 - Z) _ 

27ri7r 1-z W ; 

(99) 

which was proved by Nenciu in [3A\, implies that 

P-P''^^')xm\\cinMH-)) = ^(^')- (100) 



Step 10: There are Ui,U2 G C{T-C) fl C{V{He)) with norms hounded inde- 
pendently of e satisfying PqUiPq = and U2P0 = P0U2P0 = P0U2 such that 
U = 1 + eUi + e'^U2- In addition, \\PqUi{uY\\c{v{H]^)) ^ 1 for all I G Nq and 
mG {0,1}. 
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We define 

U^ := e-'{Po{U - 1)(1 - Pq) + (1 - Po){U - l)Po) 

and 

f/2 := e-^{Po{U - l)Po + (1 - Po){U - 1)(1 - Pq)). 

Then U = l+eUi+e'^Ui, PqUiPq = 0, and Pof/a = P0U2P0 = U2P0 are clear. 
Next we fix m G No and prove that Ui G C{T>{H^)) with norm bounded 
independent of e. The proof for U2 is similar and will be omitted. We recall 
that 

ij = {PoP + (i-Po)(i-p))5 

with 5 := (1 - (P - Pofy^^^ and that we showed [P, 5] = = [Pq, S] in 
Step O Therefore 

f/i = £-i(Pof/(l-Po) + (l-Po)f/Po) 

= £-l^(PoP(l - Po) + (1 - Po)(l - P)Po) 

= £-i^(Po(P-Po)(l-Po) - (l-Po)(P-Po)Po). (101) 
By Taylor expansion it holds 

1-^= ri(l-s)(l-s(P-Po)2)-^ds(P-Po)2. (102) 
Jo 

Let h{x) := (1 — sa;^) with s G [0, 1]. is holomorphic in Pi/2(0). Due 
to Step [H] the spectrum of P — Pq as an operator on C{V{H^)) is contained 
in Pi/4(0) for e small enough. Therefore \\Rp-Pf^{z)\\c(v{H^)) < 1 for 2; G 
9Pi/2(0) and h{P — Pq) = §qb^ ^(^0) f^{^)Pp-Po{^) dz. This allows us to 

conclude that the integral on the right hand side of fll02p is an operator 
bounded independent of e on T>{H^). This implies that the whole right 
hand side is of order e"^ in C{V{H^)) because ||(P - Pof\\c{viH^)) = 0{e^) 
by Step [61 So we get 

f/i = £-i(Po(P-Po)(l-Po) - (l-Po)(P-Po)Po) + O(£).(103) 

This yields the desired bound because \\P — Po\\ c{v{hij^)) = 0{e). We now turn 
to the claim that \\PoUi{uy\\civ{Hp)) < 1 for m G {0,1}: Using [S,Po] = 
and \\Po{i^y\\c{v(H:^)) ^ 1 due to Step [1] we obtain from fllOll) that 

WPoUiiuYWcivm-)) = \\e''SPo{P-Po){l-Po){iyy\\c(viHr)) 
< ||5-i(P-Po)(//)'|Up(i/r)) 
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We note that the decay properties of P and Pq themselves are not enough. 
Because of the we really need to consider the difference. However, it 
holds P-Po = {P- P^E+') + (pxE+i _ and via ([99]) the first difference 
can be expressed by (P^^^+i)^ — pxE+i_ Looking at the proof of Step [3] we 
see that both differences consist only of terms that carry an e with them and 
have the desired decay property. 

Step 11: For B := PqUx{,H,) and every u G {1, f/^ , f/l) 

II [-e2Ah + ^f,MPo]5||^(^) = 0{e). 

Again we restrict ourselves to the case u = f/j*. It is obvious from the 
definition of Ui in Step [H that [E{,U^Po] = 0. In view of ffTOal) . Ui (and 
thus also UI) contains, up to terms of order e, a factor P — Pq . As long as 
we commute (— £^Ah)Po with the other factors, P — Pq cancels the in the 
definition of Ui and the commutation yields the desired e by Step [H Using 
that B = PqUx = PoX + 0{e) we have 

[-£2 Ah, U:Po]B = [-^^Ah, U:Po]PoX + 0{e) 

^ [-e'A,,,e-\l-Po){P-Po)Po]PoX + 0{e) 
= (1 - Po)[-6^A^,e-\P - Po)]PoXEX + Oie) 
= {1 - Po)[-s^A^,e-\P - Po)xe]PoX + 0{e), 

The last step follows from [{—e'^Ah)PQ,XE]x = ^'(e:), which is implied by 
Lemma m because (— e^Ah)Po satisfies the assumption on A in Lemma El and 
thus 

[H,, i-e'A^)Po] X = [s'A^ + H,, (-e^A^Po] X + Oie) 

= [Vo,-e^A^]Pox - s^A4-e^A^,Po]x + Oie) 
= 0{e) 

as in Step [TJ Furthermore, due to Step [H] 

(1 - Po)[-^'Ah, e-\P - Po)xe]PoX 

= (1 - Po)[-e'A^,E-\P^-+^ - Po)xe]PoX + Oie') 

= (1 - Po)[-£%„ {PiXE+i + XE+iPiil - Xe+i))xe]PoX + 0{e) 

= (1 - Po)[-£'Ah, (T*Po + PoT,)xe]PoX + 0{e). 

On the one hand, 

{1-Po)[-6'A^,PoT,xe] = {l-Po)[-e'A^,Po]PoT,XE = 0{e) 
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by Step d] and Step [2J On the other hand, 

(l-Po)[-£'Ah,r*PoXE]PoX = (l-Po)Ti*Po[(-£'Ah),XE]PoX 

+ (l-Po)[-£'Ah,T*Po]Xi?^oX 

= (1 - Po)T*Po[(-£2Ah)Po, Xe\x + (^(^) 
+ (l-Po)[-£2Ah,r*Po]xi^PoX 

= 0{e) 

due to Step m and the above argument that [(— £:^Ah)Po, = C'(e). 
i^.- ||(f/i* + T*Po) = 0(e) for dime Nq. 

All the following estimates will be in the norm of £(7i, V^H^)) . It is easy to 
prove [Pq, Xe\x = C^(^) the same way we proved [(— £^Ah)Po, Xe]x = ^{^) 
in Step [TUJ Using again that B = PqUx = PoX + ^'(e), x = XeX: well as 
P — Po = 0(e) we obtain that 

f/*P = U*PoXEX + 0{e) 

^ e-\l - Po)(P - Po)PoXEX + Oie) 
= e-\l - Po)(P - Po)xePoX + 0{e) 

^ e-\l - Po)(P^-+^ - Po)xePoX + 0{e) 

= (1 - Po){PiXE+i + (1 - Xe+i)PiXe+i)xePoX + 0(6) 

= (1 - Po)(T*Po + Pon)xEPoX + 0{e) 

= (1 - Po)T*PoX + 0{e) 

= T*PoB + Oie) 

because (1 - Po)T*Po = T*Po by definition and PqX = B + 0{e). 

Step 13: It holds Ef e C^{C). If ipf e C^{C,nf), then 

\\Uo\\c{V{H,),V{-e^^C+Ei)) < 1, \\Uo\\c{V{~e2Ac+Ei)MHe)) ~ 1) 

and there is Xq > 1 with sup^ ||e'^°^''Vf(Q')llwt(g) ~ ^'^^ 

sup||e^"<'^>V:^,...,,,V:^^,...,,^¥.f(g)||^,(,)<l 
<? 

/or a// i/i, . . . , z/i G rb(A^C) and Ti, . . . , G rb(TC). 

We recall that Uqi/j = {{pflip)^^^ and UqiP = V'fV'- Using Lemma [TT] ii) we 
easily obtain \\i-e^Ac + E{)Uo^ < \\e-^°^''^/\V'^f^ for all ^ G V{He) 
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and \\H^U^ij\\ < \\e^Vdij\\ for all G V{-e^Ac + ^f). By dHSD the former 
estimate implies \\UQ\\c(v{He),v{-e^Ac+Ef)) ^ 1- Due to the bounded geom- 
etry of C any differential operator of second order with coefficients in 
is operator-bounded by the elliptic —Ac. So the latter estimate implies 
\\UQ\\c{v{~-e'^Ac+E{),v{Hs)) ^ 1- The other statements are true by Lemma [TTli) 
and ii). 

The results of Step [T] and Steps E] to US] together form Proposition HI □ 
Proof of Lemma \TT[ 

Because of Vq G C^(C, C^{NgC)) and [V^!, A^] = for all r due to Lemmag] 
the mapping q {Hf{q) — z)~^ is in C^{C, C{7ii)). Since Ei is a constraint 
energy band and thus separated, the projection Po(q') associated with -Ef(g) 
is given via the Riesz formula: 



where 7(g) is positively oriented closed curve encircling Ef[q) once. It can 
be chosen independent of g e C locally because the gap condition is uniform. 



Therefore (i/f(-) - z)-^ e C^{C,C{'Hi)) entails Pq e [C , C{'Hi)) . This 



means in particular that PqH is a smooth subbundle. Therefore locally it is 
spanned by a smooth section ipi of normalized eigenfunctions. By 



we see that also EfPo e C^{C,C{Hi)). Then Ef = tinfi-){EfPo) e C^(C) 



because covariant derivatives commute with taking the trace over smooth 
subbundles and derivatives of EfP^ are trace-class operators. For example 



for all r G Ti,(TC) because Pq and E^Pq are trace-class operators and the 
product of a trace-class operator and a bounded operator is again a trace- 
class operator (see e.g. [36], Theorem VI. 19). The argument that higher 
derivatives of E^Pq are trace-class operators is very similar. 



Next we will prove the statement about invariance of exponential decay under 
the application of (^f) := (1 - Po)(^f - ^f)"^! - ^o)- So let ^ eHfhe 






V,tr(PfPo) 



V.tr((EfPo)Po) 
tr((V!?EfPo)Po + (EfPo)V^!Po) 
tr((V^£;fPo)Po) + tr((PfPo)Vj^Po) < 00 
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arbitrary. The claim is equivalent to showing that there is Aq > such that 
for all A G [— Aq, Aq] 

$ := e^<^>/?Hf(^f)e~^<'^>^ 
satisfies ^Wn ^ ll^llw- The latter immediately follows from 

m-H < ||e"<'^>(^f-^f)e"'^'^^$lk (104) 

because 

< W^Wn + sup||e^<^>Poe-'<^>||£(KK.))ll^llw 

<J6C 

< II^IIh, 

where we used that E{ is a constraint energy band by assumption. We now 
turn to (11041) . We note that by the Cauchy-Schwarz inequality it suffices to 
find a Ao > such that for all A G [— Aq, Aq] 

{^\^)n < |Re($|e^<'^>(/Jf-Ef)e-^<'^>$)^| (105) 

To derive fllOSp we start with the following useful estimate, which is easily 
obtained by commuting H{ — Ef with e"'^^'^^ . 

|Re($|e^<'^>(i7f-Ef)e"^<'^><l>)| = |(<l>|(/7f - Ef)*) - A'($|(|i/|V(i/)')$)| 

> \{mHi-Ei)^)\ - A2($|$). 

Since E^ is assumed to be a constraint energy band and thus separated by a 
gap, we have 

|($|(i7f-Ef)<l>)| = |((l-Po)$|(i^f-^f)(l-Po)$>| 

> Cgap((l-Po)$|(l-Po)$> 
= Cg,p((<l>|$)-($|Po$)). 

Since Aq can be chosen arbitrary small, we are left to show that ($|Po'^') 
is strictly smaller than ($|$) independent of A G [— Ao,Ao]. Since E^ is a 
constraint energy band by assumption, we know that there are Aq > and 
C < oo independent of g G C such that We^''^"^ Po{q)e^''''''^\Hf{q) < C. Hence, 

1 = tTn,i,){PS{q)) = tr^,.(,)(e^«<''>Po(g)e^°<^>e-^o<'^)Po(g)e-^«<'^)) 

< ||e^"<'^>Po(g)e^°<'^>||«,(,)tr^,(,)(e-^«<'^)Poe-^"<^>) 

< Ctr„,(,)(e'^«<'^>Poe-^°<^>). 
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So we have that for any A with A G [— Aq, Aq] 
inf tr«,(,)(e-^<'^>Po(g)e-'<''>) > inf trw,.(,) (e-^«<'^>Po(g)e-^°<'^0 > C-\ 

Since Poe'^'^"^^ = PoRHfe'^''''^'^ = by definition of $, we have 

($|Po$) = ($|(Po-e-^<"> Poe"^<"^)$) 

< ($1$) sup tr«,(,)(Po - e-^<'^>Po(g)e-"<^>) 

g 

< ($1$) (suptr^,(,)(Po) - inftr^,(,)(e-^<'^)Po(g)e-^<^>) 

< (l-C-i)(<l>|$), 
which finishes the proof of fllOSp . 

For i) it remains to show that the derivatives of Pq produce exponential decay. 
By definition Pq satisfies 

= (Hf- E{)Po = -A.Po + VoPo - EiPo. (106) 

Let Ti, ...Tm € Fb (TC) be arbitrary. To show that the derivatives of Pq decay 
exponentially, we consider equations obtained by commutating the operator 
identity fll06p with V^^ Since Av commutes with V*^ by Lemma [H this 
yields the following hierachy of equations: 

{Hi-Ei){V'^^Po) = (V,,Ef-V^^K))Po, 

(i:ff-Pf)(v^^,,^Po) = (v.„.,i?f-v!!^,,,Vo)Po + (v.,i?f-vJ:,yo)(v!iPo) 

+ (V.,Ef-V!^^\/o)(V^,Po), 

and analogous equations for higher and mixed derivatives. Applying the 
reduced resolvent RHf{Ef) to both sides of the first equation we obtain that 

(l-Po)(V^^Po) = RHAEi){KEi-KVo)Po. 
From II e^"^"^ Poe^"^"^ \\^^^^ < 1 we conclude that 

because the derivatives of Vq and E{ are globally bounded and application of 
Rhi{E{) preserves exponential decay as we have shown above. Inductively, 
we obtain that 
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The same arguments yield || e^o^'^^ (V^^ ..^^^V';i^_...^^^Po)(l -^o)e^«<''^ H^^^^ < 1 
when we start with = Po{Hf-Ef). The assumption \\e^"^'''^ Poe^"^"'^ ||^^^^ < 1 

immediately implies || e^o<''>Po(V^i,...,^,V^^_..._^^Po)^oe^°^''^ < 1- These 

three statements together result in 

11^ \^ iyi,...,ui^ Ti,...,T,nJ^ WciH) ~ 

We now turn to ii). So we assume that ipf G C^{C,TC{{q)) for some m G Nq. 
By definition iff satisfies 

= {Hi- Ef)ipi = -A^iff + Voiff - Efifi. (107) 

for all qeC. Because of Vq G C^iC, C^{NgC)) and Ef G this is an el- 

liptic equation with coefficients in C^{C,C^{NgC)) on each fibre. Therefore 
ifi G C^{C,C^{NgC)) follows from y^f G C{|(C, 7^f (g)) and standard ellip- 
tic theory immediately. Due to iff G C^{C,T-C{{q)) we may take horizontal 
derivatives of fll07p . Using that [Av, V^] for all r by Lemma [8] ii), we end up 
with the following equations 

{Hi - Ei)V'^^^i = (V.,Ef-VX)<^f, (108) 
(iff - Ef)V^^,,,V.f = (V.„.,Ef-V:^^,,,K,)¥.f + (V.,Ef-V^^Vo)(V^,y.f) 

-F(V.,Ef-V:^,K,)(V!^^^f), 

and analogous equations up to order m. Iteratively, we see that these are all 
elliptic equations with coefficients in C^{C,C^{NgC)) on each fibre. Hence, 
we obtain iff G C™(C, C^{NgC)). So we may take also vertical derivatives of 
the above hierachy: 

{H,-E,)W:^^, = -(v;K,)¥Pf, (109) 

(iff-i?f)v:^v!?^^f = -(v;v^^K,)^f - (v:^yo)(v^^^f) 

+ {Vr,E,-V\Vo)Vl^^i) 
and so on. Since Ef is assumed to be a constraint energy band, we have that 



e 



with a constant independent of q. Choosing ip = e ^'-^^"'^(^f and taking the 
supremum over g G C we obtain the desired exponential decay of ipf. Because 
of Vo G C~(C, C^{NgC)) and E{ G C^{C) also the right-hand sides of (fTUg]) 
and (11091) decay exponentially. By i) an application of RHf{E{) preserves 
exponential decay. So we may conclude that the (/3f-orthogonal parts of Vj^^v^f 

89 



and V^^v'f decay exponentially. Together with the exponential decay of (pf 
this entails the desired exponential decay of Vj^^v^f and V^^v^f . This argument 
can now easily be iterated for the higher derivatives. 

Finally, we turn to iii). We consider a normalized trivializing section iff, in 
particular sup^g^ ll'y^fllwf is globally bounded. By assumption yjf is smooth 
as a section of PoTi.. In order to see that it is also smooth in (1 — Po)TC, 
one applies RHf{Ef) to the equations fllUSI) . which can be justified by an 
approximation argument. Hence, we only need to show boundedness of all 
the derivatives. If C is compact, this is clear. 

We recall that the eigenfunction f{{q) can be chosen real- valued for any q E C. 
If C is contractible, all bundles over C are trivializable. In particular, already 
the real eigenspace bundle PoTi. has a global smooth trivializing section iff. 
We choose a covering of C by geodesic balls of fixed diameter and take an 
arbitrary one of them called fl. We choose geodesic coordinates {x^)i=i,...^d 
and bundle coordinates {n°')a=i,...,k with respect to an orthonormal trivial- 
izing frame {i'a)a over Q as in Remark O Since iff is the only normalized 
element of the real PoTi, we have that 

= f;\ffi, (110) 

for any fixed Xq E Q and x close to it. In view of the coordinate expres- 
sion . = — T'^pu'^dno', we can split up V| _(pi into terms depending 
on .Pq, which are bounded due to i), and — r^n^9„Qy}f (xq). We already 
know that G C^{C,nf{q)). By ii) this implies G C^{C,C^{NgC)) 
with supg ||e'^"^^Vf|| ~ 1- Recalling that {u) = yT+~5^^^n^*n^ we have 
that — r^?T,^(9„c<(y9f (xo) is bounded. Noticing that all the bounds are inde- 
pendent of Q due to ([6j) (as was explained in Remark [6]) we obtain that 
V9f G Cl(C,'H{{q)) . Now we can inductively make use of flllOp and ii) to 
obtain ipi G (C, T^f (g)) . 

If Ef = inf a{Hi{q)) for all q E C, again the real eigenspace bundle is already 
trivializable. To see this we note that the groundstate of a Schrodinger 
operator with a bounded potential can always be chosen strictly positive 
(see [3B]), which defines an orientation on the real eigenspace bundle. A real 
line bundle with an orientation is trivializable. So we may argue as in the 
case of a contractable C that the derivatives are globally bounded. □ 

Proof of Lemma [12]: 

Let the assumption (lOTl) be true for / G No and m G N. The proof for — / G N 
is very similar. We fix zi, . . . , G (C \ M) fl (supp x x [^Ij 1]) and claim 
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that there is a c > independent of the Zi such that 



i=l 



an) 



< 2 



(111) 



To prove this we set $ := UT=iiH - Zi) (Az/)' [17=1 for ^ G 7^ 

and aim to show that ||\1>|| > ||$||/2. We have that 



ll^ll 



i=l 



m m 

j=l i=l 



Using the assumption flOTl) and that \zi\ < 1 for all i we have that there is a 
C < cxo independent of A and the Zi's with 



ll^ll > 



> 

> 
> 



m m 

ccix (jH^^ i[rh{zj) $11 + II n^^^(^^-) "^ii) 

i=i i=i 

m m 

CCiX\\l[HRH^{z^)<t> - CCiX\\1[Rh{z,)<I> 



i=i 

m 



cciAn(i+ 



m 



\lmzj 



||$|| - CCiX Y[\lmzj\-^ ||$| 



CCiX 



i + nr=i(k.i + iim^.i) 



m 



/2 



Now we make use of the Helffer-Sjostrand formula. We recall from the proof 
of Lemma H] that it says that 



1 

IT 



djf{z) RhXz) dz, 



where / is an arbitrary almost analytic extension of /. Here by dz we mean 
again the usual volume measure on C. By assumption x is non-negative. So 
by the spectral theorem we have x{H) = YULi X^^^{H). We choose an almost 
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analytic extension of x^^™ such that K := suppx^''™ C supp x x [—1, 1] (in 
particular the volume of K is finite) and 



Then by the Helffer-Sjostrand formula 



(112) 



X{H) = — Y\d,x'^"^{z,)T\RH{z,)dz,...dzm. 

TT .lirm . . . . 



i=l 



i=l 



We will now combine fill II) and flll2p to obtain the claimed estimate. In the 
following, we use < for 'bounded by a constant independent of H\ 



^ -^C-i^l i=l 
jj^ r. rn m 

< C[ / ni^^^^l WX{RH{z{){u)-'MdZi...dz^ 

Jk-^ i=i i=i ' 



ITT2t 



where we used that (u) (Xu) < X ~ C{ HI^i fo^' small |Imzj|. So 



{uYxmii.)-^^ 



m 



n lib I ft 

< C[ / l[\lmz,\\{Xuyl[RH{zi){u)-^^ 

•^^^ ^=l ' i=l 

„ m m m 

= ci / n iim^,i I n^/^(-2«) n^^ - '^*) 



dz\ . . . dz^^ 



dz\ . . . dz^, 



X l[RH{z,){Xu)~^{Xiyy{u)-'^ 

4 = 1 

„ m m 

■^^^ i=l ^=l 

m m 



i=l 



i=l 



C{H) 



dz\ . . . dZn 



< C'lllv^ll^, 

because of the resolvent estimate fl4UI) and (Az/)'(i/)^' < 1 for A < 1. Hence, 
||(z/)'x(if) {i^)~''\\c{H,v{Hp)) is bounded by C[ times a constant independent 
ofif. ' ' □ 
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Appendix 



Manifolds of bounded geometry 

Here we explain shortly the notion of bounded geometry, which provides the 
natural framework for this work. More on the subject can be found in [H]. 

Definition 6 Let {Ai,g) be a Riemannian manifold and let Vg denote the 
injectivity radius at q E A4. Set tm '■= infggA^rg. {M.,g) is said to be of 
bounded geometry, ifr^ > and every covariant derivative of the Riemann 
tensor 71 is bounded, i.e. 

VmGN 3C„<oo: ^(V'"7^, V'"7^) < Cm- (113) 

Here V is the Levi-Civita connection on {Ai,g) and g is extended to the 
tensor bundles T^^M. for all l,m E N in the canonical way. An open subset 
U d M. equipped with the induced metric g\u is called a subset of bounded 
geometry, if tm > and U13\) is satisfied on U. 

The definition of the Riemann tensor is given below. We note that rj^ > 
implies completeness of M.. The second condition is equivalent to postulating 
that every transition function between an arbitrary pair of geodesic coordi- 
nate charts has bounded derivatives up to any order. Finally, we note that 
the closure of a subset of bounded geometry is obviously metrically complete. 

The geometry of submanifolds 

We recall here some standard concepts from Riemannian geometry. For fur- 
ther information see e.g. [25j. 

First we give the definitions of the inner curvature tensors we use because 
they vary in the literature. We note that they contain statements about 
tensoriality and independence of basis that are not proved here! In the fol- 
lowing, we denote by T{£) the set of all smooth sections of a bundle £ and 
by Tjjj(A^) the set of all smooth (Z,m)-tensor fields over a manifold Ad. 

Definition 7 Let {A,g) be a Riemannian manifold with Levi-Civita connec- 
tion V. Let Ti, T2, Ts, r4 G r(r^). 

i) The curvature mapping R : T{TA) x T(TA) — Tl{A) is given by 
a) The Riemann tensor TZ G T^{A) is given by 

T^iTi,T2,T3,n) := g{Ti,R{T3,n)T2). 
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in) The Ricci tensor Ric e T2i^) given by 



Ric(Ti, T2) tr^ R( . , Ti)t2. 



T/ie scalar curvature k : A^M. is given by 



K :— tr^Ric. 



Here tr_4 1 means contracting the tensor t at any point q & A by an arbitrary 
orthonormal basis ofTgA. 

Remark 7 The dependence on vector fields of R, 71, and Ric can be lifted to 
the cotangent bundle TC* via the metric g. The resulting objects are denoted 
by the same letters throughout this work. The same holds for all the objects 
defined below. 

Of course, all these objects can also be defined for a submanifold once a 
connection has been chosen. There is a canonical choice given by the induced 
connection. 

Definition 8 Let C G A be a submanifold with induced metric g. Denote by 
TC and NC the tangent and the normal bundle of C. Let ti,T2,T3 G r(TC). 

i) We define V to be the induced connection on C given via 



where ti,T2 are canonically lifted to TA = TC x NC and Pt denotes the 
projection onto the first component of the decomposition. The projection 
onto the second component of the decomposition will be denoted by P± . 

a) R, Ric, and k are defined analogously with R, Ric and k from the preceding 
definition. The partial trace of R with respect to C is given by 



Here tic t means contracting the tensor t at any point q ^ C by an arbitrary 
orthonormal basis ofTqC. 

We note that V coincides with the Levi-Civita connection associated with the 
induced metric g. Now we turn to the basic objects related to the embedding 
of a submanifold of arbitrary codimension. 



VriT2 := -PtVtiT2, 



trcR := trcRicc, 
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Definition 9 Let r, Ti, T2 G r{TC), v e V{NC). 

i) The Weingarten mapping W : r{NC) — > l'J(C) is given by 

>V(i/)T := -PrVriy. 

a) The second fundamental form II( . ) : r{NC) — > T2{C) is defined by 

II(i/)(ri,r2) := ^(VnTs,!/). 

Hi) The mean curvature normal 77 e r{NC) is defined to be the unique vector 
field that satisfies 

g{7], v) = tTcW{u) V 1/ e r{NC). 

iv) We define the normal connection V"*" to be the bundle connection on the 
normal bundle given via 

V> := PxV^z/, 

where v and r are canonically lifted to TA = TC x NC. 

v) : r(TC) X r(TC) x T{NC) T{NC) denotes the normal curvature 
mapping defined by 

R\n,T2)i^ := V^^V^i. - ViV> - Vf,^,,,]!.. 

Remark 8 i) The usual relations and symmetry properties for W and II 
also hold for codimension greater than one: 

ll{i^){n,r2) = g{n,W{u)T2) = g{T2,W{u)n) =ll{u){T2,n). 

a) A direct consequence of the definitions is the Weingarten equation: 

V> = Vriy + W(i/)r. 

Hi) The normal curvature mapping K-^ is identically zero, when the dimension 
or the codimension of C is smaller than two. 



95 



Acknowledgements 



We are grateful to David Krejcirik for providing several references and for 
a careful reading of an earlier version resulting in lots of useful comments. 
We are also grateful to Luca Tenuta for helpful remarks when we began 
with this work. Furthermore, we thank Christian Loeschcke, Frank Loose, 
Christian Lubich, Olaf Post, Hans-Michael Stiepan, and Olaf Wittich for 
inspiring discussions about the topic of this paper. 



References 

D. E. Blair, Riemannian Geometry of Contact and Symplectic Manifolds, 
Progress in Mathematics 203, Birkhauser (2002). 

D. Borisov, G. Cardone, Complete asymptotic expansions for the 
eigenvalues of the Dirichlet Laplacian in thin three-dimensional rods, 
larXiv:0910.3907 lyl [math.AP]. 

D. Borisov, P. Freitas, Asymptotics of Dirichlet eigenvalues and eigen- 
functions of the Laplacian on thin domains in M'^, arXiv:0908.2327v l 
[math.AP]. 

F. Bornemann, Homogenization in time of singularly perturbed mechan- 
ical systems. Lecture Notes in Mathematics 1687, Springer (1998). 

J.-L. Brylinski, Loop Spaces, Characteristc Classes and Ceometric 
Quantization, Progress in Mathematics 107, Birkhauser (1993). 

G. Bouchitte, M. L. Mascarenhas, L. Trabucho, On the curvature and 
torsion effects in one dimensional waveguides ESAIM: Control, Optimi- 
sation and Calculus of Variations, 13, 793-808 (2007). 

J. Briining, S. Yu. Dobrokhotov, V. Nekrasov, T. Ya. Tudorovskiy, 
Quantum dynamics in a thin film. L Propagation of localized pertur- 
bations, Russ. J. Math. Phys. 15, 1-16 (2008). 

M. P. do Carmo, Differential Ceometry of Curves and Surfaces, Prentice- 
Hall (1976). 

R. C. T. da Costa, Constraints in quantum mechanics, Phys. Rev. A 
25, 2893-2900 (1982). 



96 



[10] E. B. Davies, Spectral Theory and Differential Operators, Cambridge 
studies in advanced mathematics 42, Cambridge University Press 
(1995). 

[11] G. F. Deir Antonio, L. Tenuta, Semiclassical analysis of constrained 
quantum systems, J. Phys. A 37, 5605-5624 (2004). 

[12] P. A. M. Dirac, Lectures on Quantum Mechanics, Yeshiva Press (1964). 

[13] P. Duclos, P. Exner, Curvature-induced bound states in quantum waveg- 
uides in two and three dimensions. Rev. Math. Phys. 7, 73-102 (1995). 

[14] N. Dunford, J. T. Schwartz, Linear operators part L general theory. Pure 
and apphed mathematics 7, Interscience pubhshers, inc. (1957). 

[15] P. Freitas, D. Krejcirik, Location of the Nodal Set for Thin Curved Tubes, 
Indiana Univ. Math. J. 57, 343-375 (2008). 

[16] L. Friedlander, M. Solomyak, On the spectrum of the Dirichlet Laplacian 
in a narrow infinite strip, in Spectral theory of differential operators. 
M. Sh. Birman 80th anniversary collection, AMS Translations, Series 2, 
Advances in the Mathematical Sciences 225 (2008). 

[17] R. Frocsc, I. Hcrbst, Realizing Holonomic Constraints in Classical and 
Quantum Mechanics, Commun. Math. Phys. 220, 489-535 (2001). 

[18] D. Grieser, Thin tubes in mathematical physics, global analysis and spec- 
tral geometry, in Analysis on Graphs and its Applications, Proc. Sympos. 
Pure Math. 77, Amer. Math. Soc, 565-593 (2008). 

[19] E. Hairer, Ch. Lubich, G. Wanner, Geometric numerical integration. 
Structure-preserving algorithms for ordinary differential equations. Sec- 
ond edition. Springer (2006). 

[20] B. Helffer, J. Sjostrand, Multiple wells in the semi-classical limit I, Com- 
mun. Part. Diff. Eq. 9, 337-408 (1984). 

[21] P. D. Hislop, Exponential decay of two-body eigenf unctions: A review, 
Elec. J. Diff. Eq. 04, 265-288 (2000). 

[22] M. Jaaskelainen, S. Stenholm, Localization in splitting of matter waves, 
Phys. Rev. A 68, 033607 (2003). 

[23] H. Jensen, H. Koppe, Quantum mechanics with constraints, Ann. Phys. 
63, 586-591 (1971). 



97 



[24] D. Krejcifik, Twisting versus bending in quantum wave guides, in Anal- 
ysis on Graphs and its Applications, Proc. Sympos. Pure Math. 77, 
Amer. Math. Soc, 617-636 (2008), see larXiv:0712.3371i ^2 [math-ph] for 
a corrected version. 

S. Lang, Fundamentals of Differential Geometry, Graduate Texts in 
Mathematics 191, Springer, (1999). 

M.-C. Liu, Affine maps of tangent bundles with Sasaki metric. Tensor 
(N.S.) 28, 34 -42 (1974). 

P. Maraner, A complete perturbative expansion for quantum mechanics 
with constraints, J. Phys. A 28, 2939-2951 (1995). 

P. Maraner, Monopole Gauge Fields and Quantum potentials Induced 
by the Geometry in Simple Dynamical Systems, Annals of Physics 246, 
325-346 (1996). 

R. A. Marcus, On the Analytical Mechanics of Ghemical Reactions. 
Quantum Mechanics of Linear Gollisions, J. Chem. Phys. 45, 4493- 
4499 (1966). 

A. Martinez, V. Sordoni, On the Time-Dependent Born-Oppenheimer 
Approximation with Smooth Potential, Comptes Rendus Acad. Sci. Paris 
337, 185 -188 (2002). 

A. Martinez, V. Sordoni, Twisted pseudodifferential calculus and appli- 
cation to the quantum evolution of molecules, Mem. Amer. Math. Soc. 
200, no. 936 (2009). 

K. A. Mitchell, Gauge fields and extrapotentials in constrained quantum 
systems, Phys. Rev. A 63, 042112 (2001). 

G. Nenciu, Linear Adiabatic Theory, Exponential Estimates, Comm. 
Math. Phys. 152, 479-496 (1993). 

G. Nenciu, V. Sordoni, Semiclassical limit for multistate Klein- Gordon 
systems: almost invariant subspaces and scattering theory, J. Math. 
Phys. 45, 3676-3696 (2004). 

G. Panati, H. Spohn, S. Teufel, The time- dependent Born-Oppenheimer 
approximation. Math. Modelling and Num. Anal. 41, 297-314 (2007). 

M. Reed, B. Simon, Methods of Modern Mathematical Physics L Func- 
tional Analysis, Academic Press (1978). 



98 



[37] M. Reed, B. Simon, Methods of Modern Mathematical Physics II: 
Fourier Analysis, Self-Adjointness, Academic Press (1978). 

[38] M. Reed, B. Simon, Methods of Modern Mathematical Physics IV: Anal- 
ysis of Operators, Academic Press (1978). 

[39] H. Rubin, P. Ungar, Motion under a strong constraining force, Commun. 
Pure Appl. Math. 10, 28-42 (1957). 

[40] T. Sakai, Riemannian Geometry, Translations of Mathematical Mono- 
graphs 149, AMS (1997). 

[41] M. A. Shubin, Spectral theory of elliptic operators on non-compact man- 
ifolds, Asterisque 207, 35-108 (1992). 

[42] V. Sordoni, Reduction Scheme for Semiclassical Operator-Valued 
Schrddinger Type Equation and Application to Scattering, Commun. 
Part. Diff. Eq. 28, 1221-1236 (2003). 

[43] L. Tenuta, S. Teufel, Effective dynamics for particles coupled to a quan- 
tized scalar field, Commun. Math. Phys. 280, 751-805 (2008). 

[44] S. Teufel, Adiabatic Perturbation Theory in Quantum Dynamics, Lecture 
Notes in Mathematics 1821, Springer (2003). 

[45] O. Wittich, L"^ -Homogenization of Heat Equations on Tubular Neighbor- 
hoods. larXiv:0810.5047k Al [math.AP]. 



99 



